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Several  studies  have  con^sared  performance  in  memory  search  tasks  when 
the  memorized  items  are  either  in  primary  or  secondary  memory  (see  footnote 
1).  The  focus  of  these  studies  is  on  the  dynamics  of  the  process  of 
retrieval  of  information  from  secondary  memory  to  primary  memory.  Using 
additive  factors  (Sternberg,  1969b)  and  subtraction  (Bonders,  1868/1969) 
logic,  one  can  make  inferences  about  the  processes  which  are  involved  in  the 
transfer  of  information  and  the  relative  duration  of  these  processes.  This 
can  be  achieved  by  comparing  the  regression  equation  parameters  of 
performance  in  primary  and  secondary  memory.  Differences  in  the  intercepts 
and  slopes  of  primary  and  secondary  menory  reveal  the  characteristics  of 
this  transfer  process.  Typically  these  experiments  employ  a  distractor  task 
between  the  presentation  of  the  manory  set  and  the  probe  stimulus  to  prevent 
rehearsal  of  the  memory  set  items.  Without  rehearsal,  information  in 
primary  memory  decays  rapidly  (e.g.,  Murdock,  1961;  Muter,  1980;  Peterson  & 
Peterson,  1959) ,  thus  the  information  must  be  retrieved  from  secondary 
memory  for  subjects  to  perform  the  memory  search  task. 

Sternberg  (Sternberg,  1969a,  exp.  5;  Sternberg,  Kroll,  and  Nasto,  1969) 
was  one  of  the  first  to  compare  primary  and  secondary  memory  performance  in 
the  memory  search  task.  Subjects  memorized  a  list  of  1,  3,  or  5  digits.  In 
the  secondary  memory  condition,  a  list  of  7  letters  was  sequentially 
presented  (with  a  total  duration  of  3.5  seconds)  prior  to  the  presentation 
of  tiie  probe  stimulus.  Subjects  were  instructed  to  retain  the  list  of 
letters  in  memory,  and  catch  trials  were  included  to  insure  that  subjects 
complied  with  instructions.  The  retention  of  the  list  of  letters  was 
intended  to  prevent  subjects  from  rehearsing  the  memory  set  and  require  the 
retrieval  of  the  memory  set  fran  secondary  memory.  In  the  primary  manory 
condition,  no  distractor  task  preceded  the  probe  stimulus.  Ihe  linear 
regression  equation  for  primary  memory  was  PT  =  336  +  57  (x),  while  for 
secondary  memory  RT  =  467  +  105  (x),  vrfiere  x  refers  to  memory  set  size.  Both 
the  intercept  and  the  slope  were  greater  for  the  secondary  memory  condition. 
The  increase  in  the  intercept  (131  msec)  was  taken  to  reflect  the  time  to 
locate  the  memory  set  in  secondary  memory.  The  increase  in  the  slope  (48 
msec)  was  taken  to  suggest  a  serial  transfer  of  the  entire  memory  set  into 
primary  memory.  One  caveat  in  interpreting  the  differences  between  primary 
and  secondary  memory  in  these  studies  concerns  the  duration  of  the 
distractor  task.  Evidence  frcxti  the  Brown-Peterson  paradigm  suggests  that 
the  memory  set  information  may  not  have  completely  decayed  from  primary 
menory  within  the  3.5  second  interval  (see  Flexser,  1978  for  a  similar 
argument).  If  this  were  the  case,  the  increase  in  the  slope  from  primary 
memory  to  secondary  memory  may  be  the  result  of  a  mixture  of  primary  and 
secondary  memory  performance,  and  thus  the  slope  differences  are  airbiguous. 
Furthermore,  if  the  memory  set  information  had  not  fully  decayed  from 
primary  memory,  the  difference  between  primary  and  secondary  memory  would 
underestimate  the  true  retrieval  time. 

Wickens  and  colleagues  (Coyne,  Allen,  and  Wickens,  1986;  Wickens, 

Moody,  and  Dow,  1981;  Wickens,  Moody,  and  Vidulich,  1985)  have  provided  the 
most  extensive  studies  comparing  primary  and  secondary  memory.  Wickens  et 
al.  (1981,  exp.  1)  presented  memory  sets  of  2  or  4  words  for  3  seconds. 
Following  the  memory  set,  subjects  were  presented  with  a  random  three  digit 
number  and  instructed  to  count  backwards  by  3s  for  a  duration  of  12  seconds. 
A  2  second  interval  in  which  the  random  number  was  removed  fran  the  display 
signaled  the  end  of  the  distractor  task  and  the  upcoming  probe  stimulus.  No 
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measure  of  subjects'  performance  was  obtained  for  the  distractor  task.  Ihe 
probe  was  presented  for  2  seconds,  and  the  interval  in  which  subjects  could 
respond  was  3  seconds  from  probe  onset.  A  manipulation  investigating  the 
effects  of  proactive  interference  on  primary  and  secondary  menory  was 
accomplished  by  selecting  the  words  in  the  memory  set  from  the  same 
taxonomic  category.  Three  trials  using  stimuli  from  the  same  taxonomic 
category  were  presented  (without  repetition  of  words)  before  a  new  category 
was  chosen.  The  first  memory  set  was  deemed  low  in  interference,  the  third 
high  in  interference.  Experimental  factors  were  mixed  within  the  session. 
Wickens  et  al.  (1981,  exp.  2)  was  identical  to  Wickens  et  al.  (1981,  exp. 

1),  with  the  exception  that  the  distractor  task  was  eliminated.  This 
constituted  the  primary  memory  condition.  Note  that  the  primary  -  secondary 
memory  manipulation  was  a  between  subjects  factor.  Thus  coirparisons  between 
experiments  1  and  2  provide  a  measure  of  the  time  to  transfer  the  memory  set 
from  secondary  memory  to  primary  memory. 

The  slopes  and  intercepts  obtained  in  the  experiment  were: 


Primary  Secondary 

target  RT  =  500  +  33 (x)  RT  =  595  +  39 (x) 

non-target  RT  =  501  +  36.25(x)  RT  =  630  +  37.5 (x) 


It  is  evident  that  manory  load  did  not  interact  with  the  manipulation  of 
primary  and  secondary  memory,  vAiich  is  in  contrast  to  Sternberg's  (1969a, 
exp.  5)  results.  (Dne  potential  reason  for  the  discrepancy  is  the  nature  and 
duration  of  the  distractor  task.  The  Wickens  study  prevents  subjects  from 
rehearsing  the  memory  set  for  a  longer  duration  than  Sternberg's  studies, 
and  probably  reflects  a  more  accurate  representation  of  secondary  memory 
performance.  It  is  also  apparent  that  there  is  a  substantial  difference  in 
the  intercepts  of  primary  and  secondary  memory  (112  msec).  This  estimate  is 
more  closely  in  agreement  with  the  Sternberg  studies. 

The  data  from  experiments  1  and  2  were  fitted  to  the  following 
regression  equation; 

RT  =  [  a  +  b(x)  +  t  I  +  [(  rl  +  r2)  +  q  ] 


The  first  bracketed  term  in  the  equation  describes  performance  in  primary 
memory,  where  "  a  "  refers  to  the  intercept,  "  b  "  refers  to  the  slope, 

"  X  "  refers  to  memory  set  size,  and  "  t"  refers  to  differences  between 
target  and  non- target  stimuli.  The  second  bracketed  term  describes 
performance  in  secondary  memory.  The  term  "  rl  "  refers  to  the  increase  in 
the  intercept  under  low  proactive  interference  conditions,  "  r2  "  refers  to 
the  additional  increase  under  high  proactive  interference  conditions,  and 
"  q  "  refers  to  differential  effects  for  target  and  non- target  stimuli.  The 
values  for  the  terms  in  the  equation  are; 

RT  =  [  489  +  37(x)  +  14  ]  +  [(  95  +  38)  +  17  ] 
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This  model  accounted  for  99%  of  the  variance  across  conditions.  It  is 
iii5»rtant  to  note  that  these  results  iirply  that  the  retrieval  process  is 
independent  (but  see  Sternberg,  1969b)  of  the  memory  search  process. 
According  to  this  interpretation,  there  is  a  constant  time  for  retrieving 
the  memory  set  items,  and  they  are  retrieved  in  a  chunk  (i.e.,  in  parallel). 
Furthermore,  the  effects  of  memory  load  produced  equivalent  results, 
suggesting  that  the  monory  search  process  is  the  same  for  items  in  primary 
memory  as  for  items  which  have  to  be  transferred  from  secondary  memory  to 
primary  memory. 

Wickens  (Wickens  et.  al.,  1985,  exp.  1;  Wickens,  personal  comnunication, 
March  12,  1987)  has  also  demonstrated  the  independence  of  memory  search  rate 
and  memory  retrieval  time  using  both  out  of  category  non-target  probes  and 
repeated  stimuli.  Both  effects  reduced  the  slope,  but  the  primary  - 
secondary  memory  differences  were  maintained.  Flexser  (1978)  has  also 
reported  similar  intercept  differences  between  primary  and  secondary  memory 
using  memory  sets  of  16  and  32  (which  produced  small  effects  on  search 
time).  These  results  offer  additional  support  for  the  assertion  that  the 
retrieval  process  and  the  memory  search  process  are  independent.  Wickens 
concluded  that  secondary  memory  conditions  differ  frcxn  primary  memory 
conditions  in  the  addition  of  a  retrieval  process  inserted  prior  to  the 
memory  search  operation. 

In  sum,  the  studies  which  have  coirpared  primary  and  secondary  memory 
using  memory  search  paradigms  have  universally  found  that  there  is  an 
increase  in  the  intercept  between  primary  and  secondary  memory.  Hie 
magnitude  of  this  effect  seems  to  vary  with  the  type  of  material  —  words 
produce  a  larger  effect  than  letters  or  digits  (e.g.,  Wickais  et  al.,  1985, 
exp.  2).  Furthermore,  there  is  some  evidence  that  the  magnitude  of  the 
difference  between  primary  and  secondary  memory  is  dependent  on  the  strength 
of  the  trace  in  primary  memory.  The  longer  the  interval  between  the  memory 
set  and  the  probe  stimulus,  the  greater  the  difference  between  primary  and 
secondary  memory.  This  difference  in  the  intercept  between  primary  and 
secondary  memory  has  been  interpreted  as  the  time  to  retrieve  the 
infonnation  from  secondary  memory. 

The  effect  of  memory  load  on  primary  and  secondary  memory  has  produced 
mixed  results.  The  majority  of  studies  report  parallel  curves  for  primary 
and  secondary  memory  (Wickens  et  al.,  1981;  Wickens  et  al.,  1985;  Wickens, 
personal  conmunication,  1987;  Flexser,  1978).  Forrin  and  Morin  (1969) 
reported  shallower  slopes  for  the  secondary  memory  condition,  but  subjects 
received  substantially  more  practice  (i.e.,  consistent  practice)  in  the 
secondary  memory  condition  than  in  the  primary  memory  condition.  A  few 
studies  have  reported  steeper  slopes  for  secondary  memory  conditions  (e.g., 
Sternberg,  1969;  Sternberg,  Kroll,  and  Nasto,  1969;  Peters,  1974).  However, 
in  these  studies  the  interval  between  the  memory  set  and  the  probe  stimulus 
was  such  that  it  is  questionable  whether  the  memory  set  information  had 
decc^ed  from  primary  memory.  If  one  assumes  that  rehearsal  strength's  the 
trace  in  primary  memory,  as  set  size  increases  the  average  trace  str^igth  of 
each  memory  set  item  should  decrease.  This  would  result  in  differential 
decay  rates  for  different  memory  set  sizes.  Sknaller  memory  set  sizes  should 
take  longer  to  decay,  since  the  rehearsal  cycle  time  varies  with  memory  load 
(e.g.,  Ba^eley  and  Ecob,  1973;  Cavanaugh,  1972;  Clifton  and  Bireibaum, 

1970;  Corballis,  Kirby,  and  Miller,  1972;  Monsell,  1978).  Hiis  would 
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prcxJuce  a  greater  proportion  of  primary  memory  trials  following  a  short 
distractor  task  for  smaller  set  sizes  than  for  larger  memory  set  sizes.  Ihe 
expected  outcome  would  result  in  a  steeper  slope  for  this  mixture  condition 
than  for  a  pure  primary  memory  or  a  pure  secondary  memory  condition. 

An  important  consideration  when  evaluating  the  differences  between  the 
studies  concerns  the  nature  of  the  distractor  task,  ibe  distractor  tasks 
varied  in  difficulty  and  duration,  and  none  of  the  experiments  objectively 
evaluated  subjects'  performance  in  the  distractor  task.  The  most  effective 
distractor  task  employed  a  backwards  count  task  for  a  duration  of  at  least 
12  seconds.  All  of  these  studies  reported  parallel  manory  load  functions 
for  primary  and  secondary  memory.  However,  since  performance  in  the 
distractor  task  was  not  monitored,  it  is  possible  that  s'lbjects  may  have 
rehearsed  the  memory  set  even  in  these  conditions.  Clearly  a  distractor 
task  in  which  performance  can  be  monitored  over  the  course  of  the  distractor 
interval  is  warranted. 

A  final  issue  concerns  the  interval  between  when  the  distractor  task 
was  terminated  and  when  the  probe  stimulus  was  presented.  This  duration 
varied  from  500  msec  to  2  seconds  in  the  studies  reviewed  above.  It  is 
unclear  what  subjects  were  doing  in  this  interval,  and  why  they  did  not 
transfer  the  memory  set  information  from  secondary  to  primary  memory.  Given 
the  estimates  of  retrieval  time,  subjects  should  have  had  ample  time  to 
complete  this  operation  prior  to  the  presentation  of  the  prciie  stimwlus. 

Upon  post-experimental  questioning,  subjects  reported  that  they  did  not 
retrieve  the  memory  set  information  prior  to  the  presentation  of  the  probe 
stimulus  (Wickens  et  al.,  1985).  If  subjects  had,  in  fact,  retrieved  any 
information  from  secondary  memory,  this  should  tend  to  obscure  any  true 
differences  between  primary  and  secondary  memory. 

Experiment  1;  Primary/ Secondary  Memory  Search 

Several  methodological  and  interpretative  issues  remain  which  need  to 
be  addressed  before  this  paradigm  can  be  used  to  assess  the  dynamics  of 
memory  retrieval.  The  purpose  of  experiment  1  was  to  address  these  issues 
and  provide  additional  information  concerning  the  process  of  retrieval  of 
information  fran  secondary  menory  to  primary  memory. 

One  issue  to  be  addressed  is  the  proposition  that  the  retrieval  process 
is  independent  from  the  monory  search  process.  There  was  a  trend  in  the 
literature  for  studies  using  short  distractor  intervals  (e.g.,  4  seconds)  to 
find  larger  slopes  for  secondary  memory  conditions,  while  studies  which  used 
longer  distractor  intervals  (e.g.,  12  seconds)  found  equivalent  slopes  for 
primary  and  secondary  memory.  The  present  experiment  will  evaluate  the 
contribution  of  the  duration  of  the  distractor  interval  on  the  additivity  of 
the  retrieval  and  memory  corparison  processes  by  manipulating  the  delay 
between  the  memory  set  and  the  presentation  of  the  probe  stimulus.  Three 
delay  conditions  will  be  included:  0,  4,  or  15  sec.  The  0  delay  condition 
represents  the  situation  in  which  no  distractor  task  is  presented  between 
the  memory  set  and  the  probe  (i.e.,  a  pure  primary  memory  condition).  In 
the  4  second  delay  condition,  subjects  will  perform  a  distractor  task  for  4 
seconds  prior  to  the  presentation  of  the  probe  stimulus.  The  duration  of 
the  distractor  task  in  this  condition  is  quite  close  to  the  parameters  of 
Sternberg's  experiment  in  which  an  increase  in  the  set  size  slope  was 


4 


observed  conpared  to  primary  memory.  In  the  15  second  delay  condition, 
subjects  will  perform  a  distractor  task  for  15  seconds  prior  to  the 
presentation  of  the  probe.  This  condition  is  roughly  equivalent  to  that 
employed  in  the  Wickens  studies,  and  is  intended  to  represent  a  condition  in 
which  the  memory  set  information  is  not  in  primary  memory  (i.e.,  a  pure 
secondary  memory  conditiorO.  Following  additive  factors  logic,  if  the 
memory  retrieval  process  is  independent  of  the  memory  search  process,  then 
the  effects  of  delay  should  be  additive  with  the  effects  of  memory  load. 

A  second  issue  to  be  addressed  by  this  research  is  related  to  the 
finding  in  the  Wickens  study  that  subjects  apparently  did  not  retrieve  the 
memory  set  until  the  probe  was  presented,  even  though  they  were  cued  to  stop 
the  distractor  task  and  prepare  for  the  prcfce  trial  with  a  2  second  warning 
interval.  It  is  unclear  why  subjects  did  not  retrieve  the  memory  set  during 
this  interval.  If  subjects  retrieved  any  information  prior  to  the 
presentation  of  the  probe  stimulus,  the  estimates  of  retrieval  time 
represent  an  underestimate  of  the  true  retrieval  time.  The  present  study 
will  present  a  retrieval  cue  which  signals  the  termination  of  the  distractor 
task  (if  presented).  Hie  prdae  stimulus  will  be  presented  at  three  stimulus 
onset  asynchronies  (SQAs)  following  retrieval  cue  onset:  200,  500,  and  1000 
msec.  If  subjects  are  retrieving  any  of  the  memory  set  information  prior  to 
the  presentation  of  the  probe,  then  the  retrieval  time  estimate  (the 
intercept  difference  between  primary  and  secondary  memory)  should  diminish 
as  SOA  increases.  The  three  SQA  conditions  will  be  factorial ly  combined 
with  the  three  delay  conditions.  IWo  memory  set  sizes  will  be  presented:  2 
and  4  and  target  and  non-target  trials  will  be  equiprobable.  All  conditions 
will  be  mixed  to  reduce  changes  in  bias  between  conditions.  It  is  predicted 
that  reaction  time  will  increase  as  a  function  of  delay  and  that  this 
increase  reflects  the  time  to  retrieve  the  memory  set  information  from 
secondary  memory.  It  is  also  predicted  that  if  the  transfer  of  the  memory 
set  information  is  parallel,  as  suggested  by  Wickens,  then  the  slopes  for 
the  0  and  15  sec  delay  conditions  should  be  equivalent.  If  performance  in 
the  4  sec  delay  condition  reflects  a  mixture  of  primary  and  secondary  memory 
conditions  and  this  proportion  varies  with  memory  load,  then  there  should  be 
an  increase  in  the  slope  for  this  condition  relative  to  the  0  and  15  sec 
delay  conditions.  These  effects  should  be  modulated  by  SQft.  Hie  largest 
differences  between  delay  conditions  should  be  observed  with  the  shortest 
SQA  condition.  If  subjects  retrieve  any  information  prior  to  the 
presentation  of  the  probe  stimulus,  the  difference  between  delay  conditions 
should  decrease  as  SQA  increases. 

Methods:  Experiment  1 


Subjects 

Tai  subjects  (3  female)  age  18  to  25  participated  in  experiment  1. 
Each  subject  participated  in  6  one  hour  sessions.  Subjects  were  paid  for 
their  participation. 

Stimuli 


The  stimuli  presented  in  the  Sternberg  portion  of  the  experiment  were  4 
letter,  monosyllabic  nouns  with  a  word  frequency  of  70  to  148  (Kucera, 

Hairy,  Francis,  and  Nelson,  1967).  Words  were  selected  so  as  not  to  rhyme 
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or  have  orthographic  similarity.  A  further  constraint  used  in  stimulus 
selection  was  to  examine  all  possible  word  pairs  and  elindnate  any  words 
which  resulted  in  strong  natural  associations.  Thus  the  words  formed  a 
heterogeneous  class  of  stimuli.  Appendix  E  presents  the  word-list  employed 
in  the  experiment.  The  approximate  visual  angle  of  the  words  subtended  1.2 
degrees  horizontally  and  0.5  degrees  vertically. 

The  stimuli  used  in  the  recognition  running  memory  portion  of  the 
experiment  were  the  digits  1  to  9.  Digits  were  selected  randomly.  The 
approximate  visual  angle  of  the  running  memory  stimuli  subtended  OJ  degrees 
horizontally  and  0.5  degrees  vertically. 

Apparatus 

The  experiment  was  performed  on  an  IBM  XT,  with  a  quadEGA  card  which 
permitted  cursor  control  and  synchronization.  The  stimuli  were  displayed  on 
an  IBM  monochrome  display.  Subjects  indicated  their  responses  by  pressing 
the  "Z"  key  with  the  left  index  finger  and  the  key  with  the  right  index 
finger  on  the  keyboard  of  the  IBM  computer. 

Procedure 


Each  trial  was  comprised  of  the  following  events.  A  monory  set  was 
presented  for  a  duration  of  3000  msec.  This  was  followed  by  a  1500  msec 
interval  in  which  the  display  was  blanked.  Three  delay  intervals  followed 
the  memory  set:  0,  4,  and  15  seconds.  During  the  delay  interval,  subjects 
performed  a  recognition  running  memory  task  which  prevented  rehearsal  of  the 
memory  set.  Following  the  delay  interval,  an  asterisk  was  presented  for  200 
msec  which  served  as  a  cue  for  subjects  to  retrieve  the  memory  set 
information.  Three  stimulus  onset  asynchronies  (SOAs)  were  included  between 
onset  of  the  retrieval  cue  and  the  onset  of  the  Sternberg  probe  stimulus: 

200,  500,  and  1000  msec.  The  Sternberg  probe  was  presented  for  200  msec  and 
legal  reaction  times  were  permitted  within  a  3000  msec  interval  following 
probe  onset. 

Sternberg  memory  set  sizes  of  2  and  4  words  were  used  in  the 
experiment.  Target  and  non-target  trials  were  presented  equiprobably. 

Targets  were  defined  as  items  from  the  memory  set.  Non- targets  were  items 
not  included  in  the  memory  set.  Subjects  pressed  one  key  for  target  trials 
and  another  key  for  non-target  trials.  Key  assignments  were  counterbalanced 
across  subjects.  Instructions  emphasized  both  speed  and  accuracy. 

The  recognition  running  memory  task  was  performed  in  the  interval 
between  the  presentation  of  the  memory  set  and  the  retrieval  cue.  The  task 
consisted  of  a  series  of  digits  presented  successively  for  200  msec,  with  an 
interstimulus  interval  of  1000  msec.  Ihe  subject's  task  was  to  press  one 
button  if  the  digit  oiesented  on  trial  N  matched  the  digit  presented  on 
trial  N-2  and  another  button  if  it  did  not  match  the  item.  It  is  important 
to  note  that  each  stimulus  served  as  a  probe  in  the  recognition  task,  and 
subsequently  served  as  a  tenplate  against  which  digits  presented  two  trials 
later  must  be  conpered.  For  subjects  to  successfully  perform  the  task,  it 
is  necessary  for  them  to  maintain  the  last  two  digits  in  memory.  Digits 
were  chosen  randomly,  with  the  constraint  that  mismatched  stimuli  occurred 
twice  as  often  as  matched  stimuli.  This  constraint  was  introduced  to  keep  the 
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difficulty  of  the  running  memory  task  constant  across  trials.  Subjects  were 
given  1000  msec  to  indicate  their  response.  Instructions  emphasized  both 
speed  and  accuracy.  It  should  be  re-emphasized  that  subjects  found  the 
running  memory  task  extremely  demanding,  requiring  all  their  effort  to 
maintain  performance  in  the  task.  While  subjects  performed  the  running 
memory  task,  they  were  instructed  not  to  rehearse  or  otherwise  think  about 
the  memory  set  items.  Performance  in  the  running  memory  task  provided  a 
good  index  of  subject's  compliance  to  these  instructions  (see  footnote  2). 

Design 

Memory  set  sizes  of  2  and  4  were  presented,  with  target  and  non-target 
trials  equiprobable.  Three  delay  conditions  were  included:  0,  4,  and  15 
seconds.  During  the  delay,  subjects  performed  the  recognition  running 
memory  task,  which  served  to  prevent  subjects  from  rehearsing  the  memory 
set.  In  addition,  three  SQAs  between  retrieval  cue  onset  and  prdoe  stimulus 
onset  were  employed:  200,  500,  and  1000  msec.  Conditions  were  factorial  ly 
combined  to  form  a  2  (memory  load)  X  2  (target  vs  non-target)  X  3  (delay)  X 
3  (SQA)  design.  Conditions  were  randomly  permuted  across  sessions  and 
subjects.  Subjects  participated  in  all  experimental  conditions.  A  total  of 
36  observations  per  cell  were  obtained  for  each  subject.  Uie  first  session 
of  the  experiment  served  as  practice  and  was  not  included  in  the  analyses. 

Results  and  Discussion;  Experiment  1 

Figure  1  presents  the  mean  reaction  times  obtained  in  the  Sternberg 
task.  Only  trials  in  which  subjects  were  correct  are  included.  An 
additional  constraint  employed  in  selecting  trials  for  analysis  omitted 
trials  in  which  performance  in  the  running  memory  task  (if  presented)  was 
below  chance  accuracy.  Data  are  coll^sed  across  target  and  non-target 
trials  for  the  three  different  SOA  conditions.  Line  segments  connect  memory 
set  size  means  within  a  condition.  Several  effects  are  noteworthy.  First, 
reaction  time  increased  as  a  function  of  rnemoLy  load  in  all  conditions, 
F(l,9)=36.56,  p<.001,  MSeFl7194.  Memory  load  did  not  interact  with  any 
other  variables  (all  ps  >  J.0),  suggesting  that  the  memory  search  process 
did  not  differ  between  tasks.  This  conclusion  is  borne  out  by  the  linear 
regression  slopes  fitted  to  each  of  the  conditions  in  the  experiment.  The 
average  slope  in  the  experiment  was  42.0  msec  per  item.  A  comparison  of  the 
0  and  15  second  delay  conditions  replicates  the  effects  reported  by  Wickens. 
The  l:^pothesis  that  the  difference  in  slopes  between  primary  and  secondary 
memory  reported  by  Sternberg  was  due  to  the  difference  in  the  duration  of 
the  distractor  task  was  rejected.  This  follows  given  the  equivalence  of  the 
slopes  across  all  delay  conditions.  Furthermore,  these  data  are  consistent 
with  a  nodel  in  which  the  retrieval  process  is  independent  of  the  memory 
search  process,  since  memory  load  did  not  interact  with  delay. 


A  second  noteworthy  aspect  of  the  data  is  the  effect  of  SQA  on  reaction 
time.  Ihere  was  a  general  trend  for  reaction  time  to  decrease  as  SQA 
increased,  F(2,18)=62.93,  p<.001,  MSe^SOOS.  We  interpret  this  change  as  a 
non-specific  warning  effect.  Hie  greater  the  time  between  the  onset  of  the 
retrieval  cue  and  the  onset  of  the  probe,  the  more  subjects  prepared  for  the 
upccxning  pr(^e.  This  effect  was  non-monotonic,  with  the  greatest  effect 
obtained  betweei  the  200  and  500  msec  SQA  conditions  (67  msec),  and  a  lesser 
effect  between  the  500  and  1000  msec  SQA  conditions  (42  msec). 
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Memory  Load 


Mean  reaction  time  for  each  condition  in  experiment  1.  The  data 
are  plotted  as  a  function  of  memory  set  size  for  each  delay  and 
SOA  condition.  The  left  column  in  the  figure  represents  the  200 
msec  SOA  condition,  the  center  column  represents  the  500  msec  SOA 
condition,  and  the  right  column  represents  the  1000  msec  SOA 
condition. 


FIGURE  1 
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The  manipulation  of  delay  produced  large  differences  in  performance. 

As  delay  increased,  reaction  time  increased,  F (2 ,18)  =18 .26,  p<.001, 
MSe=91,876.  Iliis  effect  produced  changes  in  the  linear  regression 
intercept,  but  not  the  slope.  Following  Sternberg's  additive  factors  logic, 
the  additivity  suggests  that  the  delay  and  memory  loads  affect  different 
processes,  llie  former  presumably  affects  a  retrieval  process,  while  the 
latter  affects  a  memory  search  process. 

The  effect  of  delay  was  modulated  by  SQA,  F(4,36)=9.81,  p<.001, 
MSe=3594.  As  SQA  increased,  the  effect  of  delay  decreased.  However,  this 
reduction  in  reaction  time  occurred  largely  between  the  0  and  4  second  delay 
conditions,  v^ile  the  difference  between  the  4  and  15  second  delay 
conditions  remained  constant  across  SOA  conditions.  The  difference  between 
the  0  and  15  second  delay  conditions  was  approximately  285  msec  in  the  200 
msec  SQA  condition.  This  decreased  to  approximately  175  msec  in  the  1000 
msec  SQA  condition. 

Hie  difference  in  delay  conditions  has  been  used  to  infer  the  duration 
of  the  retrieval  process.  Since  the  magnitude  of  the  effect  of  delay 
decreases  as  SQA  increases,  it  suggests  that  part  of  the  memory  set 
information  is  retrieved  during  the  SQA  interval.  It  is  unclear  why 
subjects  retrieved  onl^  part  of  the  information  during  the  SQA  interval. 
Given  the  estimate  of  retrieval  time  (275  msec)  and  the  time  to  process  the 
retrieval  cue  (upperbound  estimate  of  300  msec  from  siitple  KT  tasks), 
subjects  should  have  been  able  to  retrieve  the  memory  set  information  within 
600  msec.  Thus  subjects  are  not  fully  capitalizing  on  the  SQA  interval. 

One  reason  that  reaction  time  differed  as  a  function  of  delay  in  the 
1000  msec  SQA  condition  might  be  that  this  condition  represents  a  mixture  of 
primary  and  secondary  memory  trials.  Cii  some  trials,  the  memory  set 
information  may  be  retrieved,  while  in  others  the  retrieval  may  not  have 
been  caipleted.  If  the  proportion  of  primary  memory  trials  to  secondary 
memory  trials  increased  as  SQA  increased,  the  effect  of  delay  should 
diminish  as  SQA  increases,  which  is  precisely  what  is  observed  in  the  data. 
If  this  hypothesis  is  correct,  then  the  1000  msec  SQA  reaction  time 
distribution  should  represent  a  mixture  of  pure  primary  and  pure  secondary 
memory  trials. 

To  address  the  mixture  hypothesis,  the  reaction  time  distributions  for 
each  subject  and  each  condition  were  vincentized  (Vincent,  1912;  Katcliff, 
1979)  to  form  a  composite  cumulative  distribution  function  (GIF).  The 
vincentizing  procedure  is  described  in  detail  by  Ratcliff  (1979).  Briefly, 
the  reaction  times  for  each  subject  and  condition  are  sorted  into  ascending 
order  and  the  quantiles  are  calculated.  Hie  quantiles  are  then  averaged 
across  subjects  to  obtain  the  group  quantiles.  From  the  group  quantiles  a 
group  reaction  time  distribution  is  generated  which  retains  the  shape  of  the 
individual  subject  distributions.  Hiis  process  is  equivalent  to  a  simple 
linear  interpolation,  and  the  resulting  distribution  represents  the 
distribution  of  the  average  subject.  Hie  OFs  have  been  divided  into  20 
intervals,  each  containing  5%  of  the  distribution.  Hie  mean  of  each 
interval  is  cross-plotted  against  the  interval  position. 

The  CDFs  for  the  0  and  15  second  delays  are  plotted  at  each  SOA  for 
target  memory  load  4  trials  in  figure  2.  All  the  OFs  have  a  general 
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scallop  shape  which  reflects  a  positive  skew  in  the  FT  distrioutions.  Of 
central  interest  is  the  shape  of  the  distributions  for  the  different  delays 
as  a  function  of  SQA.  For  the  0  delay  condition,  the  CIFs  are  tightly 
clustered  and  parallel  over  the  entire  latency  interval.  Each  distribution 
is  shifted  with  SQA.  Ihe  OFs  for  the  15  second  delay  conditions  are  also 
plotted  in  the  figure.  For  the  200  msec  SQA  condition  the  entire  CDF  has 
shifted  out  in  time  relative  to  the  0  delay  CDFs.  In  contrast,  the  500  and 
1000  msec  SOA  conditions  produced  CDFs  which  initially  clustered  with  the  0 
delay  conditions.  At  about  the  median  of  each  CDF,  the  functions  deviate 
from  the  0  delay  conditions  and  approach  the  15  second  delay,  200  msec  SOA 
condition.  This  reflects  a  mixture  of  distributions  (Ratcliff,  1979)  and 
supports  the  hypothesis  that  the  1000  and  500  msec  SOA  conditions  reflect  a 
mixture  of  primary  and  secondary  memory  trials.  Furthermore,  if  subjects 
are  retrieving  memory  set  information  in  the  SOA  interval,  one  would  predict 
that  the  proportion  of  primary  memory  trials  would  be  larger  in  the  1000 
msec  SOA  condition.  This  in  fact  appears  to  be  the  case.  The  CDF  for  the 
500  msec  SOA  condition  deviates  from  the  0  delay  conditions  more  rapidly 
than  does  the  1000  msec  SOA  condition.  Hius,  these  data  suggest  that 
subjects  are  retrieving  information  about  the  memory  set  during  the  SOA 
interval  and  that  performance  at  longer  SOAs  reflects  a  mixture  of  primary 
and  secondary  monory  trials.  The  proportion  of  primary  to  secondary  memory 
trials  is  modulated  by  the  SOA  interval. 

To  estimate  the  proportion  of  primary  and  secondary  memory  trials  in 
the  1000  msec  SOA,  15  second  delay  condition,  hypothetical  distributions 
were  generated  by  randomly  selecting  trials  from  the  0  and  15  second  delay, 
200  msec  SOA  conditions.  Figure  3  presents  the  CDFs  for  selected 
proportions  of  primary  and  secondary  memory  trials.  The  CDF  obtained  in  the 
1000  msec  SOA,  15  second  delay  condition  (adjusted  for  SOA)  is  also  plotted 
in  figure  3.  This  condition  falls  between  50/50  and  75/25  secondary/primary 
memory  proportions  and  provides  further  support  for  the  mixture  hypothesis. 

One  effect  which  did  not  attain  significance  was  the  effect  of  response 
type  (target  vs  non-target),  F(l,9)=3.06,  p>.10,  MS^88852.  This  is  consis¬ 
tent  with  previous  research  in  which  the  memory  set  is  changed  after  every 
trial  (Sternberg,  1975).  Further,  the  probability  of  a  mismatched  trial  was 
twice  the  probability  of  a  matched  trial  in  the  running  memory  task  and  this 
may  have  offset  any  response  bias  for  target  trials  in  the  Sternberg  task. 

Conclusions;  Experiment  1 

In  sum,  the  additive  effects  of  memory  load  and  delay  provide  support 
for  the  interpretation  of  separate  retrieval  and  memory  search  processes. 

The  data  further  suggest  that  the  memory  set  information  is  retrieved  as  a 
unit.  These  results  are  in  direct  contrast  to  Sternberg's  findings.  The 
data  from  the  4  second  delay  condition  suggest  that  the  duration  of  the 
distr actor  task  was  not  responsible  for  these  differences.  Hie  manipulation 
of  SQA  provided  evidence  that  subjects  retrieve  sane  of  the  information 
during  the  SQA  interval;  however,  this  retrieval  is  not  complete  prior  to 
the  presentation  of  the  probe.  The  1000  msec  SOA  condition  should  have 
provided  anple  time  to  retrieve  the  information.  Evidence  from  the 
distributions  of  reaction  time  for  the  different  conditions  suggest  that 
performance  in  the  1000  msec  SQA  condition  may  represent  a  mixture  of 
primary  and  secondary  memory  conditions. 
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CUMULATIVE  FREQUENCY 


Target  Load  4 


The  vincentized  cumulative  distribution  functions  for  the 
0  and  15  second  delay  conditions  at  each  SOA  condition.  The 
target  memory  load  A  condition  is  presented.  The  reaction 
time  distribution  was  divided  into  57.  quantiles  and  the  mean 
latency  of  eacn  quantile  was  calculated.  The  means  are 
plotted  as  a  function  of  quantile. 


FIGURE  2 
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Experiment  2;  ERP  Analysis  of  Primary/ Secondary  Menprv  Search 


Experiment  2  examined  the  Event-Related  Brain  Potentials  (ERPs) 
elicited  in  the  primary/secondary  memory  search  task.  This  analysis  was 
conducted  to  determine  if  the  increases  in  reaction  time  as  a  function  of 
delay  in  experiment  1  were  the  result  of  increases  in  stimulus  evaluation 
processing,  response-related  processing,  or  both.  Wickens  et  al.  (1981, 
1985)  hypothesized  that  the  increase  in  reaction  time  from  primary  to 
secondary  memory  was  due  to  the  inclusion  of  a  retrieval  process  inserted 
prior  to  the  memory  conparison  process  in  secondary  memory  conditions.  If 
this  hypothesis  is  correct,  then  stimulus  evaluation  processes  should  be 
delayed  as  subjects  must  retrieve  the  memory  set  information  from  secondary 
memory.  However,  performance  in  secondary  memory  may  be  slower  than  primary 
memory  due  to  shifts  in  response  bias.  Error  rates  were  greater  in 
secondary  memory,  and  subjects  may  have  adopted  a  more  conservative  response 
bias  under  these  conditions.  Since  the  P300  coirponent  of  the  ERP  is 
sensitive  to  stimulus  evaluation  processes,  but  relatively  insensitive  to 
response-related  processes  (Magliero,  Bashore,  Coles,  and  Donchin,  1983; 
McCarthy  and  Donchin,  1981),  the  P30()  can  be  used  to  dissociate  these  two 
factors . 

The  500  msec  SQA  condition  of  experiment  1  was  selected  to  record  ERPs 
from  experiment  2.  The  experimental  parameters  were  identical  to  experiment 
1.  If  P300  latency  increases  as  a  function  of  delay,  this  will  provide 
support  for  the  hypothesis  that  a  retrieval  process  is  inserted  prior  to  the 
memory  comparison  process.  If  P300  latency  does  not  increase  as  a  function 
of  delay,  this  will  provide  support  for  the  hypothesis  that  primary  and 
secondary  memory  differences  are  the  result  of  response  criterion  shifts. 

Methods;  Experiment  2 


Subjects 

Ten  young  subjects  (age  range  18  to  30)  participated  in  the  experiment 
(see  footnote  3).  Subjects  were  right-handed  with  normal  or  corrected-to- 
normal  vision.  Subjects  were  paid  for  their  participation  in  the 
experiment. 

Stimuli  and  Apparatus 

The  stimuli  were  identical  to  experiment  1.  The  stimuli  were  presented 
on  a  Matrox  display  positioned  approximately  70  cm  from  the  subjects.  Ihe 
stimuli  were  presented  within  a  rectangle  in  the  center  of  the  display.  Hie 
rectangle  subtended  a  visual  angle  of  1.0  degrees  vertically  and  1.5  degrees 
horizontally.  Subjects  indicated  their  response  by  depressing  one  of  two 
buttons  on  a  response  box  held  with  both  hands.  Response  button  assignments 
were  counterbalanced  across  subjects. 

Procedure 


Hie  experimental  procedure  was  identical  to  experiment  1,  with  the 
exception  that  only  the  500  msec  SOft  condition  was  included. 
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CUMULATIVE  FREQUENCY 


An  examination  of  different  mixtures  of  primary  and  secondary  memory. 

The  data  are  presented  as  vincentized  CDFs.  The  solid  line  represents 
the  0  second  delay,  200  msec  SOA  condition,  the  dashed  line  represents 
the  15  second  delay,  200  msec  SOA  condition.  The  dotted  line  represents 
the  data  obtained  in  the  15  second  delay,  1000  SOA  condition  (adjusted 
for  the  main  effect  of  SOA).  Three  mixture  conditions  are  presented 
which  represent  25/75,  50/50,  and  75/25  mixtures  of  primary  and  secondary 
memory.  It  is  apparent  that  the  data  obtained  in  the  15  second  delay, 
1000  msec  SOA  condition  falls  between  the  50/50  and  25/75  mixture 
conditions.  FIGURES 
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Experimental  Design 

Hie  experiment  included  three  delay  conditions  (0,  4,  and  15  sec),  and 
two  set  sizes  (2  and  4)  in  a  Sternberg  memory  search  task.  Target  and  non¬ 
target  trials  were  equiprobable.  Thus,  the  study  was  a  3  (delay)  X  2 
(memory  load)  X  2  (response  type)  factorial  design.  The  order  of  within 
subject  conditions  was  counterbalanced  across  subjects  and  sessions. 

ERP  Recording 

The  electroencephalogram  (EEG)  was  recorded  from  three  midline  sites 
(Fz,  Cz,  and  Pz  according  to  the  International  10-20  system;  Jasper,  1958) 
and  referred  to  linked  mastoids.  The  ground  electrode  was  positioned  on  the 
left  side  of  the  forehead.  Electrooculogram  (EOC)  electrodes  were  placed 
above  and  below  the  right  eye.  Electrode  impedances  did  not  exceed  5  KOhms. 
Beckman  10  nm  diameter  Ag/AgCl  biopotential  electrodes  were  used  at  all 
electrode  sites.  Scalp  electrodes  were  affixed  with  Grass  EEG  paste. 
Reference  and  ground  electrodes  were  adhered  with  stomaseal  adhesive 
collars. 

The  EEG  and  EOG  signals  were  amplified  by  Grass  amplifiers  (model 
7P122),  and  filtered  on-line  using  a  high  frequency  cut-off  point  at  35  Hz 
and  a  time  constant  equal  to  8  sec  for  the  high  pass  filter.  Both  EEG  aix3 
EOG  were  sampled  for  2100  msec,  beginning  100  msec  prior  to  retrieval  cue 
onset.  The  data  were  digitized  every  10  msec. 

Stimulus  Generation  and  Data  Collection 

Stimulus  presentation  and  data  acquisition  were  governed  by  a  FDP  11/73 
computer  interfaced  with  a  matrox  display.  Single  trial  EEG  and  E(DG  were 
monitored  on  line.  Digitized  single  trial  data  were  stored  on  magnetic  tape 
for  subsequent  analyses.  EOG  artifacts  were  corrected  off-line  (Gratton, 
Coles,  and  Donchin,  1983). 

Results  and  Discussion;  Experiment  2 

The  reaction  time  results  replicated  the  effects  obtained  in  experiment 
1  and  will  therefore  only  be  discussed  briefly.  Reaction  time  increased  as 
a  function  of  memory  load  (F(l,7)=8.5,  p<.02,  MSeF:4998).  The  memory  set 
size  slope  was  21D.  msec/word.  Reaction  time  also  increased  as  a  function 
of  delay  (F(2,14)=9.1,  p<.01,  MSe=9.1).  The  increase  in  reaction  time  from 
the  0  to  the  15  msec  delay  condition  was  126  msec.  In  contrast  to 
experiment  1,  targets  were  responded  to  faster  than  non-targets 
(F(l,7)=138.6,  p<.001,  MSe=1403). 

Figure  4  a-f  presents  the  grand  average  ERPs  collected  in  the 
experiment.  Each  panel  presents  one  of  the  twelve  experimental  conditions. 
A  description  of  the  experimental  condition  is  presented  above  each  panel. 
Within  each  panel  the  ERPs  for  each  electrode  site  are  over-plotted.  Fz  is 
represented  by  the  solid  line,  Cz  by  the  dashed  line,  and  Pz  by  the  dotted 
line.  The  abscissa  represents  latency,  in  milliseconds,  from  presentation 
of  the  retrieval  cue  (SI).  The  recording  epoch  begins  100  msec  pre  SI.  The 
average  of  this  100  msec  interval  (within  each  electrode  site)  is  used  to 
remove  baseline  shifts.  The  dashed  vertical  line  indicates  the  onset  of  the 
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retrieval  cue  (SL).  The  dotted  vertical  line  r^resents  the  onset  of  the 
Sternberg  probe  stimulus  (S2).  The  ordinate  represents  anplitude  in 
microvolts.  Positive  changes  are  reflected  as  a  downward  deflection.  Hie 
nuirber  in  the  ujper  left-hand  corner  of  each  panel  indicates  the  number  of 
single  trials  in  the  average. 

Perusal  of  the  ERPs  reveals  several  distinct  E3?P  components.  Following 
SL  onset,  an  early  negative  deflection  appears  in  several  of  the  tracings. 
The  latency  of  this  component  is  between  100  and  110  msec  and  the  component 
is  frontal-central  maximal.  We  identify  this  component  as  the  MOO.  The 
second  EPR  component  clearly  identifiable  following  SI  is  P200,  a  positive 
deflection  occurring  afproximately  200  msec  after  SI  onset.  The  third  E3?P 
component  is  MOO,  a  negative  deflection  in  the  ERP  following  P200  with 
frontal-central  maximal  distribution  and  a  latency  from  250  to  350  msec. 

The  last  ERP  component  elicited  by  the  processing  of  the  information 
conveyed  by  SI  is  the  P300.  The  P300  carponent  has  a  parietal  maximal 
distribution  with  a  minimum  latency  of  300  msec.  P300  latency  to  SI  occurs 
between  450  to  500  msec  following  SI  onset. 

Examination  of  the  EE?P  components  elicited  by  the  processing  of  the 
Sternberg  probe  (S2)  is  carpi icated  somewhat  by  baseline  shifts  due  to  SI 
processing.  To  conpensate  for  this,  the  ERPs  were  reaveraged  and  a  new 
baseline  interval  from  onset  of  S2  to  100  msec  post  S2  onset  was  used  to 
remove  baseline  differences.  Figures  5  a-f  present  the  reaveraged  ERPs.  It 
should  be  noted  that  the  different  baselines  do  not  result  in  changes  in  the 
ERP  coirponent  structure,  nor  in  estimates  of  the  latency  of  these 
components.  However,  without  a  proper  baseline,  any  base  to  peak  amplitude 
measures  will  be  biased. 

Examination  of  the  ERP  components  following  S2  reveals  a  component 
structure  similar  to  SL,  The  coiponents  are  MOO  (a  negativity  with  a 
latency  from  100  to  110  msec  post  S2),  P200  (a  positivity  with  a  latency 
from  200  to  250  msec  which  is  somewhat  more  difficult  to  identify,  coipared 
with  the  P200s  elicited  by  M),  MOO  (a  negativity  with  a  latency  from  200 
to  350  msec),  and  P300  (a  positivity  with  a  latency  from  300  to  700  msec). 

It  should  be  noted  that  the  latency  values  are  estimates  of  the  average 
latency  within  a  condition.  The  range  of  the  P300  on  single  trials  is  from 
300  to  1200  msec. 

Because  the  focus  of  experiment  2  was  to  dissociate  stimulus  evaluation 
processing  from  response-related  processing,  our  analyses  will  focus  on  the 
P300  component  of  the  ERP.  We  have  therefore  overplotted  the  Pz  ERPs  for 
the  different  delay  conditions  in  figures  6  and  7.  The  baseline  used  for 
figure  6  a-b  is  from  100  msec  pre  M  to  SI.  The  baseline  used  for  figure  7 
a-b  is  from  S2  to  100  msec  post  S2.  A  description  of  the  experimental 
condition  appears  above  each  panel.  The  solid  line  represents  the  0  delay 
condition,  the  dashed  line  represents  the  4  second  delay  condition,  and  the 
dotted  line  r^resents  the  15  second  delay  condition.  To  examine  latency 
and  amplitude  differences  in  the  P300,  an  estimate  was  derived  for  each 
single  trial.  The  P300  was  identifi^  by  computing  the  correlation  betweei 
the  positive  segment  of  a  .5  Hz  sine  wave  and  the  Pz  electrode  (described  as 
Pz  cross  correlation  in  appendix  K).  The  correlations  were  computed  at  10 
msec  lags  within  an  interval  from  300  to  600  msec  post  SL  and  300  to  1200 
msec  post  S2  for  the  SL  and  S2  P300s,  respectively.  The  point  at  which  the 
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Grand  average  scalp  distribution  ERPs  from  experiment  2.  Ten  subjects 
are  included  in  the  average.  The  baseline  subtracted  from  each  waveform 
is  from  100  msec  pre-Sl  to  SI.  The  number  of  trials  included  in  each 
average  is  displayed  in  the  upper  left-hand  portion  of  each  plot.  The 
solid  line  represents  the  Fz  electrode,  dashed,  Cz,  and  dotted  Pz.  The 
experimental  condition  is  presented  above  each  panel.  Figures  a-f 
contain  the  12  experimental  conditions. 
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maximal  correlation  occurred  was  defined  as  the  latency  of  the  P300.  Ihe 
slope  of  the  correlation  between  the  tenplate  and  the  ERP  was  defined  as  the 
aitplitude  of  the  P300.  The  greater  the  slope,  the  larger  the  amplitude  of 
the  P300.  Note  that  this  procedure  is  unaffected  by  baseline  shifts  and  is 
preferable  to  base  to  peak  measures  in  this  instance. 

Examination  of  the  latency  of  the  P300s  elicited  by  S3,  revealed  a 
significant  response  type  X  memory  load  X  delay  interaction  (F (2 ,14) -11.5, 
p<.01,  MSe=244).  However,  the  effects  of  delay  were  non-monotonic  and  there 
was  no  coherent  pattern  in  the  data.  Because  the  interaction  included  the 
response  type  fac*or  and  the  P300s  were  elicited  prior  to  S2,  this  effect  is 
prctoably  unreliable.  No  other  SI  P300  latency  effects  attained 
significance. 

Examination  of  the  amplitude  of  the  SI  P300s  revealed  that  P300 
amplitude  increased  as  a  function  of  delay  (F(l,7)-4.2,  p<.02,  MSe=229).  Hie 
effect  was  non-mono  tonic,  with  the  smallest  P300s  elicited  in  the  0  delay 
condition.  The  difference  between  the  4  and  15  second  delay  coix3itions  did 
not  significantly  differ  (p<.05).  This  increase  in  P300  amplitude  in  the  4 
and  15  second  delay  conditions  may  reflect  an  increase  in  the  allocation  of 
processing  resources  as  subjects  switch  from  the  distractor  task  to  the 
Sternberg  task.  Additional  research  is  required  to  fully  interpret  these 
differences . 

Examination  of  the  P300s  elicited  by  S2  provides  the  critical  analysis 
of  experiment  2.  P300  latency  increased  .  j  function  of  memory  load 
(F(l,7)=6.6,  p<.05,  MSeF4070).  The  me  ory  sec  size  slope  was  165 
msec/item.  P300  latency  inc'ea^ed  as  a  function  of  delay  (F (2 ,14)  =19 .4, 
p<.001,  MSe=1718).  The  difference  between  the  0  and  15  second  delay 
conditions  was  59  msec.  This  effect  of  delay  on  P300  latency  supports  the 
hypothesis  that  the  increase  in  toacticn  time  from  primary  to  secondary 
memory  is  due  at  least  in  part  to  the  inclusion  of  a  retrieval  process 
inserted  prior  to  the  memory  comparison  process  in  secondary  memory 
conditions.  This  interpretation  is  bolstered  by  the  finding  that  P300 
latency  was  sensitive  to  both  memory  load  and  delay,  suggesting  that  the 
P300  was  elicited  after  the  memory  carparison  process.  No  other  latency 
effects  attained  significance.  In  contrast  to  reaction  time,  P300  latency 
to  targets  was  not  significantly  different  from  non-targets  (F(l,7)=25, 
p<.10,  MSe=5298). 

The  P300  latency  data  suggest  that  the  reaction  time  differences 
between  targets  and  distractors  is  due  to  response- related  processing. 
Furthermore,  the  effect  of  the  delay  was  nearly  twice  as  large  with  reaction 
time  as  with  P300  latency.  This  suggests  that  delay  may  affect  both 
stimulus  evaluation  and  response- related  processing.  To  evaluate  the 
contribution  of  response-related  processing,  a  ccmparison  of  the  single 
trial  Kr/P300  ratio  was  conducted.  By  comparing  the  relative  changes  in 
reaction  time  and  F300  latency,  it  is  possible  to  localize  the  effects  of 
different  manipulations  on  information  processing.  Increases  in  response- 
related  processing  are  identified  as  increases  in  the  Rr/P300  ratio.  A 
ratio  of  1.0  reflects  the  co-occurrence  of  reaction  time  and  the  peak  of  the 
P300.  Ratios  less  than  1.0  indicate  that  the  reaction  time  preceded  peak 
P300  latency  and  ratios  greater  than  1.0  indicate  that  the  peak  P300  latency 
preceded  the  overt  response.  It  is  important  to  note  that  absolv-te  RT/P300 
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Grand  average  scalp  distribution  ERPs  from  experiment  2.  Ten  subjects 
are  included  in  the  average.  The  baseline  subtracted  from  each  waveform 
is  from  S2  to  100  msec  post-S2.  The  number  of  trials  included  in  each 
average  is  displayed  in  the  upper  left-hand  portion  of  each  plot.  The 
solid  line  represents  the  Fz  electrode,  dashed,  Cz,  and  dotted,  Pz.  The 
experimental  condition  is  presented  above  each  panel.  Figures  a-f  contain 
the  12  experimental  conditions. 

FIGURE  5  A 
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FIGURE  5  D 


FIGURE  5  E 


FK3URE  5  F 


ratios  are  less  meaningful  than  the  relative  differences  between 
experimental  conditions.  This,  in  part,  is  due  to  the  fact  that  we  measure 
the  latency  of  the  peak  of  the  P300  rather  than  the  onset  of  the  P300. 

Examination  of  the  single  trial  Rr/P300  ratio  revealed  an  increase  from 
targets  to  non-targets  (F(l,7)=173,  p<.01,  MSe=11554),  suggesting  an 
increase  in  the  response- related  processing  for  non-targets.  This  increase 
in  response- related  processing  could  be  due  to  response  biases,  uncertainty, 
rechecking,  or  some  combination  of  these  factors.  The  Kiyp300  ratio 
increased  as  a  function  of  delay  (F(2,14)=5.3,  p<.01,  MS&=15266).  This  is 
consistent  with  the  hypothesis  that  some  of  the  increase  in  reaction  time  as 
a  function  of  delay  is  due  to  response-related  processing.  Finally,  the 
Rr/P300  ratio  increased  at  a  greater  rate  for  non-targets  than  for  targets 
(F(2,14)=73,  p<.006,  MSe=6569).  Post  hoc  comparisons  revealed  that  this 
interaction  was  due  to  a  large  increase  in  the  Rr/P300  ratio  for  the  non¬ 
target  delay  15  condition  and  may  be  the  result  of  uncertainty  about  the 
response,  or  to  a  rechecking. 

Examination  of  the  amplitude  of  the  P300s  elicited  by  S2  revealed  a 
significant  response  type  X  delay  interaction  (F(2,14)=33,  p<U)5,  MSe=101). 
Post  hoc  analysis  revealed  that  P300  amplitude  was  smaller  for  the  15  second 
delay  condition  for  targets  and  larger  for  the  15  second  delay  condition  for 
non- targets. 


Conclusion!  Experiment  2 

Experiment  2  was  conducted  to  examine  the  hypothesis  that  the  increase 
in  reaction  time  from  primary  to  secondary  memory  was  due  to  the  insertion 
of  a  retrieval  process  prior  to  the  memory  comparison  process,  to  response- 
related  processing,  or  both.  The  evidence  obtained  using  P300  latency  as  an 
index  of  stimulus  evaluation  processing  suggests  that  both  stimiulus 
evaluation  and  response- related  processing  factors  are  affected  by  delay. 
These  results  suggest  that  the  retrieval  of  information  from  secondary 
memory  may  be  more  rapid  than  the  reaction  time  estimates  suggest.  The 
results  also  suggest  that  secondary  memory  includes  a  retrieval  process 
inserted  prior  to  the  memory  comparison  process. 
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Grand  average  Pz  overplots  of  delay.  Ten  subjects  are  included  in  the  average. 
The  baseline  subtracted  from  each  waveform  is  from  S2  to  100  msec  post-S2.  The 
solid  line  represents  the  0  delay  condition,  the  dashed  line  represents  the  4 
second  delay  condition,  and  the  dotted  line  represents  the  15  second  delay 
condition.  The  experimental  condition  is  described  above  each  panel. 
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FOCTTNOTES 


Footnote  1.  The  terms  primary  and  secondary  memory  originate  with  William 
James  (1890,  p.  646-648).  Contemporary  treatments  can  be  found  by  Waugh  and 
Norman  (1965),  Craik  and  Levey  (1976),  Atkinson  and  Shiffrin  (1968), 

Baddeley  and  Hitch  (1974),  and  Baddel^  (1981).  According  to  James,  primary 
memory  reflects  the  contents  of  consciousness,  while  secondary  memory  refers 
to  information  that  no  longer  is  in  consciousness.  Following  Waugh  and 
Norman  (1965)  primary  memory  can  be  defined  in  terms  of  its  limited  capacity 
and  the  rapid  decay  of  information  if  rehearsal  is  not  permitted.  Rehearsal 
serves  to  maintain  information  in  primary  memory  and  transfer  the 
information  to  secondary  memory.  Secondary  memory  is  characterized  as  a 
large  capacity  storage  system  with  a  low  forgetting  rate. 

Footnote  2.  Care  was  ernployed  selecting  the  words  so  as  not  to  form  any 
obvious  categories  and  subjects  were  instructed  not  to  use  elaborative/ 
chunking  strategies  when  memorizing  the  memory  set.  These  manipulations 
were  intended  to  minimize  any  effectiveness  of  the  conpression  of  the  menory 
set  into  subgroups  of  categories.  Furthermore,  since  this  was  a  VM  task,  it 
would  be  quite  difficult  for  subjects  to  form  associations  for  every 
pairwise  and  four-way  comparison. 

Footnote  3.  TVo  subjects  (subjects  8  and  9,  see  ajpendix  N)  were  excluded 
from  the  analysis  because  of  excessive  error  rate,  resulting  in  a  total 
sanple  size  rt  r 
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Appendix  A 


PMSM. PAS  This  is  the  Pascal  source  code  for  the  primary/secondary  memory 
experiment.  The  program  is  written  in  Turbo  Pascal  and  assumes  an  enhanced 
graphics  adapter  (EGA)  card.  The  program  should  be  compiled  prior  to 
initial  execution.  Two  INCLUDE  files  are  included  at  compile  time:  CRT. PAS 
and  TIMER. PAS.  The  source  for  these  procedures  can  be  found  in  appendix  B 
and  C,  respectively.  The  inclusion  of  these  procedures  is  accomplished  by 
the  {$!  CRT. PAS}  and  {$!  TIMER. PAS}  commands  at  the  beginning  of  the 
program. 

A  description  of  the  experimental  procedure  can  be  found  in  the  methods 
section  of  the  paper.  Instructions  to  subjects  are  included  as  appendix  M. 
To  run  the  program  an  input  stimulus  file  and  an  output  data  file  must  be 
specified  e.g.,  Type:  PMSM  <input  file>  <output  file>.  The  program  will  not 
proceed  unless  the  stimulus  input  file  exists  on  disc  and  the  output  data 
file  does  not  exist  on  disc.  A  sample  stimulus  input  file  is  included  as 
appendix  F.  This  file  contains  20  word  stimuli,  the  response  assignment  keys 
for  'Yes'  and  'No'  responses,  and  96  primary/secondary  memory  trial  blocks. 
The  program  presents  the  96  trial  blocks  (run  time  approximately  20  minutes) 
and  then  writes  108  identification  parameters  for  each  block  to  the  output 
file.  Program  PMSM  pauses  prior  to  the  first  block  and  after  every  16  trial 
blocks  to  provide  rest  breaks.  Subjects  continue  the  program  by  pressing 
the  enter  key.  Several  options  may  be  specified  by  pressing  designated  keys 
during  run  time: 

'Q'  will  quit  the  program  and  write  the  output  to  disc. 

'P'  will  pause  the  program.  To  continue,  press  enter. 

'F'  will  provide  trial  by  trial  feedback  to  subjects.  Once  this  option  is 
specified,  it  remains  in  effect  throughout  the  96  trial  blocks. 
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{  Modification  date:  Jan  12,  1988 
Programmer:  David  Stray er 
University  of  Illinois 

Program  to  control  primary/secondary  memory  study 

ids  output  by  the  program: 

id(l)  -  subject  number 
id(2)  -  session 

id(3)  -  block 

id(4)  -  delay 

id(5)  -  soa  of  retrieval  cue 

id(6)  -  size  of  memory  set 

id(7)  -  constant  (1)  or  varied  (2)  mapping 

id(8)  -  ASCII  code  for  Yes_button 

id(9)  -  ASCII  code  for  No__button 

Running  memory  ids 


Trial  1 

id(lO)  - 
Id(ll)  - 
id(12)  - 
ld(13)  - 
ld(14)  - 


Trial  2 

id(15)  - 
id(16)  - 
ld(17)  - 
id(18)  - 
id(19)  - 


Trial  3 

ld(20)  - 
id(21)  - 
ld(22)  - 
id(23)  - 
ld(24)  - 


Trial  4 

ld(25)  - 
id(26)  - 


running  memory  stimulus 
rt  (for  running  memory  trial) 

stim  type  for  running  memory  (l=target,  2-nontarget) 
response  for  running  memory  (0-none,  l*=target,  2-nontarget) 
correctness  for  running  memory:  0-no  response,  1-hit, 
2=correct  rejection,  3-false  alarm,  4=miss 


running  memory  stimulus 
rt  (for  running  memory  trial) 

stim  type  for  running  memory  (l=target,  2-nontarget) 
response  for  running  memory  (0=none,  l-target,  2-nontarget) 
correctness  for  running  memory:  0-no  response,  1-hit, 
2-correct  rejection,  3-false  alarm,  4-miss 


running  memory  stimulus 
rt  (for  running  memory  trial) 

stim  type  for  running  memory  (l-target,  2-nontarget) 
response  for  running  memory  (0-none,  l-target,  2-nontarget) 
correctness  for  running  memory:  0-no  response,  1-hlt, 
2-correct  rejection,  3-false  alarm,  4-mi88 


running  memory  stimulus 
rt  (for  running  memory  trial) 
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id(27) 

id(28) 

id(29) 


Trial  5 

id(30) 
id(31) 
id(32) 
id(33) 
id( 34 ) 


Trial  6 

id(35) 

id(36) 

id(37) 

id(38) 

id(39) 


Trial  7 

id(40) 
id(41) 
id(42) 
id(43) 
id(44 ) 


Trial  8 

id(45) 

id(46) 

ld(47) 

id(48) 

id(49) 


Trial  9 

id(50) 

id(51) 

id(52) 

id(53) 

ld(54) 


stim  type  for  running  memory  (l=target,  2=nontarget) 
response  for  running  memory  (0=none,  l=target,  2=nontarget) 
correctness  for  running  memory:  0=no  response,  l=hit, 
2=correct  rejection,  3=false  alarm,  4=miss 


running  memory  stimulus 
rt  (for  running  memory  trial) 

stim  type  for  running  memory  (l=target,  2=nontarget) 
response  for  running  memory  (0=none,  l=target,  2=nontarget) 
correctness  for  running  memory:  0=no  response,  l=hit, 
2=correct  rejection,  3=false  alarm,  4=miss 


running  memory  stimulus 
rt  (for  running  memory  trial) 

stim  type  for  running  memory  (l=target,  2=nontarget) 
response  for  running  memory  (0=none,  l=target,  2»=nontarget ) 
correctness  for  running  memory:  0=no  response,  l=hlt, 
2=correct  rejection,  3=false  alarm,  4=miss 


running  memory  stimulus 
rt  (for  running  memory  trial) 

stim  type  for  running  memory  (l=target,  2=nontarget) 
response  for  running  memory  (0=none,  l=target,  2=nontarget) 
correctness  for  running  memory:  0=no  response,  l=hit, 
2=correct  rejection,  3=false  alarm,  4=miss 


running  memory  stimulus 
rt  (for  running  memory  trial) 

stim  type  for  running  memory  (l=target,  2=nontarget) 
response  for  running  memory  (O=none,  l=target,  2=nontarget) 
correctness  for  running  memory:  0=no  response,  l^hit, 
2=correct  rejection,  3=false  alarm,  4=miss 


running  memory  stimulus 
rt  (for  running  memory  trial) 

stim  type  for  running  memory  (l=target,  2-nontarget) 
response  for  running  memory  (0=none,  l=target,  2=nontarget) 
correctness  for  running  memory:  0“no  response,  l»hit, 
2=correct  rejection,  3=false  alarm,  4=mi6S 


Trial  10 


id(55)  - 
id(56)  - 
ld(57)  - 
id(58)  - 
ld(59)  - 


Trial  11 

id(60)  - 
id(61)  - 
id(62)  - 
id(63)  - 
id(64)  - 


Trial  12 

id(65)  - 
id(66)  - 
ld(67)  - 
id(68)  - 
ld(69)  - 


Trial  13 

id(70)  - 
id(71)  - 
id(72)  - 
id(73)  - 
ld(74)  - 


Trial  14 

id(75)  - 
ld(76)  - 
id(77)  - 
id(78)  - 
id(79)  - 


Trial  15 

id(80)  - 
id(81)  - 
id(82)  - 
ld(83)  - 
id(84)  - 


running  memory  stimulus 
rt  (for  running  memory  trial) 

stim  type  for  running  memory  (l=target,  2=nontarget) 
response  for  running  memory  (0*none,  l*target,  2“nontarget) 
correctness  for  running  memory:  0“no  response,  l“hit, 
2=correct  rejection,  3=false  alarm,  4“mlss 


running  memory  stimulus 
rt  (for  running  memory  trial) 

stim  type  for  running  memory  (l^target,  2-nontarget) 
response  for  running  memory  (0=none,  1-target,  2-nontarget) 
correctness  for  running  memory:  0-no  response,  1-hit, 
2=correct  rejection,  3=false  alarm,  4-mi8s 


running  memory  stimulus 
rt  (for  running  memory  trial) 

stim  type  for  running  memory  (1-target,  2-nontarget) 
response  for  running  memory  (0=none,  1-target ,  2-nontarget) 
correctness  for  running  memory:  0»no  response,  l-hit, 
2=correct  rejection,  3-false  alarm,  4-miss 


running  memory  stimulus 
rt  (for  running  memory  trial) 

stim  type  for  running  memory  ( 1-target ,  2-nontarget) 
response  for  running  memory  (0=none,  1-target,  2-nontarget) 
correctness  for  running  memory:  0-no  response,  1-hit, 
2-correct  rejection,  3-false  alarm,  4-miss 


running  memory  stimulus 
rt  (for  running  memory  trial) 

stim  type  for  running  memory  (1-target,  2-nontarget) 
response  for  running  memory  (0-none,  1-target,  2-nontarget) 
correctness  for  running  memory:  0-no  response,  1-hit, 
2-correct  rejection,  3-false  alarm,  4-ml8S 


running  memory  stimulus 
rt  (for  running  memory  trial) 

stim  type  for  running  memory  (1-target,  2-nontarget) 
response  for  running  memory  (0-none,  1-target,  2-nontarget) 
correctness  for  running  memory;  0-no  response,  1-hlt, 
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2=correct  rejection,  3=false  alarm,  4=mlss 
Sternberg  ids 

id(85)  -  memory  set  word  //I,  letter  #1 

id(86)  -  memory  set  word  #1,  letter  //2 

ld(87)  -  memory  set  word  #1,  letter  //3 

id(88)  -  memory  set  word  #1,  letter  #4 

id(89)  -  memory  set  word  #2,  letter  #1 

id(90)  -  memory  set  word  i'/2,  letter  y/2 

ld(91)  -  memory  set  word  #2,  letter  #3 

id(92)  -  memory  set  word  #2,  letter  //4 

ld(93)  -  memory  set  word  ?i*3,  letter  #1 

id(94)  -  memory  set  word  #3,  letter  #2 

id(95)  -  memory  set  word  //3,  letter  //3 

id(96)  -  memory  set  word  #3,  letter  //4 

id(97)  -  memory  set  word  #4,  letter  #1 

id(98)  -  memory  set  word  #4,  letter  #2 

id(99)  -  memory  set  word  #4,  letter  #3 

Id(lOO)  -  memory  set  word  //4,  letter  #4 

id(lOl)  -  rt  (for  Sternberg  trial) 

id(102)  -  stim  type  (l=target,  2=nontarget) 

id(103)  -  response  (0=none,  l=target,  2=nontarget) 

id(104)  -  correctness:  0=no  response,  l=hit,  2=correct  rejection, 

3=false  alarm,  4=miss 

id(105)  -  stimulus  presented,  letter  //I 
id(106)  -  stimulus  presented,  letter  ill 
id(107)  -  stimulus  presented,  letter  #3 
id(108)  -  stimulus  presented,  letter  /i*4 
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Program  pmsm  ; 

{$C-}  {disable  cntrl  c  option} 

($1  TIMER. PAS}  (include  timer  routines} 

($1  CRT. PAS}  (CRT  control  routines  (synch  and  cursor  manipulation)} 


var 


ID  :  array [ 1 . .96 , 1 .. 108]  of  integer; 
Condition_codes  :  array} 1 . .96 ,1 . .7]  of  integer; 

CM_targ:  array[ 1 . .5 , 1 . .4]  of  char; 

CM_dist:  array [ 1 . .5 , 1 . .4 ]  of  char; 

VM_stimuli:  array] 1 . . 10 , 1 . .4]  of  char; 

Targ:  array [ 1 . . 5 , 1 . .4]  of  char; 

Dist:  array [ 1 . .5 , 1 . .4]  of  char; 


stmdur,  (duration  of  stimulus} 

nback,  (number  back  in  running  memory  task} 

block,  (block  number} 

subject,  (subject  number} 

session,  (session  number} 

soa,  (stimulus  onset  asynchrony} 

cmvm,  (consistent  (1)  or  varied  (2)  mapping} 

memsiz,  (memory  set  size} 

stype,  (stim  type  (p^  or  neg)  for  Stern  trial} 
delay:  Integer;  (  delay  i.rom  memory  set  to  probe} 


Yes_button;  char;  (key  to  press  for  positive  response} 
No_button:  char;  (key  to  press  for  negative  response} 

s  :  TimeStamp;  (Timer  variable  for  TimerMllli  calls} 


Numf lie, 

Newf lie. 

Exists , 
scrnon, 
pause, 
fback, 

done  :  Boolean; 


(true 

if 

(true 

if 

(true 

if 

(true 

if 

(true 

if 

(true 

if 

(true 

if 

number  of  files  is  2} 

output  file  does  not  exist  on  disc} 

input  file  exists  on  disc} 

stim  on,  false  if  stim  off} 

P  key  pressed} 

F  key  pressed} 

Q  key  pressed} 


Infile, 

Outfile:  Text; 


41 


Function  RandomTwixt( low, high:  Integer):  Integer  ; 


{  Procedure  to  return  a  random  integer  between  the  values  of  low  and  high  } 
begin 

RandomTwixt  :=  Random(high-low+l )+low; 
end ; 
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P-ocedure  Init_all;  {  intializes  variables  } 


var 

1. 

j:  integer; 

begin 

CrtCursorOff ;  {turn  off  the  cursor} 

clrscr;  {clear  the  screen} 

stmdur  :=  200;  {stimulus  duration} 

nback  :=  2;  {number  back  in  r-memory} 

block  :*  0;  {block  counter} 

Timerlnlt;  {initializes  the  clock} 

done  :=  false;  {true  if  Q  key  pressed} 

fback  :=  false;  {true  if  F  key  pressed} 
pause  :=  true;  {true  if  P  key  pressed} 

for  1  :=  1  to  96  do 
for  J  :=  1  to  108  do 
begin 

ld[iij]  0;  {zero  id  array} 

end; 

end; 
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Procedure  Read_stim; 
var 

i, 

j:  integer; 
begin 


for  i  :=  1  to  5  do  {read  in  CM  target  stimuli} 

begin 

readln(Inf ile,CM_targ [ i] ) ; 
end ; 

for  i  :=  1  to  5  do  (read  in  CM  distractor  stimuli} 

begin 

readlndnf  ile,CM_dist[il ) ; 
end ; 

for  i  ;=  1  to  10  do  (read  in  VM  stimuli} 

begin 

readln(Infile,VM_stimuli[i] ); 
end; 


readln( Inf ile,Yes _ button, No_button) ;  (read  key  assignments} 

for  j :=  1  to  96  do  {read  in  condition  codes} 

begin 

for  i:=  1  to  7  do  read( Inf ile, condition  codes[j,i]); 
if  j  <  96  then  readln(lnf ile) ; 

end; 

end;  {read  stim} 
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Procedure  Stlm_Sel;  {shuffle  stimuli  and  choose  targ  and  dlst  stimuli} 
var 

tempi  :  char; 
loop,  (loop  counter} 

1, 

j, 

k  :  Integer;  {variable  equal  to  a  randomly  selected  array  position} 


begin 

block  :=  block+1;  {Increment  block  counter} 

If  block  =  96  then  done  true;  {exit  program  after  this  block} 

If  (block  mod  16)  =  0  then  pause  ;=  true; 


subj  ect 

session 

soa 

delay 

cmvm 

memslz 

stype 


condition 

condition 

condltlon_ 

condltlon_ 

condltlon_ 

condition 

condition 


codes 

codes 

codes 

codes 

codes 

codes 

codes 


[block, 1] 
[block, 2] 
[block, 3] 
[block, 4] 
[block, 5] 
[block, 6] 
[block,?] 


{store  Id  values} 


Id [ block, 1] 
Id [block, 2] 
ld[block,3] 
ld[block,4] 
ld[block,5i 
Id [block, 6] 
ld[block,7] 
ld[block,8] 
ld[block,9] 


*  subject; 

»  session; 

=  block; 

=  delay; 

“  soa; 

>  memslz; 

»  cmvm; 

“ord(Yes_button) ; 
“ord(No  button); 


for  k  :*  1  to  10  do 
begin 

for  loop  1  to  5  do  [permute  CM  targets} 
begin 


j  RandomTwlxt( 1 ,5) ; 
for  1  :■  1  to  4  do 
begin 

tempi  CM_targ[j,l]; 
CM_targ[j,l]  :■  CM_targ[loop,l] ; 
CM_targ[loop,l]  tempi; 
end;  {for} 
end;  {for} 
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for  loop  :=  1  to  5  do  {permute  CM  distractors} 
begin 

j  :=  RandoraTwlxt( 1 ,5) ; 
for  i  :=  1  to  4  do 
begin 

tempi  :=  CM_dist [ j , i ] ; 

CM_dist[j,i]  :=  CM_dist [ loop,i ] ; 

CM_dist'(loop,i]  :=  tempi; 
end;  {for} 
end;  {for} 

for  loop  :=  1  to  10  do  {permute  VM  array} 
begin 

j  :=  RandomTwixt( 1 , 10) ; 
for  i  :=  1  to  4  do 
begin 

tempi  :=  VM_stimuli [ j ,i] ; 

VM_stimuli[ j ,i]  ;=  VM_6timulitloop,i] ; 
VM_stimuli[loop,i]  :=  tempi; 
end;  {for} 
end;  {for} 

end;  {for} 

if  cmvm  =  1  then  {choose  CM  stimuli} 
begin 

for  j  :=  1  to  memsiz  do 
for  i  :=  1  to  4  do 
begin 

targ[j,i]  :=  CM_targ [ j ,i] ;  {target  array} 
end; 

for  j  :=  1  to  5  do 
for  i  ;=  1  to  4  do 
begin 

dist[j,i]  :=  CM_dist[j ,i] ;  {distractor  array} 
end; 

end; 

if  cmvm  =  2  then  {choose  VM  stimuli} 
begin 

for  j  ;»  1  to  memsiz  do 
for  i  ;*  1  to  4  do 
begin 

targ[j,i]  ;»  VM_stimulit j ,i] ;  {target  array} 
end; 


for  j  ;=  1+memslz  to  S+memslz  do 
for  i  ;=  1  to  4  do 
begin 

dist [ j-memsiz , i 1  :=  VM_8timuli [ j , 1 ] ;  {dlstractor  array} 
end; 

end; 


end;  {stimulus  select} 
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Procedure  Memory_Set;  {procedure  to  present  Memory  Set} 

var 

j  . 

i:  integer; 
word:  string[4] ; 
response:  char; 
begin 

if  keypressed  then 
begin 

read(kbd, response) ; 
if  upcase( response) 
if  upcase( response) 
if  upcase( response) 
end;  {if/then} 

for  i  :=  1  to  memsiz  do 

begin  {presentation  of  memory  set} 

if  memsiz  =  2  then  gotoXY(27+(i*8) ,8) ;  {position  cursor} 

if  memsiz  =  4  then  gotoXY( l9+( i*8) ,8) ;  {position  cursor} 

word  :=  targ[i,l]+targ[i,2]+targ[i,3]+targti,4] ; 
writeln(trm,word) ;  {turn  on  memory  set} 
for  j  :=1  to  4  do  id(block,80+(4*i)+j ]  :=ord( targ[ i,j ] ) ; 

end; 

TimerStamp(s) ;  {start  the  clock} 

while  TlmerMilli(s)  <  3000  do; 

clrscr;  {turn  off  memory  set} 


while  TlmerMilli(s)  <4500  do; 
end;  {memory  set  procedure} 


{clear  the  Kbd  buffer} 

=  'Q'  then  done  :=  true; 

=  'F'  then  fback  :■*  true; 
=  'F'  then  pause  :=  true; 
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Procedure  Run  Trial; 


{procedure  to  present  Run  trial  } 


var 

RT, 

pos_neg , 
itrl, 

jstim:  integer; 

Istim:  array[1..3]  of  integer; 
response  ;  char; 
begin 

if  delay  <>  0  then  (if  delay  is  0,  bypass  running  memory  trials} 
begin 

for  itrl;=  1  to  Abs(delay)  do 
begin 

if  keypressed  then  {clear  the  Kbd  buffer} 

begin 

read(kbd, response) ; 

if  upcaseC response)  =  'Q'  then  done  :=  true; 
if  UpcaseC  response)  =  'F'  then  fback  :=»  true; 
if  upcase(response)  =  then  pause  :=  true; 
end;  {if /then} 

RT  :=  -1;  {set  RT  =  -1} 

TimerStamp(s) ;  {start  the  clock} 

{  prepare  stimuli  during  pre-stimulus  interval  (5  msec)} 

pos_neg  RandomTwixtC  1 ,3)  ; 
if  pos_neg  =  3  then  pos_neg  ;=  2; 

{  if  a  positive  trial,  repeat  stimulus  nback 

if  a  negative  trial,  select  a  new  stimulus} 


if  pos_neg  =  1  then  {positive  trial} 
begin 

if  itrl  <=  nback  then  jstim  RandomTwixtCl ,9) ; 
if  itrl  >  nback  then  Jstim  ;■  l8tim[nback] ; 
end; 

if  pos_neg  ■»  2  then  {negative  trial} 
begin 

jstim  RandomTwixtC  1 ,9) ; 
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while  jstim  =  lstim[nback]  do 
begin 

jstim  :=  RandomTwixt( 1 ,9) ; 
end ; 

end;  {if/then} 

{store  successive  stimuli} 

lstim[3]  :=  lstim[2]; 
lstim[2]  :=  lstim[l]; 
lstim[ 1]  :=  jstim; 
id [ block, ( it rl*5)+5 ]  ;=  jstim; 


while  TimerMllli(s)  <  5  do;  (wait  out  prestimulus  interval} 


gotoXY(40,ll);  (pos  ition  cursor} 

CrtSyncEnd;  (sync  CRT} 

TimerStamp(s) ;  (start  the  clock} 

if  delay  >  0  then 

writeln(trm,jstim) ;  {turn  on  stimulus} 

scrnon  :=  true;  {set  screen  flag} 

while  TimerMilli(s)  <  1000  do 

begin 

if  timerMilli(s)  >=  stmdur  then 
if  scrnon  then 
begin 

clrscr;  {turn  screen  off} 

scrnon  :=  false;  {reset  screen  flag} 

end ; 


if  RT  -  -1  then  {if  no  response  has  been  made} 

if  keypressed  then  {check  for  keypress} 

begin 

RT  :=  TimerMilli(s);  {set  RT  =  elapsed  time} 

readCkbd, response);  {read  which  key  was 

pressed} 

end;  {if /then} 

end ; 

{  classify  trial  and  store  RT  and  stim  type  (l=hit,  2“cr,  3=fa,  4=miss)} 
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id[block,(itrl*5)+6]  rt; 
id[block,(ltrl*5)+7]  :■  pos_neg; 

id[block,(itrl*5)+8]  :•=  0; 

id[block,(itrl*5)+9]  :»  0; 


if  rt  >  0  then 
begin 

if  upcaseC response)  «*  Yes_button  then 
begin 

id[block,(itrl*5)+8]  1; 

if  id[block,(itrl*5)+7j  ■  1  then  id[block,(ltrl*5)+9]  :«  1; 

if  id[block,(itrl*5)+7]  ■  2  then  id[block,(itrl*5)+9j  :■  3; 

end; 

if  upcase( response)  =  No_button  then 
begin 

id[block,(itrl*5)+8]  2; 

if  id[block,(itrl*5)+71  -  2  then  idtblock,(itrl*5)+9)  :•  2; 

if  idfblock,(itrl*5)+7]  »  1  then  id(block,(itrl*5)+9i  ;■  4; 

end; 

end;  {if /then} 
end;  {for} 
end;  {if/ then} 

end;  {running  memory  procedure} 
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Procedure  Stern_Trial;  {procedure  to  present  Stern  trial  } 
var 


RT, 

j, 

i:  integer; 
word:  string[4]; 

response:  char; 
begin 


if  keypressed  then  (clear  the  Kbd  buffer} 

begin 

read(kbd , response) ; 

if  upcase(response)  =»  'Q'  then  done  :=  true; 
if  upcase( response)  =  'F'  then  fback  :=  true; 
if  upcase( response)  =  'P'  then  pause  :=»  true; 
end;  {if/then} 

RT  :=  -1;  (set  RT  *  -1} 


{  deliver  retrieval  cue  if  appropriate  } 


if  soa  >  0  then 
begin 

gotoXY(40 , 11 ) ; 
CrtSyncEnd; 
TimerStamp( s ) ; 
writeln( trm,'*' ) ; 
gotoXY(40,ll) ; 


(set  position  on  screen} 
{sync  CRT} 

{start  the  clock} 

{turn  on  retrieval  cue} 
{set  position  on  screen} 


while  TimerMilli(s)  <  stmdur  do;  {wait  for  end  of  stmdur} 
writeln(tnn,'  ');  {turn  off  retrieval  cue} 


while  TimerMilli(s)  <  soa  do; 


end;  {if/then} 

{  prepare  stimuli  during  pre-stimulus  interval  (5  msec)} 

TimerStamp(s) ;  {start  the  clock} 

gotoXY( 39 , 13) ;  {set  position  on  screen} 

if  stype  =  1  then  {positive  trial} 
begin 

i  :=■  RandomTwlxt(  1  ,memsiz)  ; 

word  :=  targ [ i , 1 ] +targ [ i ,2 ]+targ [ i ,3]+targ [ i ,4 ] ; 
for  j  :  =  1  to  4  do  id  [  block,  104+J  ]  :“«ord(  targ  [  1 ,  j  ] ) ; 
end; 

if  stype  =  2  then  {negative  trial} 
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begin 

1  :=  RandomTwlxtC 1 ,5) ; 

word  :»  dist[i, l]+dlst[i,2]+dl8t[i,3]+dlst[i,4] ; 
for  j  :-l  to  4  do  id (block, 104+j ]  ;-ord(di8t(i,j ] ) ; 
end; 

while  TlmerMllll(s)  <  5  do;  {wait  out  prestimulus  interval} 


CrtSyncEnd; 
TlmerStamp(s) ; 
writeln(tnn,word) ; 
scrnon  :•  true; 


{sync  CRT} 

{start  the  clock} 
{turn  on  stimulus} 
{set  screen  flag} 


while  TlmerMilll(s)  <  3000  do 
begin 

if  timerMilli(s)  >••  stmdur  then 
if  scrnon  then 
begin 

clrscr;  {turn  screen  off} 

scrnon  :*  false;  {reset  screen  flag} 

end; 


if  RT  -  -1  then 


pressed} 


if  keypressed  then 
begin 

RT  TlmerMilli(s); 
read(kbd, response) ; 


{if  no  response  has  been  made} 
{check  for  keypress} 


{set  RT  -  elapsed  time} 
{read  which  key  was 


end;  {if /then} 


{  classify  trial  and  store  RT  and  stim  type  (I»hlt,  2“cr,  3-fa,  4-mls8)} 


ld[block,101] 
id [block, 102] 
ld[block,103] 
id [block, 104] 


-  rt; 

-  stype; 

-  0; 

-  0; 


if  rt  >  0  then 


begin 
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if  upcaseC response)  =  Yes_button  then 
begin 

idfblock, 103]  :=  1; 

if  id [block, 102]  =  1  then  id[block, 104]  ;=  1; 

if  id[block,1021  =■  2  then  id [ block, 104 ]  :=  3; 

end ; 

if  upcase( response)  =  No__button  then 
begin 

id [ block, 103]  :=  2; 

if  id  I  block , 102 )  =  2  then  id [ block , 104]  :=  2; 

if  id[block,102]  =  1  then  id[block,104]  :=  4; 

end ; 

end;  {if/then} 
end;  {sternberg  procedure] 
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Procedure  Pause_program;  {procedure  to  pause  program) 
begin 

if  pause  Chen 
begin 
clrscr; 

gotoXY(20,20); 

if  block  -  0  then  write(trm,'  Hit  ENTER  to  initiate  program  '); 
if  block  >  0  then  wrlCe(trm, 'Program  paused:  Hit  ENTER  Co  continue'); 
while  not  keypressed  do; 
clrscr; 

pause  false; 
end; 
end; 


55 


Procedure  Feed  back; 


{procedure  Co  present  feedback} 


var 
error, 
corr , 

i:  integer; 
begin 

error  :=  0; 
corr  :=  0; 

if  fback  then 
begin 

if  delay  <>  0  then 
begin 

for  i  :=  3  to  Abs(delay)  do 
begin 

case( id [ block, ( i*5)+9 ] )  of 
0,3,4:  error  :=  error+1; 

1,2:  corr  :=  corr+1 ; 
end ; 
end; 

gotoXY(25,23) ; 

write(trm,'  Run  —  ',corr,'  correct  trials  out  of  ' ,corr+error) 
end; 

gotoXY (25,24); 


case 

( id [block , 104 

] )  of 

0: 

write 

( t  rm , ' 

Stern  — 

No  response  detected') 

1,2: 

write 

(trm,' 

Stern  — 

Correct  response  ') 

3,4: 
end ; 
end ; 

write 

( tnn,' 

Stern  — 

Incorrect  response  ') 

end ; 
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Procedure  Sesslon__end;  {operations  following  session  termination} 


var 

lines, 

1, 

j :  integer; 
begin 

tlmerklll;  {kill  the  clock} 

for  1  ;=  1  to  96  do 
begin 

lines  :*  1; 
for  j  :=  1  to  108  do 
begin 

write(Outflle,id(l,j ] :5);  {write  data  to  disc} 
if (lines  mod  16)  ■  0  then  wrlteln(0utf lie) ; 
lines  llnes+l; 
end; 

writeln(0utf lie) ; 
end; 

CrtCursorOn;  {turn  the  cursor  back  on} 

Close  (Infile); 

Close  (Outfile); 

end; 
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Begin  {Main  program  loop} 

Numfile  :=  false; 

if  ParamCount  <>  2  then 
begin 

writeln(trm, 'ERROR:  Incorrect  number  of  parameters'); 

writeln( trm) ; 

writeln(trm, 'Correct  format:  pmsra  <stimulus  file>  <output  file>') 
end ; 

if  ParamCount  =  2  then 
begin 

Numfile  :=  true; 

Assign( Inf lie  jParamStrC 1 ) ) ; 

{$!-}  Reset(Infile)  {$!+}; 

Exists  :=  (lOresult  =  0); 
if  Not  Exists  then 

writeln( trm, 'ERROR:  Input  file  does  not  exist'); 

if  Exists  then 
begin 

Asslgn(0utf ile,ParamStr(2)) ; 

{$!-}  Reset(Outfile)  {$!+}; 

Newfile  :=  (’lOresult  <>  0); 
if  Not  Newfile  then 

writeln(trm, 'ERROR:  Output  file  already  exists'); 

if  Newfile  then 
Rewrite(Outf ile) ; 
end; 


end ; 


If  Numfile  and  Exists  and  Newfile  then 
begin 

Init_All;  {  initialize  counters  and  loop  variables} 

Read_Stim;  {  read  in  the  stimuli  and  conditions  from  disc} 


repeat 

Pause_program; 
Stim_Sel ; 
Memory_Set; 
Run_Trlal; 
Stern_Trial ; 
Feed  back; 


{  holds  program  until  key  is  struck} 
{  selects  targ  and  dist  stimuli} 

{  present  Memory  set  } 

{  present  Run  trial  } 

{  present  Stern  trial  } 

{  present  feedback  if  option  set} 


until  done; 


Session_end;  {Write  trials  to  disc  and  terminate  session} 

end; 
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END.  {pmsm.pas} 
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Appendix  B 


CRT. PAS  This  is  the  include  file  for  CRT  routines  which  permit 
synchronization  of  the  display  and  control  of  the  cursor.  The  file  is 
automatically  included  with  the  {$!  CRT. PAS}  command  in  program  PMSM.PAS. 
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procedure  CrtModule  ;  external  'CRT. BIN'  ; 
procedure  CrtSyncEnd  ;  external  CrtModule [ $3]  ; 
procedure  CrtSyncStart  ;  external  CrtModule[$6) 
procedure  CrtCursorOn  ;  external  CrtModule [$9] 
procedure  CrtCursorOff  ;  external  CrtModule [$C] 
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Appendix  C 


TIMER. PAS  This  is  the  timer  procedure  routines  included  with  the  {$I 
TIMER. PAS}  command  In  program  PMSM.PAS.  This  code  enables  a  clock  used  for 
millisecond  timing  of  events. 
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type 

TimeStamp  = 
record 

count , low, high:  Integer 
end  ; 

procedure  TimerModule;  external  'TIMER.BIN'; 
procedure  Timerlnlt  ;  external  TlmerModule[$3]  ; 
procedure  TimerKlll  ;  external  TimerModule [ $6]  ; 

procedure  TlmerStamp(var  stamp:  TimeStamp)  ;  external  TlmerModule[$9]  ; 
function  TimerMilli(stamp:  TimeStamp):  Integer  ;  external  TimerModule [$C] 

function  UnsignedIntegerToReal(u:  Integer):  Real  ; 
var 

x:  Real  ; 

begin  {  UnsignedIntegerToReal  } 
if  u>=0  then  x  :=  u 
else  X  :=  u+65536.0  ; 

UnsignedIntegerToReal  :=  x  ; 
end  {  UnsignedIntegerToReal  }  ; 

function  TimeStampToReal( stamp:  TimeStamp):  Real  ; 
const 

ShiftLeftl6  =  65536.0; 
var 

rLow,rHlgh,rCount:  Real  ; 
total:  Real  ; 

begin  {  TimeStampToReal  } 

rCount  :*  Un8ignedIntegerToReal(stamp. count)  ; 
rLow  :»  UnsignedIntegerToReal(stamp. low)  ; 
rHigh  :=  UnsignedlntegerToReal(stamp.high)  ; 
total  :=■  rHigh  ; 

total  :*■  total*Shif tLef tl6+rLow  ; 
total  :“  total*Shif tLef tl6+rCount  ; 

TimeStampToReal  :=  total  ; 
end  {  TimeStampToReal  }  ; 

procedure  TimeOfDay(var  hour, min, sec, milli:  Integer)  ; 
const 

TimeBase  =  14.31818e+06  ; 
mllliDiv  =  1.193180e+03  ; 
secDiv  ■“  1.193180e+06  ; 
mlnDiv  -  7.159080e+07  ; 
hourDiv  =  4.295448e+09  ; 

ShiftLeftl6  -  65536.0  ; 

var 

stamp:  TimeStamp  ; 
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total:  Real  ; 

begin  {  TimeOfDay  } 

TimerStampC stamp)  ; 
total  :=  TimeStampToRealC stamp ) 
hour  :=  Trunc( total/hourDiv )  ; 
total  :=  total-total*hourDiv  ; 
min  :=  Trunc(total/minDiv)  ; 
total  :=  total-total*minDiv  ; 
sec  :=  Trunc(total/secDiv)  ; 
total  :=  total-total*secDiv  ; 
mill!  :=  Round( total/milliDiv ) 

end  {  TimeOfDay  }  ; 


Appendix  D 


STMGEN.PAS  This  turbo  Pascal  program  generates  the  stimulus  input  files  used 
by  program  PMSM.  The  program  reads  20  four  letter  words  from  file  STM. DAT 
(see  appendix  E),  then  prompts  the  user  for  the  'Yes"  and  'No'  Key 
assignments  (which  are  to  be  entered  on  the  same  line  without  separation, 
e.g.,  Z/  would  specify  the  'Yes'  key  is  Z  and  the  'No'  key  is  /)•  The 
program  will  generate  96  trial  blocks.  Appendix  F  illustrates  the  output  of 
STMGEN.  The  files  generated  by  this  program  can  also  be  used  by  program 
PMSM.FLX,  the  fortran  version  of  PMSM  which  collects  ERP  data  (see  appendix 
I). 
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{  Modification  date:  Jan  16,  1988} 


Program  Stmgen; 

Function  RandomTwixt ( low, high :  Integer):  Integer  ; 

{  Procedure  to  return  a  random  integer  between  the  values  of  low  and  high  } 
beg  in 

RandomTwixt  :=  Raadom(high-low+l )+] ow; 
end ; 


var 

i . 
j  . 
k. 

loop , 
temp , 
itrial, 
rep, 
soa , 
cmvm , 
delay , 
roemsiz , 
stype, 
session , 

subject:  integer; 

tempi , 
yes_button, 
no_button:  char; 

stimuli:  array [ 1 . . 20 , 1 . . 4 ]  of  char; 
trial  :  array [ 1 . . 96 , 1 . . 7 ]  of  integer; 
block  :  array} 1.. 96]  of  integer; 

sub :  string [2] ; 
ses :  string [ 2] ; 

infile  :  text; 
outfile  :  text; 

begin 

write('Enter  subject  number  '); 

read( subj  ect ) ; 

writeln; 

write('Enter  keys  for  yes  and  no  '); 
read(yes_button,no_button) ; 
writeln; 
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{assign  stimulus  Input  file} 

assign(inf ile,'stm.dat') ; 
reset(lnf lie) ; 

{read  words  from  disc} 

for  1  1  to  20  do  {read  In  stimuli} 

begin 

readln(inflle,stimuli[i] ) ; 
end; 

{close  words  file} 
cl ciseCinf  lie) ; 


for  session  :=  1  to  12  do 
begin 

itrial  ;■  0; 
cmvm  :  •»  2 ; 
soa  :»  500; 

for  rep  :=  1  to  8  do 
for  delay  :»  1  to  3  do 
for  memsiz  ;«  1  to  2  do 
for  stype  :■  1  to  2  do 
begin 

itrial  itrlal+1; 
block[itrial]  :=  itrial; 
trial[itrial,l]  :•«  subject; 
trial[itrial,2]  :■  session; 
trial[itrial,3j  soa; 

If  delay  -  1  then  trlal[itrial,4]  :■  0; 

if  delay  »  2  then  trlal[ Itrial, 4]  :■  4; 

if  delay  3  then  trlal[itrlal,4)  15 

trial [itrial ,5]  :■  cmvm; 
trial[itrial,6]  :■  memsiz*2; 
trial[itrial ,7]  ;■  stype; 
end; 

{permute  condition} 

for  k  1  to  40  do 

for  loop  :■  1  to  96  do 
begin 

j  :■  RandomTwlxt( 1 ,96) ; 
temp  :■  block[j]; 
blockij}  block! loop}; 
block[loop]  temp; 
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end ; 


{assign  disc  file} 
str(subj  ect ,sub) ; 
str( session , ses ) ; 

writelnC trm, 'Generating  file;  Stm' , sub, ,ses , ' .dat') ; 
assignC  outf ile stm'+sub+'_'+ses+' .dat' ) ; 
rewrite(outf ile) ; 

{write  stimulus  files  to  disc} 

for  i  :=  1  to  20  do 
begin 

wri teln( outf ile, stimul i [ i ] ) ; 
end ; 

wri teln( outf ile, yes_but ton, no_but ton) ; 

for  i  :=  1  to  96  do 
begin 

for  j  ;=  1  to  7  do  write(outf  ile, trial  [block.[i]  ,j  ]  :5) 
writeln(outfile) ; 
end ; 

close(outf ile) ; 

end ; 
end. 
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Appendix  E 


File  STM. DAT  Input  file  used  by  program  STMGEN  to  generate  stimulus  sequence 
files  for  program  PMSM.  The  file  contains  20  four  letter  words. 
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SIGN 

CLAY 

WISH 

GAME 

FARM 

DEAL 

HOUR 

CLUB 

HAIR 

SHIP 

FILE 

RAIN 

DUST 

TONE 

EDGE 

TERM 

MARK 

BANK 

PuSE 

fcASE 


Appendix  F 


Zii£  STM1--1.DAT  Stimulus  Input  file  used  by  program  PMSM.  This  file  Is  the 
output  of  program  STMGEN.  This  file  contains  20  word  stimuli,  the  response 
assignment  keys  for  Yes'  and  'No'  responses,  and  96  primary/secondary 
memory  trial  parameters.  Each  of  the  96  trials  has  seven  parameters:  They 
are:  subject,  session,  stlmllus  onset  asynchrony  (SOA)  between  the  retrieval 
cue  and  the  probe,  number  of  running  memory  trials  between  the  memory  set 
and  the  probe  (Delay),  consistent  (1)  or  varied  mapping  of  stimulus  to 
response,  memory  set  size,  and  stimulus  response  type  (target*!, 
non-target*2). 
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SIGN 

CLAY 

WISH 

GAME 

FARM 

DEAL 

HOUR 

CLUB 

HAIR 

SHIP 

FILE 

RAIN 

DUST 

TONE 

EDGE 

TERM 

MARK 

BANK 

ROSE 

BASE 

Z/ 

1 

1 

1 

1 

1 

1 

I 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

I 

1 

1 

1 

1 

1 

1 


I  500  4  2  2  1 

1  500  0  2  2  1 

1  500  15  2  4  2 

1  500  15  2  4  2 

1  500  4  2  4  2 

1  500  4  2  2  2 

1  500  15  2  2  1 

1  500  0  2  4  2 

1  500  15  2  4  2 

1  500  15  2  2  1 

1  500  4  2  4  1 

1  500  4  2  4  1 

1  500  0  2  2  1 

1  500  0  2  2  1 

1  500  0  2  4  2 

1  500  0  2  2  2 

1  500  15  2  4  1 

1  500  0  2  4  2 

1  500  0  2  4  2 

1  500  4  2  2  2 

1  500  0  2  2  2 

1  500  4  2  4  1 

1  500  4  2  2  1 

1  500  4  2  2  1 

1  500  0  2  4  1 

1  500  4  2  2  1 

1  500  15  2  4  1 

1  500  4  2  2  1 

1  500  4  2  4  1 

1  500  0  2  2  1 
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1  1  500  4  2  4  2 

1  1  500  4  2  2  2 

i  1  500  0  2  4  1 

1  1  500  0  2  2  1 

1  1  500  4  2  2  2 

1  1  500  4  2  4  2 

1  1  500  0  2  2  2 

1  1  500  4  2  4  2 

1  1  500  4  2  4  2 

1  1  500  15  2  4  2 

1  1  500  4  2  4  I 

1  1  500  15  2  2  2 

1  1  500  0  2  4  2 

1  1  500  4  2  2  2 

1  1  500  0  2  2  2 
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Appendix  G 


AVE . PAS  Program  to  examine  output  file  created  by  PMSM  and  give  quick 
summary  statistics.  This  program  will  need  to  be  customized  for  the 
experimental  conditions. 


{$g512}  {redirect  input  from  disc} 

{$p512}  {redirect  output  to  disc} 

program  ave; 

{Modification  date:  Mar  20,  1987} 


var 

id  :  array [1.. 108]  of  integer; 

ave  :  array  [ 1 . .2 , 1 . .2 , 1 .  .2 , 1 .  .3 , 1 . . 3 , 1 • .2 ]  of  real 

rmem  :  array  {1..2,1..2,1.«3,1»«A]  of  real; 

subject, 

session, 

lines, 

delay, 

soa , 

memsiz , 

cmvm, 

posneg , 

corr , 

point , 

i. 

j:  integer; 
rt:  real; 
begin 

{blank  output  array} 

for  cmvm  :=  1  to  2  do 

for  posneg  :=  1  to  2  do 

for  memsiz  :=  1  to  2  do 

for  soa  :=  1  to  3  do 

for  delay  :=  1  to  3  do 

begin 

ave [ cmvm, posneg , memsiz , soa, delay , 1 ]  0.0; 

ave[cmvm, posneg, memsiz, soa, delay, 2]  0.0; 

rmem[ cmvm, memsiz , delay , 1 ]  :=  0.0; 

rmem[cmvm, memsiz, delay, 2]  :=  0.0; 

rmem[ cmvm, memsiz, delay, 3]  :*  0.0; 

rmem[cmvm, memsiz, delay, 4]  :=  0.0; 

end ; 

{loop  through  data  file  —  increment  appropriate  bins} 

while  not  EOF  do 

begin 

lines  :»1; 

for  j  :=  1  to  108  do 
begin 

read ( id [j  ] ) ; 
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lf(llnes  mod  16)  =0  then  readln; 

lines  :=  lines+1; 

end; 

readln; 

{assign  values  read  from  disc} 


if 

id[4] 

s 

0  then  delay 

I  ^ 

1; 

(1 

s 

delay  0} 

if 

id[4] 

Z5 

4  then  delay 

:= 

2; 

(2 

delay  4} 

if 

id[4] 

= 

15  then  delay 

•  s 

3; 

(3 

= 

delay  15} 

if 

id[5] 

200  then  soa 

*  m 

1; 

(1 

soa  200} 

if 

id[5] 

= 

500  then  soa 

:- 

2; 

(2 

- 

soa  500} 

if 

ld[5] 

= 

1000  then  soa 

:  = 

3; 

(3 

BS 

soa  1000} 

posneg  :=  id [102];  {l”pos,  2-neg} 

memsiz  :=  id[6]  div  2;  {l=*memsi22,  2“mem8iz4} 

cmvm  :=  id[7];  {l=cm,  2“vm} 

rt  ;=  id [101];  (reaction  time} 

corr  :=  id[104];  (correctness:  l=hit,  2"cr,  3“fa,  4“miss} 

subject  :=  id[l];  (subject  number} 

session  :=  id[2];  (session  number} 

case(corr)  of 
1,2: 
begin 

ave[ cmvm, posneg , memsiz ,soa, delay , 1]  :■ 
ave[cmvm,posneg,memsiz,soa,delay, l]+rt; 
ave[ cmvm, posneg ,memsiz ,8oa,delay,2]  :- 
av e [cmvm, posneg, memsiz,soa,delay, 2] +1.0; 
end; 
end; 

for  i  :=  3  to  id[4]  do 
begin 

rt  :*  id[(i*5)+6] ; 
case(id[ (i*5)+9] )  of 
1,2:  begin 

rmem[ cmvm, memsiz, delay, 1] 
rmem[ cmvm,memslz,delay , 1 ]+l .0 ; 

rmem[ cmvm, memsiz , delay, 4] 
rmem[ cmvm, memsiz, d elay, 4 ]+rt; 
end; 

3,4:  rmen[ cmvm, memsiz, delay, 2]  :■ 

rmem[ cmvm, memsiz, d elay, 2] +1 .0; 

0  :  rmem[cmvm, memsiz, delay, 3]  :■ 

rmem[ cmvm, memsiz, delay, 3] +1 .0; 

end; 
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end ; 
end ; 

{calculate  averages  and  write  to  output  file} 

writeln( 'Subj ect  '.subject); 
writeln('Session  '.session); 

for  cmvm  :=  1  to  2  do 
begin 

if  cmvm  =  1  then  writeln('Consistent  Mapping'); 
if  cmvm  =  2  then  writeln( 'Varied  Mapping'); 
for  posneg  :=  1  to  2  do 
begin 

if  posneg  =  1  then  writeln('Positive  trials'); 
if  posneg  =  2  then  writeln( 'Negaitve  trials'); 
writeln; 

wri teln( ' Delay  0  4  15  '); 

for  memsiz  :=  1  to  2  do 
for  soa  :=  1  to  3  do 
begin 

if  memsiz  =  I  then  write('M  =  2:  SOA  =  ',soa:l); 
if  memsiz  =  2  then  write('M  =  4:  SOA  =  ',soa:l); 
for  delay  :=  I  to  3  do 
begin 

if  ave[ cmvm, posneg .memsiz ,soa, delay ,2]  >  0  then 

ave [ cmvm .posneg .memsiz. soa .delay . 1 ]  ;= 

ave[cmvm. posneg  .merisiz  .soa. delay.  1  ]  / 

ave [ cmvm . posneg . memsiz. soa .delay . 2 ] ; 

wri te(ave [ cmvm .posneg .memsiz . soa. delay. 1] :5:0); 

write(ave[ cmvm . posneg .memsiz .soa .delay . 2] :4 :0) ; 

end ; 

writeln; 
end ; 

writeln; 
end ; 

writeln; 

writeln('Running  memory  performance'); 
writeln; 

write( 'Memsiz  Delay  #corr  //err  //nrsp  RT'); 
writeln; 

for  memsiz  :=•  1  to  2  do 
for  delay  :=  2  to  3  do 
begin 

if  rmem[ cmvm .memsiz .delay . 1 ]  >  0  "hen 
rmem[ cmvm.memsiz .delay .4 )  := 
rmera[ cmvm .memsiz .delay .4 ] / 
rmem [ cmvm , mems 1 z . d  el  ay . 1 ) ; 

rmem [ cmvm.memsiz .delay .4 ]  :=  rmem[cmvm. memsiz. delay. 4] ; 
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write( (memsiz*2) ; 6) ; 

if  delay  =  2  then  write(' 

if  delay  =  3  then  write('  15'); 

for  point  :=  1  to  4  do  write(rmem[cmvm,memsiz , delay, point] ;6 ;0) 
writeln; 

end; 
end ; 
end. 
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Appendix  H 


PMSM. FLX  Source  code  for  primary/secondary  memory  experiment  used  to 
collect  ERPs.  The  file  is  written  using  a  flees  preprocessor.  See  appendix 
I  for  the  Fortran  Source  Code. 

The  experimental  procedure  is  identical  to  the  Pascal  version  of  the 
program.  After  compilation,  to  execute  the  program  type  "RU  PMSM".  The 
program  prompts  the  user  for  a  stimulus  input  file  (which  should  be  in 
identical  format  to  appendix  F).  A  menu  will  be  presented  with  the 
following  options; 

(R)un  (T)ape  (C)alib  (Q)ult.  "R"  initiates  the  experiment.  "Q"  Quits  the 
program.  "C"  gathers  5  calibration  trials.  "T"  enters  the  tape  manipulation 
menu.  The  tape  menu  consists  of  the  following  options:  (E)xit  the  menu  and 
return  to  main  menu,  (R)ewind  the  tape,  (L)abel  the  tape  header,  (S)klp  a 
specified  number  of  EOFs,  (V)erify  the  data  written  to  tape,  (W)rite  an  EOF. 

Two  special  real-time  options  are  provided:  If  the  user  types  "P"  the 
program  will  pause  until  the  user  presses  carriage  return.  If  the  user 
types  "Q",  the  run-time  program  will  abort  and  return  to  the  main  menu. 

The  program  assumes  5  channels  of  data  are  to  be  recorded.  The  channels  are 
1)  EOC,  2)  Fz,  3)  Cz,  4)  Pz,  and  5)  stimulus  marker.  Each  trial  block  is 
written  to  tape  with  108  identification  parameters  and  5  multiplexed  data 
channels. 
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PMSM.FLX 


C 


MODIFICATION  DATE:  JAN  16,  1988 

PROGRAM  TO  CONTROL  PRIMARY/ SECONDARY  MEMORY  STUDY 

IDS  OUTPUT  BY  THE  PROGRAM: 

ID(1)  -  SUBJECT  NUMBER 

ID(2)  -  SESSION 

ID(3)  -  BLOCK 

ID(4)  -  DELAY 

ID(5)  -  SOA  OF  RETRIEVAL  CUE 

ID(6)  -  SIZE  OF  MEMORY  SET 

ID(7)  -CONSTANT  (1)  OR  VARIED  (2)  MAPPING 
ID(8)  -  CODE  FOR  YES  BUTTON 

ID(9)  -  CODE  FOR  NO  BUTTON 

RUNNING  MEMORY  IDS 

TRIAL  1 

ID(IO)  -  RUNNING  MEMORY  STIMULUS 
ID(ll)  -  RT  (FOR  RUNNING  MEMORY  TRIAL) 

ID(12)  -  STIM  TYPE  FOR  RUNNING  MEMORY  (1=TARGET,  2=NONTARGET) 

ID(13)  -  RESPONSE  FOR  RUNNING  MEMORY  (0=NONE,  1=TARGET,  2-NONTARGET) 
ID(14)  -  CORRECTNESS  FOR  RUNNING  MEMORY:  0=NO  RESPONSE,  1=HIT, 
2-CORRECT  REJECTION,  3=FALSE  ALARM,  4-MISS 

TRIAL  2 

C  ID(15)  -  RUNNING  MEMORY  STIMULUS 

C  ID(16)  -  RT  (FOR  RUNNING  MEMORY  TRIAL) 

C  ID(17)  -  STIM  TYPE  FOR  RUNNING  MEMORY  (1-TARGET,  2-NONTARGET) 

C  ID(18)  -  RESPONSE  FOR  RUNNING  MEMORY  (O-NONE,  1-TARGET,  2-NONTARGET) 

C  ID(19)  -  CORRECTNESS  FOR  RUNNING  MEMORY:  O-NO  RESPONSE,  1-HIT, 

C  2-CORRECT  REJECTION,  3-FALSE  ALARM,  4-MISS 

C 

C  TRIAL  3 
C 

C  ID(20)  -  RUNNING  MEMORY  STIMULUS 

C  ID(21)  -  RT  (FOR  RUNNING  MEMORY  TRIAL) 

C  ID(22)  -  STIM  TYPE  FOR  RUNNING  MEMORY  (1-TARGET,  2-NONTARGET) 

C  ID(23)  -  RESPONSE  FOR  RUNNING  MEMORY  (O-NONE,  1-TARGET,  2-NONTARGET) 

C  ID(24)  -  CORRECTNESS  FOR  RUNNING  MEMORY:  O-NO  RESPONSE,  1-HIT, 

C  2-CORRECT  REJECTION,  3-FALSE  ALARM,  4-MISS 

C 

C  TRIAL  4 

C 

C  ID(25)  -  RUNNING  MEMORY  STIMULUS 
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C  ID(26)  -  RT  (FOR  RUNNING  MEMORY  TRIAL) 

C  ID(27)  -  STIM  TYPE  FOR  RUNNING  MEMORY  (1=TARGET,  2=N0NTARGET) 

C  ID(28)  -  RESPONSE  FOR  RUNNING  MEMORY  (0=NONE,  1=TARGET,  2=N0NTARGET) 

C  ID(29)  -  CORRECTNESS  FOR  RUNNING  MEMORY:  0=NO  RESPONSE,  1=HIT, 

C  2=CORRECT  REJECTION,  3=FALSE  ALARM,  4=MISS 

TRIAL  5 

ID(30)  -  RUNNING  MEMORY  STIMULUS 
ID(31)  -  RT  (FOR  RUNNING  MEMORY  TRIAL) 

ID(32)  -  STIM  TYPE  FOR  RUNNING  MEMORY  (1=TARGET,  2=NONTARGET) 

ID(33)  -  RESPONSE  FOR  RUNNING  MEMORY  (0=N0NE,  1=TARGET,  2=NONTARGET) 
ID(34)  -  CORRECTNESS  FOR  RUNNING  MEMORY:  0=NO  RESPONSE,  1=HIT, 
2=CORRECT  REJECTION,  3=FALSE  ALARM,  4=MISS 

TRIAL  6 
C 

C  ID(35)  -  RUNNING  MEMORY  STIMULUS 

C  ID(36)  -  RT  (FOR  RUNNING  MEMORY  TRIAL) 

C  ID(37)  -  STIM  TYPE  FOR  RUNNING  MEMORY  (1=TARGET,  2-NONTARGET) 

C  ID(38)  -  RESPONSE  FOR  RUNNING  MEMORY  (0=NONE,  1=TARGET,  2=NONTARGET) 

C  ID(39)  -  CORRECTNESS  FOR  RUNNING  MEMORY:  0=NO  RESPONSE,  1=HIT, 

C  2=CORRECT  REJECTION,  3=FALSE  ALARM,  4=MISS 

C 

C  TRIAL  7 

C 

C  ID(40)  -  RUNNING  MEMORY  STIMULUS 

C  ID(41)  -  RT  (FOR  RUNNING  MEMORY  TRIAL) 

C  ID(42)  -  STIM  TYPE  FOR  RUNNING  MEMORY  (1=TARGET,  2-NONTARGET) 

C  ID(43)  -  RESPONSE  FOR  RUNNING  MEMORY  (0=NONE,  1-TARGET,  2-NONTARGET) 

C  ID(44)  -  CORRECTNESS  FOR  RUNNING  MEMORY:  0=N0  RESPONSE,  1-HIT, 

C  2-CORRECT  REJECTION,  3-FALSE  ALARM,  4-MISS 

C 

C  TRIAL  8 

C 

C  ID(45)  -  RUNNING  MEMORY  STIMULUS 

C  ID(46)  -  RT  (FOR  RUNNING  MEMORY  TRIAL) 

C  ID(47)  -  STIM  TYPE  FOR  RUNNING  MEMORY  (1-TARGET,  2-NONTARGET) 

C  ID(48)  -  RESPONSE  FOR  RUNNING  MEMORY  (O-NONE,  1-TARGET,  2-NONTARGET) 

C  ID(49)  -  CORRECTNESS  FOR  RUNNING  MEMORY:  O-NO  RESPONSE,  1-HIT, 

C  2-CORRECT  REJECTION,  3-FALSE  ALARM,  4-MISS 

C 

C  TRIAL  9 

C 

C  ID(50)  -  RUNNING  MEMORY  STIMULUS 

C  ID(51)  -  RT  (FOR  RUNNING  MEMORY  TRIAL) 

C  ID(52)  -  STIM  TYPE  FOR  RUNNING  MEMORY  (1-TARGET,  2-NONTARGET) 

C  ID(53)  -  RESPONSE  FOR  RUNNING  MEMORY  (O-NONE,  1-TARGET,  2-NONTARGET) 

C  ID(54)  -  CORRECTNESS  FOR  RUNNING  MEMORY:  O-NO  RESPONSE,  1-HIT, 

C  2-CORRECT  REJECTION,  3-FALSE  ALARM,  4-MISS 

C 
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TRIAL  10 

C  ID(55)  -  RUNNING  MEMORY  STIMULUS 

C  ID(56)  -  RT  (FOR  RUNNING  MEMORY  TRIAL) 

C  ID(57)  -  STIM  TYPE  FOR  RUNNING  MEMORY  (1=TARGET,  2=N0NTARGET) 

C  ID(58)  -  RESPONSE  FOR  RUNNING  MEMORY  (0=NONE,  1-TARGET,  2-NONTARGET) 

C  ID(59)  -  CORRECTNESS  FOR  RUNNING  MEMORY:  O-NO  RESPONSE,  1-HIT, 

C  2-CORRECT  REJECTION,  3-FALSE  ALARM,  4-MISS 

C 

C  TRIAL  1 1 

C 

C  ID(60) 

C  ID(61) 

C  ID(62) 

C  ID(63) 

C  ID(64) 

C 
C 

C  TRIAL  12 

C 

C  ID(65) 

C  ID(66) 

C  ID(67) 

C  ID(68) 

C  ID(69) 

C 
C 

C  TRIAL  13 

C 

C  ID(70) 

C  ID(71) 

C  ID(72) 

C  ID(73) 

C  ID(74) 

C 
C 

C  TRIAL  14 

C 

C  ID(75) 

C  ID(76) 

C  ID(77) 

C  ID(78) 

C  ID(79) 

C 
C 

C  TRIAL  15 

C 

C  ID(80)  -  RUNNING  MEMORY  STIMULUS 

C  ID(81)  -  RT  (FOR  RUNNING  MEMORY  TRIAL) 

C  ID(82)  -  STIM  TYPE  FOR  RUNNING  MEMORY  (1-TARGET,  2-NONTARGET) 

C  ID(83)  -  RESPONSE  FOR  RUNNING  MEMORY  (O-NONE,  1-TARGET,  2-NONTARGET) 


-  RUNNING  MEMORY  STIMULUS 

-  RT  (FOR  RUNNING  MEMORY  TRIAL) 

-  STIM  TYPE  FOR  RUNNING  MEMORY  (1-TARGET,  2-NONTARGET) 

-  RESPONSE  FOR  RUNNING  MEMORY  (O-NONE,  1-TARGET,  2-NONTARGET) 

-  CORRECTNESS  FOR  RUNNING  MEMORY:  O-NO  RESPONSE,  1-HIT, 
2-CORRECT  REJECTION,  3-FALSE  ALARM,  4-MISS 


-  RUNNING  MEMORY  STIMULUS 

-  RT  (FOR  RUNNING  MEMORY  TRIAL) 

-  STIM  TYPE  FOR  RUNNING  MEMORY  (1-TARGET,  2-NONTARGET) 

-  RESPONSE  FOR  RUNNING  MEMORY  (O-NONE,  1-TARGET,  2-NONTARGET) 

-  CORRECTNESS  FOR  RUNNING  MEMORY:  O-NO  RESPONSE,  1-HIT, 
2-CORRECT  REJECTION,  3-FALSE  ALARM,  4-MISS 


-  RUNNING  MEMORY  STIMULUS 

-  RT  (FOR  RUNNING  MEMORY  TRIAL) 

-  STIM  TYPE  FOR  RUNNING  MEMORY  (1-TARGET,  2-NONTARGET) 

-  RESPONSE  FOR  RUNNING  MEMORY  (O-NONE,  1-TARGET,  2-NONTARGET) 

-  CORRECTNESS  FOR  RUNNING  MEMORY:  O-NO  RESPONSE,  1-HIT, 
2-CORRECT  REJECTION,  3-FALSE  ALARM,  4-MISS 


-  RUNNING  MEMORY  STIMULUS 

-  RT  (FOR  RUNNING  MEMORY  TRIAL) 

-  STIM  TYPE  FOR  RUNNING  MEMORY  (1-TARGET,  2-NONTARGET) 

-  RESPONSE  FOR  RUNNING  MEMORY  (O-NONE,  1-TARGET,  2-NONTARGET) 

-  CORRECTNESS  FOR  RUNNING  MEMORY:  O-NO  RESPONSE,  1-HIT, 
2-CORRECT  REJECTION,  3-FALSE  ALARM,  4-MISS 
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c 

c 

c 

c 

c 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c- 

c 


ID(84)  -  CORRECTNESS  FOR  RUNNING  MEMORY:  O-NO  RESPONSE,  1-HIT, 
2-CORRECT  REJECTION,  3-FALSE  ALARM,  4-MISS 

STERNBERG  IDS 


c 

ID(85) 

- 

MEMORY 

SET 

WORD 

n. 

LETTER 

#1 

c 

ID(86) 

- 

MEMORY 

SET 

WORD 

n. 

LETTER 

in 

c 

ID(87) 

- 

MEMORY 

SET 

WORD 

//i. 

LETTER 

in 

c 

1D(88) 

- 

MEMORY 

SET 

WORD 

#1, 

LETTER 

//4 

c 

c 

ID(89) 

_ 

MEMORY 

SET 

WORD 

//2, 

LETTER 

in 

c 

ID(90) 

- 

MEMORY 

SET 

WORD 

(/2, 

LETTER 

in 

c 

ID(91) 

- 

MEMORY 

SET 

WORD 

#2, 

LETTER 

a  3 

c 

r* 

ID(92) 

- 

MEMORY 

SET 

WORD 

n. 

LETTER 

iiU 

C 

c 

ID(93) 

— 

MEMORY 

SET 

WORD 

#3. 

LETTER 

ii\ 

c 

ID(94) 

- 

MEMORY 

SET 

WORD 

LETTER 

in 

c 

ID(95) 

- 

MEMORY 

SET 

WORD 

//3, 

LETTER 

in 

c 

fy 

ID(96) 

- 

MEMORY 

SET 

WORD 

//3, 

LETTER 

#4 

L 

c 

ID(97) 

MEMORY 

SET 

WORD 

#4, 

LETTER 

n 

c 

ID(98) 

- 

MEMORY 

SET 

WORD 

H, 

LETTER 

in 

c 

ID(99) 

- 

MEMORY 

SET 

WORD 

#4, 

LETTER 

in 

ID(IOO)  -  MEMORY  SET  WORD  #4,  LETTER  //4 

ID(lOl)  -  RT  (FOR  STERNBERG  TRIAL) 

ID(102)  -  STIM  TYPE  (1-TARGET,  2-NONTARGET) 

ID(103)  -  RESPONSE  (0=^NONE,  1-TARGET,  2-NONTARGET) 

ID(104)  -  CORRECTNESS:  O-NO  RESPONSE,  1-HIT,  2-CORRECT  REJECTION, 

3-FALSE  ALARM,  4-MISS 

ID(105)  -  STIMULUS  PRESENTED,  LETTER  //I 
IDU06)  -  STIMULUS  PRESENTED,  LETTER  #2 
ID(107)  -  STIMULUS  PRESENTED,  LETTER  //3 
ID(108)  -  STIMULUS  PRESENTED,  LETTER  #4 


C 

IMPLICIT  INTEGER  (A-Z) 

INTEGER  COND( 96,7), LSTIM( 3 ) , ISEED( 2 ) 

BYTE  NAME( 17 ), ERR, TEMPI ,WORD(5) ,NUM(9) ,DIGIT(2) , YES, NO, 
+CMTARG( 5,4), CMDI ST( 5 , 4 ) , VMSTIM( 10,4) ,TARG( 5 , 4 ) , DI ST( 5 , 4 ) 
REAL  SEED, RAN 

EQUIVALENCE  (SEED, ISEED) , ( ISEED( 1 ) ,S1 ) , ( ISEED( 2) ,S2 ) 

COMMON  /RAN/  S1,S2 

COMMON  /DATA/  ID( 108) ,DATA( 5 ,210) 

C 

DATA  STMDUR,NBACK/200,2/ 

DATA  NCHAN ,NPTS ,DRATE , NIDS , BASE/ 5 , 2 10 , 10 , 108 , 100/ 

DATA  NUM/'1','2','3','4','5','6','7','8','9'/ 
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GET  A  RANDOM  SEED  AND  WARM  UP  RANDOM  NUMBER  GENERATOR 

SEED=SECNDS(0) 

DO  (RANDOM=l ,2000) 

CALL  RAN(S1,S2) 

FIN 

INITIALIZE  MATROX 
CALL  MCINIT(l) 

SET  UP  DISPLAY  ASSIGNMENTS 
CALL  GQINIT 

CALL  GQTTL('  PM/SM  EXPERIMENT') 

CALL  GQLABL(1,'FZ  ELECTRODE') 

CALL  GQLABL(2,'CZ  ELECTRODE') 

CALL  GQLABL(3,'PZ  ELECTRODE') 

CALL  GQLABL(4,'EOG  ELECTRODE') 

CALL  GQLABL(5,'STM  MARKER') 

ASSIGN  STIMULUS  FILE  FROM  DISC 

TYPE  INPUT  STIMULUS  FILE  FOR  THIS  SUBJECT' 

CALL  GETSTR(5,NAME,16,ERR) 

CALL  ASSIGN(40,NAME,0,'0LD') 

READ  STIMULUS  FILE 

DO  (1=1,5) 

READ(40,1)  (CMTARG(I,J),J=1,4) 

FIN 

DO  (1=1,5) 

READ(40,1)  (CMDIST(I,J),J=1,4) 

FIN 

DO  (1-1,10) 

READ(40,1)  (VMSTIM(I,J),J-1,4) 

FIN 

READ(40,2)  YES, NO 
IF(YES.EQ.'/') 

YES BUT- 1 

NOBUT-0 

FIN 

IF(YES.EQ.'Z') 

YESBUT-0 

NOBUT-1 

FIN 
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DO  (J=l,96) 

READ(40,3)  (C0ND(J,I),I=1,7) 
FIN 


C 

1  FORMAT(4Al) 

2  F0RMAT(2A1) 

3  FORMAT(7I5) 

C 

C  ===============  get  the  option  FROM  THE  USER  ================ 

C 

C  GET  THE  OPTION 

C 

100  CONTINUE 

TYPE  110 

no  FORMAT(//'$  MENU:  (R)UN,  (T)APE,  (C)ALIB,  (Q)UIT->') 

ACCEPT  120,I0PT 
120  FORMAT(Al) 

C 

C  JUMP  ON  OPTION 

C 

IF(IOPT.EQ.'C')  CALL  CALIB(NCHAN,NPTS,DRATE,NIDS) 
IF(IOPT.EQ.'Q')  GOTO  2000 
IF(IOPT.EQ.'R')  GOTO  500 

IF(IOPT.EQ.'T')  CALL  TAPE(NCHAN,NPTS,NIDS) 

GOTO  100 
C 

C  ==========================  RUN  TIME  SETUP  ====================== 

C 

C  TRAP  "Q"  FOR  ABORT. 

C 

500  ISWFLG=0 

CALL  SCCA(ISWFLG) 

CALL  KEYON  ITURN  OFF  ECHO,  CHECK  FOR  KEY  PRESSES 

C 

C  INITIALIZE  NBYTES  FOR  TAPE  WRITE 

C 

NBYTES=(NIDS+NCHAN*NPTS)*2 

C 

C  INITIALIZE  DEVICES 

C 

CALL  LINOFF 
CALL  CLINIT 
CALL  DIINIT 
CALL  DO IN IT 
DO  (ICLOCK=2,5) 

CLOCK* I CLOCK 
CALL  PCINITC CLOCK) 

FIN 

C 

C==-==========  tlAIN  RUN  TIME  LOOP-=— ===-=«-==»== 
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C 

DO  (BLOCK=l,96) 

ZERO  ID  ARRAY 
CALL  ZER0(ID,NIDS*2) 

CHECK  FOR  "Q"  AND  "P" 

CALL  KEYS(KEYFLG)  ! CHECK  FOR  KEY  PRESSES 

IF(KEYFLG.EQ.l) 

CALL  CLKILL 
CALL  MCERAS 
CALL  DIKILL 
CALL  DOKILL 
CALL  ADKILL 
CALL  KEYOFF 
DO  (ICLOCK=2,5) 

CLOCK* I CLOCK 
CALL  PCKILL( CLOCK) 

FIN 

CALL  EOF(O) 

CALL  EOF(O) 

CALL  SKIP(0,-2) 

GOTO  100 
FIN 


IBLK=BLOCK 

IF(((IBLK/16)*16).EQ.IBLK.AND.IBLK.NE.96) 

TYPE  Program  paused' 

CALL  MCTEXT( 20 ,200 Pause  —  Press  button  to  continue  ') 
CALL  MCDISP(l) 

WAIT  FOR  SUBJECT  TO  PRESS  BUTTON 


C 

C 

C 


CALL  DISET(0,' SINGLE', ISTATO) 

CALL  DISETd, 'SINGLE', ISTATl) 

CALL  WAIT (ISTATO, I ST ATI) 

CALL  MCDISP(O) 

CALL  MCERAS 
FIN 

READ  CONDITION,  STORE  IDS  AND  SELECT  STIMULI 

SUB  *C0ND(BL0CK,1) 

SESS  -COND(BLOCK,2) 

SOA  =C0ND(BL0CK,3) 

DELAY  -COND(BLOCK,4) 

CMVM  *COND(BLOCK,5) 

MEMSIZ-COND( BLOCK , 6 ) 

STYPE  -COND(BLOCK,7) 
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ID(1)=  SUB 
II><2)=  SESS 
ID(3)=  BLOCK 
ID(4)=  DELAY 
ID(5)=  SOA 
ID(6)=  MEMSIZ 
ID(7)=  CMVM 
ID(8)=  YESBUT 
ID(9)=  NOBUT 

INFORM  EXPERIMENTER  OF  CONDI TLC-N 
TYPE  30 , BLOCK , DELAY , CMVM , MEMS I Z , ST YPE 

FORMAT($'  Block  ',13,'  Delay  ',12,'  CMVM  ',11,'  Memsiz  ',Ii, 
+'  PosNeg  ',11) 

DO  (K=l,10) 

PERMUTE  CM  TARGETS 

DO  (LOOP  =1,5) 

J=ITWIXT(1 ,5) 

DO  (1=1,4) 

TEMPI  =  CMTARG(J,I) 

CMTARG(J,I)  =  CMTARG(L00P,1) 

CMTARG(LOOP,I)  =  TEMPI 

FIN 

FIN 

PERMUTE  CM  DISTRACTORS 

DO  (LOOP  =1,5) 

J=ITWIXT( 1,5) 

DO  (I=i;4) 

TEMPI  =  CMDIST(J,I) 

CMDIST(J,I)  =  CMDIST(LOOP,I) 

CMDIST(LOOP,I)  =  TEMPI 

FIN 

FIN 

PERMUTE  VM  STIMULI 

DO  (LOOP  =1,10) 

J=ITWIXT(1,10) 

DO  (1=1,4) 

TEMPI  =  VMSTIM(J,I) 

VMSTIM(J,I)  =  VMSTIM(LOOP,I) 

VMSTIM(LOOP,I)  =  TEMPI 

FIN 

FIN 
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FIN 


C 

CHOOSE  CM  STIMULI 

IF(CMVM.EQ.l) 

DO  (J=1,MEMSIZ) 

DO  (1=1,4) 

TARG(J,I)=CMTARG(J,I) 
FIN 
FIN 

DO(J=l ,5) 

DO(I=l ,4) 

DISK J,I)=CMDIST( J,I) 
FIN 
FIN 

FIN 

CHOOSE  VM  STIMULI 


IF(CMVM.EQ.2) 

DO  (J=1,MEMSIZ) 

DO  (1=1,4) 

TARG(J,I)=VMSTIM(J,I) 

FIN 
FIN 
C 

DO( J=MEMSIZ+ 1 ,MEMSIZ+5 ) 

DO(I=l ,4) 

DIST( J-MEMSIZ,I)=VMSTIM(J,I) 

FIN 

FIN 

C 

FIN 

C 

C  ============1======  present  MEMORY  SET  »== 

C 

C 

C  WRITE  STIMULI  TO  MATROX  AND  STORE  IDS 
C 

DO(I=l ,MEMSIZ) 

C 

IF(MEMSIZ.EQ.2) 

XPOS=60+(I*35) 

YPOS-lOO 

FIN 

C 

IF(MEMSIZ.EQ.4) 

XPOS-27+(I*35) 
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YP0S=100 

FIN 

C 

DO  (J=l,4) 

WORD(J)=TARG(I ,J) 

FIN 

C 

D0( J=1  ,4) 

ID(80+(4*I)+J)=TARG(I,J) 

FIN 

WORD(5)=0 

CALL  MCTEXT(XPOS,YPOS,WORD) 

FIN 

C 

C  TURN  ON  STIMULUS 
C 

CALL  MCSYNC  !  SUNCH  DISPLAY 

CALL  MCDISP(l)  !  TURN  ON  STIMULUS 

CALL  PCSET( 1+1 , 'SINGLE', 3, 3000, IFLAGl) 

CALL  WAIT( IFLAGl)  !  WAIT  FOR  3  SECONDS 
CALL  MCDISP(O)  !  TURN  OFF  STIMULUS 

CALL  MCERAS  !  ERASE  MATROX 

WAIT  FOR  1500  MSEC  BEFORE  PRESENTING  PROBES 

CALL  PCSET( 1+ 1 , 'SINGLE' ,3,1500 , IFLAGl ) 

CALL  WAIT( IFLAGl) 

===============  running  memory  TRIALS  ==========. 

IF(DELAY.NE.O) 

DO  (ITRL=1,IABS( DELAY)) 

PREPARE  STIMULI  DURING  PRE-STIMULUS  INTERVAL 
C 

CALL  PCSET( 1+1 ,' SINGLE', 3, 5, IFLAGl) 
POSNEG=ITWIXT( 1,3) 

IF(POSNEG.EQ.3)  P0SNEG=2 
C 

C  IF  A  POSITIVE  TRIAL,  REPEAT  STIMULUS  NBACK 

C  IF  A  NEGATIVE  TRIAL,  SELECT  A  NEW  STIMULUS 

C 

IF(POSNEG.EQ.l)  !  POSITIVE  TRIAL 
IFdTRL.LE. NBACK)  JSTIM=ITWIXT(  1  ,9 ) 
IFdTRL.GT. NBACK)  JSTIM=LSTIM( NBACK) 

FIN 

C 

IF(POSNEG.EQ.2)  !  NEGATIVE  TRIAL 
10  JSTIM-ITWIXTd  ,9) 

IF( JSTIM.EQ.LSTIM(NBACK))  GOTO  10 
FIN 
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c 

C  STORE  SUCCESSIVE  STIMULI 
C 

LSTIM(3)  =  LSTIM(2) 

LSTIM(2)  =  LSTIM(l) 

LSTIM(l)  =  JSTIM 
ID((ITRL*5)+5)  =  JSTIM 
C 

XPOS=120 

YPOS=120 

DIGIT(1)=NUM( JSTIM) 

DIGIT(2)=0 

IF  ( DELAY. GT.O)  CALL  MCTEXT(XPOS,YPOS, DIGIT) 

C 

CALL  WAIT(IFLAGl)  !  WAIT  OUT  PRESTIMULUS  INTERVAL 
C 

CALL  MCSYNC  !  SYNCH  DISPLAY 

CALL  MCDISP(l)  !  TURN  ON  STIMULUS 

CALL  PCSET( 1+1 .'SINGLE', 3, lOOO.IFLAGl) 

CALL  DIRESP(RESP, 'SINGLE', 2, 1+1 ,RT,RSTAT) 

CALL  PCSET( 1+2 , 'SINGLE' ,3,200, IFLAG2 ) 

CALL  WAIT(IFLAG2) 

CALL  MCDISP(O)  !  TURN  OFF  STIMULUS 

CALL  MCERAS  !  ERASE  MATROX 

CALL  WAIT(IFLAGl)  !WAIT  FOR  END  OF  1000  MSEC  INTERVAL 
C 

C  CLASSIFY  TRIAL  AND  STORE  RT  AND  STIM  TYPE  (1=HIT,  2<=CR,  3=FA,  4*MISS) 
C 

ID((ITRL*5)+6)  =  RT 
ID((ITRL*5)+7)  =  POSNEG 
ID((ITRL*5)+8)  =  0 
ID((ITRL*5)+9)  =  0 
C 

IF(RSTAT.NE.O) 

C 

IF(RESP.EQ.YESBUT) 

ID((ITRL*5)+8)  -  1 

IF  (ID((ITRL*5)+7).EQ.l)  ID((ITRL*5)+9)  =  1 
IF  (ID((ITRL*5)+7).EQ,2)  ID((ITRL*5)+9)  -  3 
FIN 
C 

IF(RESP.EQ.NOBUT) 

ID((ITRL*5)+8)  =  2 

IF  (ID((ITRL*5)+7).EQ.2)  ID((ITRL*5)+9)  =  2 
IF  (ID((ITRL*5)+7).EQ.l)  ID((ITRL*5)+9)  =•  4 
FIN 
C 

FIN 

C 

FIN 

FIN 
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C 

C  ======»==*=.====  STERNBERG  PROBE  TRIAL  ——————— 

C 

C 

C  START  THE  DIGITIZER 
C 

CALL  ADSET( DATA, 2 , 1 .NCHAN.NPTS , 3 ,DRATE , DSTAT) 

wait  out  baseline  clock 

CALL  PCSET( 1+ 1 , ' SINGLE' , 3 , BASE , IFLAGl ) 

CALL  WAIT( IFLAGl) 

DELIVER  RETRIEVAL  CUE 

IF(SOA.GE.200) 

XP0S=120 
YPOS=120 

CALL  MCCHAR(XPOS,YPOS,'*') 

CALL  MCSYNC  !  SYNCH  DISPLAY 

CALL  DOSET(O)  !TURN  ON  STIMULUS  MARKER 

CALL  MCDISP(l)  !  TURN  ON  MATROX 

CALL  PCSET( 1+ 1 , 'SINGLE' , 3 , SOA , IFLAG 1 ) 

CALL  PCSET( 1+2 ,' S INGLE' ,3,200, IFLAG2 ) 

CALL  WAIT(IFLAG2) 

CALL  MCDISP(O)  !  TURN  OFF  MATROX 

CALL  DOCLR(O)  !  TURN  OFF  STIMULUS  MARKER 

CALL  MCERAS  !  ERASE  MATROX 

CALL  WAIT( IFLAGl)  !  WAIT  FOR  END  OF  SOA  INTERVAL 
FIN 

PREPARE  STIMULI  DURING  PRE-STIMULUS  INTERVAL 

CALL  PCSET( 1+ 1 , 'SINGLE ',3,5, IFLAGl ) 

XPOS=112 
YP0S=140 
C 

IF(STYPE.EQ.l)  ! POSITIVE  TRIAL 
I=ITWIXT(1,MEMSIZ) 

DO(J-l,4) 

WORD(J)=TARG(I,J) 

ID(104+J)=TARG(I,J) 

FIN 

FIN 

C 

IF  (STYPE.EQ.2)  'NEGATIVE  TRIAL 
I-ITWIXT(1,5) 

DO(J-l,4) 

WORD(J)-DIST(I,J) 

ID(104+J)-DIST(I,J) 

FIN 
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FIN 

WORD(5)=0 

CALL  MCTEXT(XPOS,YPOS,WORD) 

CALL  WAIT(IFLAGl) 

CALL  MCSYNC  !  SYNCH  DISPLAY 

CALL  DOSET(O)  !  TURN  ON  STIMULUS  MARKER 

CALL  MCDISP(l)  !  TURN  ON  DISPLAY 

CALL  PCSET( 1+1, 'SINGLE', 3, 3000»IFLAG1) 

CALL  DIRESP(RESP,' SINGLE' ,2 ,1+1 ,RT,RSTAT) 

CALL  PCSET( 1+2 , ' SINGLE' , 3 , 200 , IFLAG2 ) 

CALL  WAIT(IFLAG2) 

CALL  MCDISP(O)  !  TURN  OFF  DISPLAY 

CALL  DOCLR(O)  !  TURN  OFF  STIMULUS  MARKER 

CALL  MCERAS  !  CLEAR  MARTOX 

CALL  WAIT(IFLAGl)  !  WAIT  FOR  END  OF  3  SECOND  INTERVAL 

CLASSIFY  TRIAL  AND  STORE  RT  AND  STIM  TYPE  (1-HIT,  2-CR,  3-FA,  4-MISS) 

ID(lOl)  -  RT 
ID(102)  =  STYPE 
ID(103)  =  0 
ID(104)  -  0 

IF(RSTAT,NE.O) 

IF(RESP.EQ.YESBUT) 

ID(103)  -  1 

IF  (ID(102).EQ.l)  ID(104)  -  1 

IF  (ID(102).EQ.2)  ID(104)  -  3 

FIN 

IF(RESP.EQ.NOBUT) 

ID(103)  -  2 

IF  (ID(102).EQ.2)  ID(104)  -  2 

IF  (ID(102).EQ.l)  ID(104)  -  4 

FIN 

FIN 

CHECK  IF  ANY  POINTS  ARE  OUT  OF  BOUNDS 
ERROR-0 

DO  (IPTS-1,210) 

DO  (ICHAN-1,4) 

IF( IABS( DATA( ICHAN, IPTS ) ) .GE, 2047 )  ERROR-ERROR+1 

FIN 

FIN 

WRITE  DATA  TO  TAPE 

CALL  PUT(0,ID,NBYTES,MSTAT) 
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TYPE  20,  ID(101),ID(104) 

FORMATC'  RT  ',14,'  Corr  ',11) 

IF  (ERROR. NE.O)  TYPE  *, 'Error  in  last  trial  of  digitized  EEC' 
CALL  DISP(NCHAN,NPTS) 

CALL  WAIT(MSTAT) 

C 

FIN 

C 

C  KILL  THE  DEVICES 
C 

CALL  CLKILL 
CALL  MCERAS 
CALL  DIKILL 
CALL  DOKILL 
CALL  ADKILL 
CALL  KEYOFF 
DO  (ICLOCK=2,5) 

CL0CK=ICL0CK 
CALL  PCKILLC CLOCK) 

FIN 

C 

C  WRITE  DOUBLE  EOF  AND  SKIP  BACK  I 
C 

CALL  EOF(O) 

CALL  EOF(O) 

CALL  SKIP(0,-2) 

C 

C  RETURN  TO  MAIN  MENU 
C 

GOTO  100 
C 

C  «===========a=====  QUIT  SECTION 

C 

C  TURN  OFF  THE  DISPLAY  AND  EXIT  PROGRAM 
C 

2000  CONTINUE 

CALL  CLKILL 
CALL  MCERAS 
CALL  GQINIT 
CALL  DIKILL 
CALL  ADKILL 
CALL  KEYOFF 
STOP 
END 
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C  * 

C  SUBROUTINES  * 

C  * 

CA4tA4cAVt*4t:fcAAAA4tA4cAA4t*AAAA4r**:kA4iAA4EA*A*AAAAA*:kilt***4tAAAA4t4t4t*4tA4tAA 

CCC 

CCC  SUBROUITNE  CALIB 

CCC 

C 

C  GATHERS  FIVE  CALIBRATION  TRIALS 
C 

SUBROUTINE  CALIB( NCHAN ,NPTS ,DRATE , NIDS ) 

IMPLICIT  INTEGER  (A-Z) 

C 

COMMON  /DATA/  ID( 108) ,DATA(5,210) 

NBYTES-(NIDS+NCHAN*NPTS)*2 

C 

CALL  LINOFF 
CALL  CLINIT 
WRITE(7,10) 

10  FORMAT(/////'  TURN  ON  CALIBRATION  PULSE') 

PAUSE  'HIT  <CR>  TO  CONTINUE' 

C  DIGITIZE  5  TRIALS  OF  CALIBRATION  DATA 
DO  (LOOP-1,5) 

CALL  ADSET( DATA, 2 , 1 , NCHAN ,NPTS , 3 ,DRATE , DSTAT) 

CALL  WAIT( DSTAT)  !WAIT  FOR  END  OF  DIGITIZING 

CALL  PUT(0,ID,NBYTES,MSTAT) 

CALL  WAIT(MSTAT)  !WAIT  FOR  END  OF  WRITTING  TO  TAPE 

C 

C  DISPLAY  DATA  ON  GT 

C 

CALL  DISP( NCHAN, NPTS) 

C 

FIN 

CALL  EOF<0) 

CALL  EOF(O) 

CALL  SKIP(0,-2) 

RETURN 

END 

^  ssaMSBSsaMsassaaasiaiaiaiaaaisssaaaasssssaissBSKSKBBSMaHSB 
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ccc 

CCC  DISP 
CCC 

SUBROUTINE  DISP( NCHAN.NPTS) 

COMMON  /DATA/  ID( 108) ,DATA(5 ,210) 

C  SUBROUTINE  TO  DISPLAY  WAVEFORM  TO  GT  DISPLAY 

C 

C  DISPLAY  DATA  ON  GT44 

C 

CALL  DEBASE(DATA, NCHAN.NPTS, 2, 1 ,NPTS,XBAR) 
CALL  GQDISP( 1,2, DATA ,NCHAN ,NPTS ,-2 ) 

CALL  DEBASE(DATA, NCHAN.NPTS, 3, 1 ,NPTS,XBAR) 
CALL  GQDISP( 2,3, DATA.NCHAN ,NPTS ,-2) 

CALL  DEBASE(DATA,NCHAN,NPTS,4,1 ,NPTS,XBAR) 
CALL  GQDISP( 3 , 4 .DATA, NCHAN , NPTS ,-2 ) 

CALL  DEBASE( DATA, NCHAN, NPTS.l ,1 ,NPTS,XBAR) 
CALL  GQDISP(4,1 .DATA.NCHAN , NPTS , -2) 

CALL  DEBASEC DATA, NCHAN, NPTS, 5, I , NPTS, XBAR) 
CALL  GQDISP(5,5,DATA,NCHAN,NPTS,-2) 

RETURN 

END 
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CCC 

CCC  ITWIXT 
CCC 

RETUEINS  AN  INTEGER  IN  THE  RANGE  GIVEN,  RANDOMLY  CHOSEN. 

INTEGER  FUNCTION  ITWIXT(MIN,MAX) 

IMPLICIT  INTEGER  (A-Z) 

REAL  PAN 

COMMON  /RAN/  S1,S2 
100  ITWIXT=RAN(S1 ,S2)*(MAX-MIN+1)+MIN 

IFdTWIXT.LT. MIN. OR. ITWIXT. GT. MAX)  GOTO  100 

RETURN 

END 
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ccc 

CCC  SUBROUTINE  KEYS 

CCC 

C 

C  CHECK  FOR  STOP  OR  QUIT  COMMANDS  (KEY  SUBROUTINE) 

C 

SUBROUTINE  KEYS(KEYFLG) 

C 

INTEGER  KEYFLG 
BYTE  KEY 
KEYFLG=0 
KEY=KEYCHK(IDUM) 

IF(KEY.NE.'P')  GOTO  10 
CALL  KEYOFF 

PAUSE  HIT  <CR>  TO  CONTINUE 
CALL  KEYON 

10  IF(KEY.NE.'Q')  GOTO  20 

CALL  KEYOFF 

WRITE(7,*)  'BLOCK  ABORTED;  BACK  TO  MENU  (2  EOFS  WRITTEN) 
KEYFLG=1 
20  RETURN 
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ccc 

CCC  TAPE 

CCC 

C 

C  SUBROUTINE  TO  MANIPULATE  TAPE. 

C 

SUBROUTINE  TAPE(NCHAN,NPTS,NIDS) 
IMPLICIT  INTEGER  (A-Z) 

BYTE  lOPTl 

COMMON  /DATA/  ID( 108) ,DATA(5,210) 


C 

C  FIND  NBYTES  FOR  TAPE  READ,  SKIP  A  LINE 
C 


NBYTES=(NIDS+NCHAN*NPTS)*2 
TYPE  ' 

C 

C  GET  THE  OPTION  FROM  THE  USER 
C 

5  imTE(7,10) 

10  FORMAT ('$TAPE:  (E)XIT,  (R)EWIND,', 

+  '  (L)ABEL,  (S)KIP,  (V)ERIFY,  (W)EOF>  ') 

READ(5,15)  lOPTl 
15  FORMAT(Al) 

C 

C  JUMP  ACCORDING  TO  OPTION 
C 


IF(IOPTl.EQ.'E') 
IF(IOPTl.EQ.'R') 
IF(IOPTl.EQ.'W') 
IFCIOPTI.EQ.'S') 
IF(IOPTl.EQ.'V') 
IF(IOPTl.EQ.'L') 
GOTO  5 


RETURN 

CALL  SKIP (0,0) 
CALL  EOF(O) 
GOTO  50 
GOTO  100 
GOTO  300 


C 

C  SKIP  EOF  OPTION 
C 

50  WRITE(7,55) 

55  FORMATC'S NUMBER  OF  EOFS  TO  SKIP  >') 

READ(5,*)  I FILES 
CALL  SKIP(0,IFILES) 

GOTO  5 


C 

C  VERIFY  TAPE  OPTION  -  READ  ONLY  IDS 
C 

100  DO  200  NT»1,9000 

CALL  GET(0, ID, NBYTES, MSTAT) 
CALL  WAIT( MSTAT) 
IF(MSTAT.NE.l)  GOTO  210 
IF(ID(1).NE. "177777)  GOTO  150 
WRITE(7,125)  (ID(L),L“2,50) 
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i23  FORMAT('  LABEL ,49A1 ) 

GOTO  200 

i'^0  lvTlITE(7,175)  NT,(ID(K),K=1 ,108) 

175  FORMAT('  RECORD:  ' , 16 , / , 1 1 ( 1016/ ) ) 

C  DISPLAY  DATA  ON  GT 

CALL  DISP(NCHAN,NPTS) 

C 

200  CONTINUE 

2i0  NT=NT-1 

GOTO  5 
C 

C  WRITE  A  LABEL 

C 

300  CALL  ZERO(ID,100) 

ID(1)="177777 

WRITE(7,311) 

311  FORMAT('$ LABEL  (LESS  THAN  50  CHARS)  >') 

READ(5,3i0)  (ID(L),L=2,50) 

310  F0RMAT(49A1) 

CALL  PUT(0,1D,100,MSTAT) 

CALL  WAIT(MSTAT) 

CALL  ZERO(ID,100) 

GOTO  5 
END 
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Appendix  I 


PMSM.FTX  Source  code  for  primary/secondary  memory  experiment  used  to 
collect  ERPS.  This  is  the  FORTRAN  source  code  (the  output  of  the  flees 
preprocessor)  and  must  be  compiled  using  the  FORTRAN  compiler* 
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PMSM.FTX 


C 
C 
C 

C  MODIFICATION  DATE:  JAN  16,  1988 
C 

C  PROGRAM  TO  CONTROL  PRIMARY/ SECONDARY  MEMORY  STUDY 
C 

C  IDS  OUTPUT  BY  THE  PROGRAM: 

C 


c 

ID(1) 

- 

SUBJECT  NUMBER 

c 

ID(2) 

- 

SESSION 

c 

ID(3) 

- 

BLOCK 

c 

ID(4) 

- 

DELAY 

c 

ID(5) 

- 

SOA  OF  RETRIEVAL  CUE 

c 

1D(6) 

- 

SIZE  OF  MEMORY  SET 

c 

ID(7) 

- 

CONSTANT  (1)  OR  VARIED  (2)  MAPPING 

c 

1D(8) 

- 

CODE  FOR  YES  BUTTON 

c 

ID(9) 

- 

CODE  FOR  NO  BUTTON 

C 

C  RUNNING  MEMORY  IDS 

C 

C  TRIAL  1 

C 

C  ID(IO)  -  RUNNING  MEMORY  STIMULUS 

C  ID(ll)  -  RT  (FOR  RUNNING  MEMORY  TRIAL) 

C  ID(12)  -  STIM  TYPE  FOR  RUNNING  MEMORY  (1-TARGET,  2-NONTARGET) 

C  ID(13)  -  RESPONSE  FOR  RUNNING  MEMORY  (O-NONE,  1-TARGET,  2-NONTARGET) 

C  ID(14)  -  CORRECTNESS  FOR  RUNNING  MEMORY:  O-NO  RESPONSE,  1-HIT, 

C  2-CORRECT  REJECTION,  3-FALSE  ALARM,  4-MISS 

C 

C  TRIAL  2 

C 

C  ID(15) 

C  ID(16) 

C  ID(17) 

C  ID(18) 

C  ID(19) 

C 
C 

C  TRIAL  3 

C 

C  ID(20) 

C  ID(21) 

C  ID(22) 

C  ID(23) 

C  ID(24) 

C 
C 

C  TRIAL  4 

C 

C  ID(25)  -  RUNNING  MEMORY  STIMULUS 


-  RUNNING  MEMORY  STIMULUS 

-  RT  (FOR  RUNNING  MEMORY  TRIAL) 

-  STIM  TYPE  FOR  RUNNING  MEMORY  (1-TARGET,  2-NONTARGET) 

-  RESPONSE  FOR  RUNNING  MEMORY  (O-NONE,  1-TARGET,  2-NONTARGET) 

-  CORRECTNESS  FOR  RUNNING  MEMORY:  O-NO  RESPONSE,  1-HIT, 
2-CORRECT  REJECTION,  3-FALSE  ALARM,  4-MISS 


-  RUNNING  MEMORY  STIMULUS 

-  RT  (FOR  RUNNING  MEMORY  TRIAL) 

-  STIM  TYPE  FOR  RUNNING  MEMORY  (1-TARGET,  2-NONTARGET) 

-  RESPONSE  FOR  RUNNING  MEMORY  (O-NONE,  1-TARGET,  2-NONTARGET) 

-  CORRECTNESS  FOR  RUNNING  MEMORY:  O-NO  RESPONSE,  1-HIT, 
2-CORRECT  REJECTION,  3-FALSE  ALARM,  4-MISS 
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C  ID(26)  -  RT  (FOR  RUNNING  MEMORY  TRIAL) 

C  ID(27)  -  STIM  TYPE  FOR  RUNNING  MEMORY  (1=TARGET,  2=NONTARGET) 

C  ID(28)  -  RESPONSE  FOR  RUNNING  MEMORY  (0=N0NE,  1-TARGET,  2-NONTARGET) 

C  ID(29)  -  CORRECTNESS  FOR  RUNNING  MEMORY:  O-NO  RESPONSE,  1-HIT, 

C  2-CORRECT  REJECTION,  3-FALSE  ALARM,  4-MISS 

C 

C  TRIAL  5 

C 

C  ID(30)  -  RUNNING  MEMORY  STIMULUS 

C  ID(31)  -  RT  (FOR  RUNNING  MEMORY  TRIAL) 

C  ID(32)  -  STIM  TYPE  FOR  RUNNING  MEMORY  (1-TARGET,  2-NONTARGET) 

C  ID(33)  -  RESPONSE  FOR  RUNNING  MEMORY  (O-NONE,  1-TARGET,  2-NONTARGET) 

C  ID(34)  -  CORRECTNESS  FOR  RUNNING  MEMORY:  O-NO  RESPONSE,  1-HIT, 

C  2-CORRECT  REJECTION,  3-FALSE  ALARM,  4-MISS 

C 

C  TRIAL  6 

C 

C  ID(35)  -  RUNNING  MEMORY  STIMULUS 

C  ID(36)  -  RT  (FOR  RUNNING  MEMORY  TRIAL) 

C  ID(37)  -  STIM  TYPE  FOR  RUNNING  MEMORY  (1-TARGET,  2-NONTARGET) 

C  ID(38)  -  RESPONSE  FOR  RUNNING  MEMORY  (O-NONE,  1-TARGET,  2-NONTARGET) 

C  ID(39)  -  CORRECTNESS  FOR  RUNNING  MEMORY:  O-NO  RESPONSE,  1-HIT, 

C  2-CORRECT  REJECTION,  3-FALSE  ALARM,  4-MISS 

C 

C  TRIAL  7 

C 

C  ID(40)  -  RUNNING  MEMORY  STIMULUS 

C  ID(41)  -  RT  (FOR  RUNNING  MEMORY  TRIAL) 

C  ID(42)  -  STIM  TYPE  FOR  RUNNING  MEMORY  (1-TARGET,  2-NONTARGET) 

C  ID(43)  -  RESPONSE  FOR  RUNNING  MEMORY  (O-NONE,  1-TARGET,  2-NONTARGET) 

C  ID(44)  -  CORRECTNESS  FOR  RUNNING  MEMORY:  O-NO  RESPONSE,  1-HIT, 

C  2-CORRECT  REJECTION,  3-FALSE  ALARM,  4-MISS 

C 

C  TRIAL  8 

C 

C  ID(45)  -  RUNNING  MEMORY  STIMULUS 

C  ID(46)  -  RT  (FOR  RUNNING  MEMORY  TRIAL) 

C  ID(47)  -  STIM  TYPE  FOR  RUNNING  MEMORY  (1-TARGET,  2-NONTARGET) 

C  ID(48)  -  RESPONSE  FOR  RUNNING  MEMORY  (O-NONE,  1-TARGET,  2-NONTARGET) 

C  ID(49)  -  CORRECTNESS  FOR  RUNNING  MEMORY:  O-NO  RESPONSE,  1-HIT, 

C  2-CORRECT  REJECTION,  3-FALSE  ALARM,  4-MISS 

C 

C  TRIAL  9 

C 

C  ID(50)  -  RUNNING  MEMORY  STIMULUS 

C  ID(51)  -  RT  (FOR  RUNNING  MEMORY  TRIAL) 

C  ID(52)  -  STIM  TYPE  FOR  RUNNING  MEMORY  (1-TAPGET,  2-NONTARGET) 

ID(53)  -  RESPONSE  FOR  RUNNING  MEMORY  (O-NONE,  1-TARGET,  2-NONTARGET) 
ID(54)  -  CORRECTNESS  FOR  RUNNING  MEMORY:  O-NO  RESPONSE,  1-HIT, 
2-CORRECT  REJECTION,  3-FALSE  ALARM,  4-MISS 
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C  TRIAL  10 
C 

C  ID(55)  -  RUNNING  MEMORY  STIMULUS 

C  ID(56)  -  RT  (FOR  RUNNING  MEMORY  TRIAL) 

C  ID(57)  -  STIM  TYPE  FOR  RUNNING  MEMORY  (1-TARGET,  2-NONTARGET) 

C  ID(58)  -  RESPONSE  FOR  RUNNING  MEMORY  (O-NONE,  1-TARGET,  2-NONTARGET) 

C  ID(59)  -  CORRECTNESS  FOR  RUNNING  MEMORY:  O-NO  RESPONSE,  1-HIT, 

C  2-CORRECT  REJECTION,  3-FALSE  ALARM,  4-MISS 

C 

C  TRIAL  11 

C 

C  ID(60)  -  RUNNING  MEMORY  STIMULUS 

C  ID(61)  -  RT  (FOR  RUNNING  MEMORY  TRIAL) 

C  ID(62)  -  STIM  TYPE  FOR  RUNNING  MEMORY  (1-TARGET,  2-NONTARGET) 

C  ID(63)  -  RESPONSE  FOR  RUNNING  MEMORY  (O-NONE,  1-TARGET,  2-NONTARGET) 

C  ID(64)  -  CORRECTNESS  FOR  RUNNING  MEMORY:  O-NO  RESPONSE,  1-HIT, 

C  2-CORRECT  REJECTION,  3-FALSE  ALARM,  4-MISS 

C 

C  TRIAL  12 

C 

C  ID(65) 

C  ID(66) 

C  ID(67) 

C  ID(68) 

C  ID(69) 

C 
C 

TRIAL  13 

C  ID(70) 

C  ID(71) 

C  ID(72) 

C  ID(73) 

C  ID(74) 

C 
C 

C  TRIAL  14 

C 

C  ID(75) 

C  ID(76) 

C  ID(77) 

C  ID(78) 

C  ID(79) 

C 
C 

C  TRIAL  15 

C 

C  ID(80)  -  RUNNING  MEMORY  STIMULUS 

C  ID(81)  -  RT  (FOR  RUNNING  MEMORY  TRIAL) 

C  ID(82)  -  STIM  TYPE  FOR  RUNNING  MEMORY  (1-TARGET,  2-NONTARGET) 

C  ID(83)  -  RESPONSE  FOR  RUNNING  MEMORY  (O-NONE,  1-TARGET,  2-NONTARGET) 


-  RUNNING  MEMORY  STIMULUS 

-  RT  (FOR  RUNNING  MEMORY  TRIAL) 

-  STIM  TYPE  FOR  RUNNING  MEMORY  (1-TARGET,  2-NONTARGET) 

-  RESPONSE  FOR  RUNNING  MEMORY  (O-NONE,  1-TARGET,  2-NONTARGET) 

-  CORRECTNESS  FOR  RUNNING  MEMORY:  O-NO  RESPONSE,  1-HIT, 
2-CORRECT  REJECTION,  3-FALSE  ALARM,  4-MISS 


-  RUNNING  MEMORY  STIMULUS 

-  RT  (FOR  RUNNING  MEMORY  TRIAL) 

-  STIM  TYPE  FOR  RUNNING  MEMORY  (1-TARGET,  2-NONTARGET) 

-  RESPONSE  FOR  RUNNING  MEMORY  (O-NONE,  1-TARGET,  2-NONTARGET) 

-  CORRECTNESS  FOR  RUNNING  MEMORY:  O-NO  RESPONSE,  1-HIT, 
2-CORRECT  REJECTION,  3-FALSE  ALARM,  4-MISS 


-  RUNNING  MEMORY  STIMULUS 

-  RT  (FOR  RUNNING  MEMORY  TRIAL) 

-  STIM  TYPE  FOR  RUNNING  MEMORY  (1-TARGET,  2-NONTARGET) 

-  RESPONSE  FOR  RUNNING  MEMORY  (O-NONE,  1-TARGET,  2-NONTARGET) 

-  CORRECTNESS  FOR  RUNNING  MEMORY:  O-NO  RESPONSE,  1-HIT, 
2-CORRECT  REJECTION,  3-FALSE  ALARM,  4-MISS 
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ID(84)  -  CORRECTNESS  FOR  RUNNING  MEMORY:  0=NO  RESPONSE,  1=HIT, 
2=CORRECT  REJECTION,  3=FALSE  ALARM,  4=MISS 

STERNBERG  IDS 
C 

C  ID(85)  -  MEMORY  SET  WORD  //I,  LETTER  #1 

C  ID(86)  -  MEMORY  SET  WORD  //I,  LETTER  //2 

C  ID(87)  -  MEMORY  SET  WORD  //I,  LETTER  #3 

C  ID(88)  -  MEMORY  SET  WORD  //I,  LETTER  //4 

C 

C  ID(89)  -  MEMORY  SET  WORD  #2,  LETTER  //I 

C  ID(90)  -  MEMORY  SET  WORD  #2,  LETTER  #2 

C  ID(91)  -  MEMORY  SET  WORD  //2,  LETTER  //3 

C  ID(92)  -  MEMORY  SET  WORD  #2,  LETTER  //4 

C 

C  ID(93)  -  MEMORY  SET  WORD  //3,  LETTER  #1 

C  ID(94)  -  MEMORY  SET  WORD  #3,  LETTER  #2 

C  ID(95)  -  MEMORY  SET  WORD  #3,  LETTER  if 3 

C  ID(96)  -  MEMORY  SET  WORD  #3,  LETTER  #4 

C 

C  ID(97)  -  MEMORY  SET  WORD  ifk ,  LETTER  if  I 

C  ID(98)  -  MEMORY  SET  WORD  ifh ,  LETTER  if  2 

C  ID(99)  -  MEMORY  SET  WORD  /M ,  LETTER  if  3 

C  ID(IOO)  -  MEMORY  SET  WORD  #4,  LETTER  if  A 
C 

C  ID(lOl)  -  RT  (FOR  STERNBERG  TRIAL) 

C  ID(102)  -  STIM  TYPE  (1=TARGET,  2=NONTARGET) 

C  ID(103)  -  RESPONSE  (0=NONE,  1*TARGET,  2=NONTARGET) 

C  ID(104)  -  CORRECTNESS:  0=NO  RESPONSE,  l=HIT,  2»CORRECT  REJECTION, 

C  3=FALSE  ALARM,  4-MISS 

C 

C  ID(105)  -  STIMULUS  PRESENTED,  LETTER  if  I 

C  ID(106)  -  STIMULUS  PRESENTED,  LETTER  if2 

C  ID(107)  -  STIMULUS  PRESENTED,  LETTER  if 3 

C  ID(108)  -  STIMULUS  PRESENTED,  LETTER  //4 

C 

C - 

C 

C 

IMPLICIT  INTEGER  (A-Z) 

INTEGER  COND( 96,7), LSTIM( 3 ) , ISEED( 2 ) 

BYTE  NAME( 17), ERR, TEMPI ,WORD(5) ,NUM(9) ,DIGIT(2) ,YES,NO, 
+CMTARG(5,4),CMDIST(5,4),VMSTIM(10,4),TARG(5,4),DIST(5,4) 

REAL  SEED, RAN 

EQUIVALENCE  ( SEED , ISEED) , ( ISEED( 1 ) , S 1 ) , ( ISEED( 2 ) , S2 ) 

COMMON  /RAN/  S1,S2 

COMMON  /DATA/  ID( 108) ,DATA(5,210) 

C 

DATA  STMDUR,NBACK/200,2/ 

DATA  NCHAN ,NPTS , DRATE , NIDS ,BASE/ 5 , 2 10 , 10 , 108 , 100/ 

DATA  NUM/'1','2','3','4','5','6','7','8','9'/ 
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c 

C  GET  A  RANDOM  SEED  AND  WARM  UP  RANDOM  NUMBER  GENERATOR 

(' 

Kj 

SEED=SECNDS(0) 

DO  32765  RANDOM=l ,2000 
CALL  RAN(S1 ,S2) 

32765  CONTINUE 
C 

C  INITIALIZE  MATROX 

C 

32766  CALL  MCINIT(l) 

C  SET  UP  DISPLAY  ASSIGNMENTS 

C 

CALL  GQINIT 

CALL  GQTTL('  PM/3M  EXPERIMENT') 

CALL  GQLABL(1,'FZ  ELECTRODE') 

CALL  GQLABL(2,'CZ  ELECTRODE') 

CALL  GQLABL(3,'PZ  ELECTRODE') 

CALL  GQLABL(4,'EOG  ELECTRODE') 

CALL  GQLABL(5,'STM  MARKER') 

C 

C  ASSIGN  STIMULUS  FILE  FROM  DISC 

C 

TYPE  INPUT  STIMULUS  FILE  FOR  THIS  SUBJECT' 

CALL  GETSTR(5,NAME,16,ERR) 

CALL  ASSIGN(40,NAME,0,'OLD') 

C 

C  READ  STIMULUS  FILE 

C 

DO  32763  1=1,5 

READ(40,1)  (CMTARG(I,J),J=1,4) 

32763  CONTINUE 
C 

32764  DO  32761  1=1,5 

READ(40,1)  (CMDIST(I,J),J=1,4) 

32761  CONTINUE 
C 

32762  DO  32759  1=1,10 

READ(40,1)  (VMSTIM(I,J),J=1,4) 

32759  CONTINUE 
C 

32760  READ(40,2)  YES, NO 

IF( .NOT.(YES.EQ.'/'))  GO  TO  32758 

YESBUT=1 

NOBUT=0 

32758  IF(.NOT.(YES.EQ.'Z'))  GO  TO  32757 
YESBUT=0 
NOBUT=l 
C 

32757  DO  32755  J=l,96 
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READ(40,3)  (C0ND(J,I),I-1,7) 

3275:)  CONTINUE 
C 

32756  CONTINUE 

1  FORMAT(AAl) 

2  FORMAT ( 2 Al) 

3  FORMAT(7I5) 

C 

C  «-==«-»=.===.==  get  the  option  FROM  THE  USER  ——«======== 

C 

C  GET  THE  OPTION 

C 

100  CONTINUE 

TYPE  110 

no  FORMAT(//'$  MENU:  (R)UN,  (T)APE,  (C)ALIB,  (Q)UIT->') 
ACCEPT  120,IOPT 
120  FORMAT(Al) 

C 

C  JUMP  ON  OPTION 

C 


IF(IOPT.EQ.'C') 
IF(IOPT.EQ.'Q') 
IF(IOPT.EQ.'R') 
IF(IOPT.EQ.'T') 
GOTO  100 


CALL  CALIB(NCHAN,NPTS,DRATE,NIDS) 
GOTO  2000 
GOTO  500 

CALL  TAPE(NCHAN,NPTS,NIDS) 


C 

C  RUN  TIME  SETUP 

C 

C  TRAP  "Q"  FOR  ABORT. 

C 

500  ISWFLG-0 

CALL  SCCA(ISWFLG) 

CALL  KEYON  ITURN  OFF  ECHO,  CHECK  FOR  KEY  PRESSES 

C 

C  INITIALIZE  NBYTES  FOR  TAPE  WRITE 

C 

NBYTES-(NIDS+NCHAN*NPTS)*2 

C 

C  INITIALIZE  DEVICES 

C 

CALL  LINOFF 

CALL  CLINIT 

CALL  DIINIT 

CALL  DOINIT 

DO  32753  ICLOCK-2,5 

CLOCK- I CLOCK 

CALL  PCINIT( CLOCK) 

32753  CONTINUE 
C 

C— — — — —  MAIN  RUN  TIME  LOOP— —————— 

C 
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32754  DO  32751  BLOCK-1,96 
C 

C  ZERO  ID  ARRAY 

C 

CALL  ZERO(ID,NIDS*2) 

C 

C  CHECK  FOR  "Q"  AND  "P" 

C 

CALL  KEYS(KEYFLG)  ! CHECK  FOR  KEY  PRESSES 

IF(.NOT.(KEYFLG.EQ.l))  GO  TO  32750 

CALL  CLKILL 

CALL  MCERAS 

CALL  DIKILL 

CALL  DOKILL 

CALL  ADKILL 

CALL  KEYOFF 

DO  32748  ICLOCK-2,5 

CLOCK- I CLOCK 

CALL  PCKILLC CLOCK) 

32748  CONTINUE 

32749  CALL  EOF(O) 

CALL  EOF(O) 

CALL  SKIP (0,-2) 

GOTO  100 

C 

32750  IBLK-BLOCK 

iF(.NOT.(((IBLK/16)*16).EQ.IBLK.AND.lBLK.NE.96))  GO  TO  32747 
TYPE  * Program  paused'’ 

CALL  MCTEXT( 20, 200, 'Pause  —  Press  button  to  continue  ') 

CALL  MCDISP(l) 

C 

C  WAIT  FOR  SUBJECT  TO  PRESS  BUTTON 
C 

CALL  DISET(0, 'SINGLE', ISTATO) 

CALL  DISETd, 'SINGLE', ISTATl) 

CALL  WAIT( ISTATO , ISTATl ) 

CALL  MCDISP(O) 

CALL  MCERAS 
C 

C  READ  CONDITION,  STORE  IDS  AND  SELECT  STIMULI 
C 

32747  SUB  -COND ( BLOCK, 1 ) 

SESS  -COND(BLOCK,2) 

SOA  -COND(BLOCK,3) 

DELAY  -COND(BLOCK,4) 

CMVM  -COND(BLOCK,5) 

MEMSI Z-COND( BLOCK , 6 ) 

STYPE  -COND (BLOCK, 7) 

C 

ID(1)-  SUB 
ID(2)-  SESS 
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1D(3)-  BLOCK 
ID(4)*  DELAY 
ID(5)-  SOA 
ID(6)-  MEMSIZ 
ID(7)»  CMVM 
ID(8)“  YESBUT 
ID(9)-  NOBUT 
C 

C  INFORM  EXPERIMENTER  OF  CONDITION 

C 

TYPE  30 , BLOCK , DELAY , CMVM , MEMSIZ , STYPE 
30  FORMAT($'  Block  ',13,'  Delay  ',12,'  CMVM  ',11, 
+'  PosNeg  ',11) 

C 

DO  32745  K-1,10 
C 

C  PERMUTE  CM  TARGETS 

C 

DO  32743  LOOP  «1,5 
J-ITWIXT(1,5) 

DO  32741  1-1,4 
TEMPI  -  CMTARG(J,I) 

CMTARG(J,I)  -  CMTARG(L00P,1) 

CMTARG(LOOP,I)  -  TEMPI 

32741  CONTINUE 

32742  CONTINUE 

32743  CONTINUE 
C 

C  PERMUTE  CM  DISTRACTORS 
C 

32744  DO  32739  LOOP  -1,5 
J-ITWIXT(1,5) 

DO  32737  1-1,4 
TEMPI  -  CMDIST(J,I) 

CMDIST(J,I)  -  CMDIST(LOOP,I) 

CMDIST(L00P,I)  -  TEMPI 

32737  CONTINUE 

32738  CONTINUE 

32739  CONTINUE 
C 

C  PERMUTE  VM  STIMULI 
C 

32740  DO  32735  LOOP  -1,10 
J-ITWIXT(1,10) 

DO  32733  1-1,4 
TEMPI  -  VMSTIM(J,1) 

VMSTIM(J,I)  -  VMSTIM(LOOP,I) 

VMSTIM(LOOP,I)  -  TEMPI 

32733  CONTINUE 

32734  CONTINUE 

32735  CONTINUE 


Memsiz  ', 
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32736  CONTINUE 

32745  CONTINUE 
C 

C  CHOOSE  CM  STIMULI 
C 

32746  IF(.NOT.(CMVM.EQ.l))  GO  TO  32732 
DO  32730  J=1,MEMSIZ 

DO  32728  1=1,4 
TARG(J,I)=CMTARG(J,I) 

32728  CONTINUE 

32729  CONTINUE 

32730  CONTINUE 
C 

32731  DO  32726  J=l,5 
DO  32724  1=1,4 

DIST( J,I)=CMDIST(J,I) 

32724  CONTINUE 

32725  CONTINUE 

32726  CONTINUE 
C 

32727  CONTINUE 
C 

C  CHOOSE  VM  STIMULI 
C 

32732  IF(.NOT.(CMVM.EQ.2))  GO  TO  32723 
DO  32721  J=1,MEMSIZ 

DO  32719  1=1,4 
TARG( J , I ) =VMSTIM( J , I ) 

32719  CONTINUE 

32720  CONTINUE 

32721  CONTINUE 
C 

32722  DO  32717  J=MEMSIZ+1 ,MEMSIZ+5 
DO  32715  1=1,4 

DIST( J-MEMSIZ , I)=VMSTIM( J , I ) 

32715  CONTINUE 

32716  CONTINUE 

32717  CONTINUE 
C 

32718  CONTINUE 
C 

C  ==================1  PRESENT  MEMORY  SET  ===■ 

C 

C 

C  WRITE  STIMULI  TO  MATROX  AND  STORE  IDS 
C 

32723  DO  32713  I=1,MEMSIZ 
C 

IF(.NOT.(MEMSIZ.EQ.2))  GO  TO  32712 
XP0S=60+(I*35) 

YPOS-lOO 


no 


c 

32712  IF(.N0T.(MEMSIZ.EQ.4))  GO  TO  32711 
XPOS-27+(I*35) 

YPOS-lOO 

C 

32711  DO  32709  J-1,4 

WORD(J)-TARG(I,J) 

32709  CONTINUE 
C 

32710  DO  32707  J-1,4 
ID(80+(4*I)+J)-TARG(1,J) 

32707  CONTINUE 

32708  WORD(5)"0 

CALL  MCTEXT(XPOS,YPOS,WORD) 

32713  CONTINUE 
C 

C  TURN  ON  STIMULUS 
C 

32714  CALL  MCSYNC  !  SUNCH  DISPLAY 

CALL  MCDISP(l)  !  TURN  ON  STIMULUS 

CALL  PCSET( 1+1 SINGLE' ,3 ,3000,IFLAG1 ) 

CALL  WAIT(IFLAGl)  !  WAIT  FOR  3  SECONDS 
CALL  MCDISP(O)  !  TURN  OFF  STIMULUS 

CALL  MCERAS  !  ERASE  MATROX 

C 

C  WAIT  FOR  1500  MSEC  BEFORE  PRESENTING  PROBES 
C 

CALL  PCSET( 1+1 SINGLE' , 3 , 1500 .IFLAGl ) 

CALL  WAIT( IFLAGl) 

C 

C  RUNNING  MEMORY  TRIALS  — — 

C 

IF(. NOT. (DELAY. NE.O))  GO  TO  32706 
DO  32704  ITRL-1,IABS(DELAY) 

C 

C  PREPARE  STIMULI  DURING  PRE-STIMULUS  INTERVAL 
C 

CALL  PCSET( 1+1 ,' SINGLE', 3, 5, IFLAGl) 
P0SNEG-ITWIXT(1,3) 

IF(P0SNEG.EQ.3)  POSNEG-2 
C 

C  IF  A  POSITIVE  TRIAL,  REPEAT  STIMULUS  NBACK 

C  IF  A  NEGATIVE  TRIAL,  SELECT  A  NEW  STIMULUS 

C 

IF(.NOT.(POSNEG.EQ.l))  GO  TO  32703 
IF( ITRL . LE . NBACK)  JSTIM-ITWIXT( 1,9) 

IF( ITRL. GT. NBACK)  JSTIM-LSTIM(NBACK) 

C 

32703  IF(.NOT.(POSNEG.EQ.2))  GO  TO  32702 
10  JSTIM-ITWIXT(1,9) 

IF(JSTTM.EQ.LSTIM( NBACK))  GOTO  10 


c 

■;  STORE  SUCCESSIVE  STIMULI 

12702  LSTIM(3)  =  LSTIM(2) 

LSTIM(2)  =  LSTIM(l) 

LSTIM(l)  =  JSTIM 

ID((ITRL*5)+5)  =  JSTIM 

XP0S=120 

YP0S=120 

DIGIT( 1)=NUM( JSTIM) 

DIGIT(2)=0 

IF  (DELAY. GT.O)  CALL  MCTEXT( XPOS ,YPOS .DIGIT) 

(J 

CALL  WAIT(IFLAGl)  !  WAIT  OUT  PRESTIMULUS  INTERVAL 

r 

CALL  MCSYNC  !  SYNCH  DISPLAY 

CALL  MCDISP(l)  !  TURN  ON  STIMULUS 

CALL  PCSET( 1+1 , 'SINGLE' , 3 , 1000 .IFLAGl ) 

CALL  DIRESP(RESP, 'SINGLE', 2, 1+1 .RT.RSTAT) 

CALL  PCSET( 1+2 , 'SINGLE' ,3,200, 1FLAG2 ) 

CALL  WAIT(IFLAG2) 

CALL  MCDISP(O)  !  TURN  OFF  STIMULUS 

CALL  MCERAS  !  ERASE  MATROX 

CALL  WA1T( IFLAGl)  IWAIT  FOR  END  OF  1000  MSEC  INTERVAL 

C 

C  CLASSIFY  TRIAL  AND  STORE  RT  AND  STIM  TYPE  (1-HIT,  2-CR,  3-FA,  4-MlSS) 

C 

ID((ITRL*5)+6)  =  RT 
ID((ITRL*5)+7)  =  POSNEG 
ID((ITRL*5)+8)  =  0 
ID((ITRL*5)+9)  =  0 

i'j 

IF(.NOT.(RSTAT.NE.O))  GO  TO  32701 

C 

IF(.NOT.(RESP.EQ.YESBUT))  GO  TO  32700 
ID((ITRL*5)+8)  =  1 

IF  (ID((ITRL*5)+7).EQ.l)  ID(( ITRL*5)+9)  =  1 
IF  (ID((ITRL*5)+7).EQ.2)  ID( (ITRL*5)+9)  =  3 

C 

J2700  IF(.NOT.(RESP.EQ.NOBUT))  GO  TO  32699 
ID((ITRL*5)+8)  =  2 

IF  (IU((ITRL*5)+7) .EQ.2)  ID(( ITRL*5)+9)  =  2 
IF  (ID((ITRL*5)+7).EQ.l)  ID( ( ITRL*5)+9)  -  4 

C 

32699  CONTINUE 

C 

32701  CONTINUE 
12704  CONTINUE 
■>2  703  CONTINUE 
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c  STERNBERG  PROBE  TRIAL 

C 

C 

C  START  THE  DIGITIZER 

C 

32706  CALL  ADSET(DATA,2,1 ,NCHAN,NPTS,3,DRATE,DSTAT) 

C 

C  WAIT  OUT  BASELINE  CLOCK 

C 

CALL  PCSET( 1+1 , 'SINGLE' ,3 .BASE, IFLAGl) 

CALL  WAIT(IFLAGl) 

C 

C  DELIVER  RETRIEVAL  CUE 

C 

IF(.NOT.(SOA.GE.200))  GO  TO  32698 

XPOS-120 

YP0S=120 

CALL  MCCHAR(XPOS,YPOS,'*') 

CALL  MCSYNC  !  SYNCH  DISPLAY 

CALL  DOSET(O)  !TURN  ON  STIMULUS  MARKER 

CALL  MCDISP(l)  !  TURN  ON  MATROX 

CALL  PCSET( 1+1 , 'SINGLE' ,3 .SOA.IFLAGl) 

CALL  PCSET( 1+2, 'SINGLE', 3, 200, IFLAG2) 

CALL  WAIT(IFLAG2) 

CALL  MCDISP(O)  !  TURN  OFF  MATROX 

CALL  DOCLR(O)  !  TURN  OFF  STIMULUS  MARKER 

CALL  MCERAS  !  ERASE  MATROX 

CALL  WAIT(IFLAGl)  !  WAIT  FOR  END  OF  SOA  INTERVAL 
C 

C  PREPARE  STIMULI  DURING  PRE-STIMULUS  INTERVAL 
C 

32698  CALL  PCSET( 1+1 ,'SINGLE' , 3 ,5 .IFLAGl ) 

XPOS-112 

YPOS-140 

C 

IF(.NOT.(STYPE.EQ.l))  GO  TO  32697 
I-ITWIXT(1 .MEMSIZ) 

DO  32695  J-1,A 
WORD(J)-TARG(I,J) 

ID(104+J)-TARG(I,J) 

32695  CONTINUE 

32696  CONTINUE 
C 

32697  IF(.N0T.(STYPE.EQ.2))  GO  TO  32694 
I-ITWIXT(1,5) 

DO  32692  J-1,4 
WORD(J)-DIST(I,J) 

ID(104+J)-DIST(I,J) 

32692  CONTINUE 

32693  CONTINUE 

32694  WORD(5)-0 
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CALL  MCTEXT(XPOS,YPOS,WORD) 

CALL  WAIT(IFLAGl) 

C 

CALL  MCSYNC  !  SYNCH  DISPLAY 

CALL  DOSET(O)  !  TURN  ON  STIMULUS  MARKER 

CALL  MCDISP(l)  !  TURN  ON  DISPLAY 

CALL  PCSET( 1+1 , 'SINGLE', 3, 3000, IFLAGl) 

CALL  DIRESP(RESP, 'SINGLE', 2, 1+1, RT.RSTAT) 

CALL  PCSET( 1+2 , ' S INGLE' ,3,200 , IFLAG2 ) 

CALL  WAIT(IFLAG2) 

CALL  MCDISP(O)  !  TURN  OFF  DISPLAY 

CALL  DOCLR(O)  !  TURN  OFF  STIMULUS  MARKER 

CALL  MCERAS  !  CLEAR  MARTOX 

CALL  WAIT( IFLAGl)  !  WAIT  FOR  END  OF  3  SECOND  INTERVAL 
C 

C  CLASSIFY  TRIAL  AND  STORE  RT  AND  STIM  TYPE  (1=HIT,  2-CR,  3=FA,  4=MISS) 
C 

ID(lOl)  =  RT 
ID(102)  =  STYPE 
ID(103)  =  0 
ID(104)  =  0 
C 

IF(.NOT.(RSTAT.NE.O))  GO  TO  32691 
IF(.NOT.(RESP.EQ.YESBUT))  GO  TO  32690 
ID(103)  *  1 

IF  (ID(102).EQ.l)  ID(104)  =  1 

IF  (ID(102).EQ.2)  ID(104)  =  3 

C 

32690  IF(.NOT.(RESP.EQ.NOBUT))  GO  TO  32689 
ID(103)  =  2 

IF  (ID(102).Eq.2)  ID(104)  =  2 

IF  (ID(102).EQ.l)  ID(104)  =  4 

C 

32689  CONTINUE 
C 

C  CHECK  IF  ANY  POINTS  ARE  OUT  OF  BOUNDS 
C 

32691  ERROR^O 

DO  32687  IPTS=1,210 
DO  32685  lCHAN-1,4 

I F( I ABS ( DATA( ICHAN , IPTS ) ) . GE . 204 7 )  ERROR-ERROR+1 

32685  CONTINUE 

32686  CONTINUE 

32687  CONTINUE 
C 

C  WRITE  DATA  TO  TAPE 

C 

32688  CALL  PUT(0,ID,NBYTES,MSTAT) 

TYPE  20,  ID(101),ID(104) 

20  FORMAT('  RT  ',14,'  Corr  ',11) 

IF  (ERROR.NE.O)  TYPE  *,'Eiror  in  last  trial  of  digitized  EEG' 


CALL  DISP(NCHAN,NPTS) 

CALL  WAIT(MSTAT) 

C 

32751  CONTINUE 
C 

C  KILL  THE  DEVICES 
C 

32752  CALL  CLKILL 
CALL  MCERAS 
CALL  DIKILL 
CALL  DOKILL 
CALL  ADKILL 
CALL  KEYOFF 

DO  32683  ICLOCK-2,5 
CLOCK- I CLOCK 
CALL  PCKILLC CLOCK) 

32683  CONTINUE 
C 

C  WRITE  DOUBLE  EOF  AND  SKIP  BACK  1 
C 

32684  CALL  EOF(O) 

CALL  EOF(O) 

CALL  SKIP(0,-2) 

C 

C  RETURN  TO  MAIN  MENU 
C 

GOTO  100 
C 

C  QUIT  SECTION 

C 

C  TURN  OFF  THE  DISPLAY  AND  EXIT  PROGRAM 
C 

2000  CONTINUE 

CALL  CLKILL 
CALL  MCERAS 
CALL  GQINIT 
CALL  DIKILL 
CALL  ADKILL 
CALL  KEYOFF 
STOP 
END 
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c  * 

C  SUBROUTINES  * 

C  * 

ccc 

CCC  SUBROUITNE  CALIB 

CCC 

C 

C  GATHERS  FIVE  CALIBRATION  TRIALS 
C 

SUBROUTINE  CALI B( NCHAN , NPTS , DRATE , NIDS ) 

IMPLICIT  INTEGER  (A-Z) 

C 

COMMON  /DATA/  ID( 108) ,DATA(5 ,210) 

NBYTES=(NIDS+NCHAN*NPTS)*2 

C 

CALL  LINOFF 
CALL  CLINIT 
WR1TE(7,10) 

10  FORMAT( /////'  TURN  ON  CALIBRATION  PULSE') 

PAUSE  'HIT  <CR>  TO  CONTINUE' 

C  DIGITIZE  5  TRIALS  OF  CALIBRATION  DATA 

DO  32765  L00P=1,5 

CALL  ADSETC DATA, 2,1, NCHAN, NPTS, 3, DRATE, DSTAT) 

CALL  WAIT( DSTAT)  '.WAIT  FOR  END  OF  DIGITIZING 

CALL  PUT(0,ID,NBYTES,MSTAT) 

CALL  WAIT(MSTAT)  !WAIT  FOR  END  OF  WRITTING  TO  TAPE 

C 

C  DISPLAY  DATA  ON  GT 

C 

CALL  DISP( NCHAN, NPTS) 

C 

32765  CONTINUE 

32766  CALL  EOF(O) 

CALL  EOF(O) 

CALL  SKIP(0,-2) 

RETURN 

END 

Q  sassassasasssassssssssBsassasssasassssBsssssrBSBsasEaas 
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ccc 

CCC  DISP 
CCC 

SUBROUTINE  DISP(NCHAN,NPTS) 

COMMON  /DATA/  ID( 108) ,DATA(5 ,210) 

C  SUBROUTINE  TO  DISPLAY  WAVEFORM  TO  GT  DISPLAY 
C 

C  DISPLAY  DATA  ON  GT44 
C 

CALL  DEBASE ( DATA , NCHAN , NPTS , 2 , 1 , NPTS , XBAR) 
CALL  GQDI SP ( 1 , 2 , DATA , NCHAN , NPTS , -2 ) 

CALL  DE  BASE ( DATA , NCHAN , NPTS , 3 , 1 , NPTS , XBAR) 
CALL  GQDISP(2,3,DATA,NCHAN,NPTS,-2) 

CALL  DEBASE ( DATA , NCHAN , NPTS , 4 , 1 , NPTS , XBAR) 
CALL  GQDISP( 3 ,4 .DATA, NCHAN ,NPTS ,-2) 

CALL  DEBASE( DAT A, NCHAN, NPTS ,1,1 ,NPTS ,XBAR) 
CALL  GQDISP(4,1 ,DATA,NCHAN,NPTS,-2) 

CALL  DEBASE( DATA, NCHAN , NPTS ,5,1 ,NPTS , XBAR) 
CALL  GQDISP' 5 , 5 .DATA , NCHAN , NPTS ,-2 ) 

RETURN 

END 
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ccc 

CCC  ITWIXT 

CCC 

c 

C  RETURNS  AN  INTEGER  IN  THE  RANGE  GIVEN,  RANDOMLY  CHOSEN. 

C 

INTEGER  FUNCTION  ITWIXT( MIN ,MAX) 

IMPLICIT  INTEGER  (A-Z) 

REAL  RAN 

COMMON  /RAN/  SI ,S2 

100  ITWIXT=RAN(S1 ,S2)*(MAX-MIN+1)+MIN 

IF( ITWIXT. LT. MIN. OR. ITWIXT. GT. MAX)  GOTO  100 

RETURN 

END 
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CCC 

CCC  SUBROUTINE  KEYS 
CCC 

CHECK  FOR  STOP  OR  QUIT  COMMANDS  (KEY  SUBROUTINE) 

SUBROUTINE  KEYS(KEYFLG) 

C 

INTEGER  KEYFLG 
BYTE  KEY 
KEYFLG-O 
KEY-KEYCHK(IDUM) 

IF(KEY.NE.'P')  GOTO  10 
CALL  KEYOFF 

PAUSE  HIT  <CR>  TO  CONTINUE  — ' 

CALL  KEYON 

10  IF(KEY.NE.'Q')  GOTO  20 

CALL  KEYOFF 

WRITE(7,*)  'BLOCK  ABORTED:  BACK  TO  MENU  (2  EOFS  WRITTEN)' 
KEYFLG- 1 
20  RETURN 

END 


119 


ccc 

CCC  TAPE 

CCC 

C 

C  SUBROUTINE  TO  MANIPULATE  TAPE. 

C 

SUBROUTINE  TAPE(NCHAN,NPTS,NIDS) 

IMPLICIT  INTEGER  (A-Z) 

BYTE  lOPTl 

COMMON  /DATA/  ID( 108) ,DATA( 5 ,210) 

C 

C  FIND  NBYTES  FOR  TAPE  READ,  SKIP  A  LINE 
C 

NBYTES=(NIDS+NCHAN*NPTS)*2 
TYPE  ' 

C 

C  GET  THE  OPTION  FROM  THE  USER 
C 

5  WRITE(7,10) 

10  FORMAT("$TAPE:  (E)XIT,  (R)EVn:ND,', 

+  '  (L)ABEL,  (S)KIP,  (V)ERIFY,  (W)EOF>  ') 

READ(5,15)  lOPTl 
15  FORMAT(Al) 

C 

C  JUMP  ACCORDING  TO  OPTION 
C 

IFCIOPTl.EQ.'E")  RETURN 
IF(IOPTl.EQ.'R')  CALL  SKIP(0,0) 
IF(IOPTl.EQ.'W')  CALL  EOF(O) 

IF(IOPTl.Eq.'S')  GOTO  50 
IF(IOPTl.EQ.'V')  GOTO  100 
IF(IOPTl.EQ.'L')  GOTO  300 
GOTO  5 
C 

C  SKIP  EOF  OPTION 
C 

50  WRITE(7,55) 

55  FORMAT('$ NUMBER  OF  EOFS  TO  SKIP  >') 

READ(5,*)  I FILES 
CALL  SKIP(0,IFILES) 

GOTO  5 
C 

C  VERIFY  TAPE  OPTION  -  READ  ONLY  IDS 
C 

100  DO  200  NT=1,9000 

CALL  GET(0, ID, NBYTES, MSTAT) 

CALL  WAIT( MSTAT) 

IF(MSTAT.NE.l)  GOTO  210 
IF(ID(1).NE. ”177777)  GOTO  150 
WRITE(7,125)  (ID(L),L-2,50) 

125  FORMAT('  LABEL; ' ,49A1 ) 
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GOTO  200 

150  WR1TE(7,175)  NT,(ID(K) ,K-1 ,108) 

175  FORMAK'  RECORD:  ',I6,/,11(  1016/)) 

C 

C  DISPLAY  DATA  ON  GT 
C 

CALL  DISP(NCHAN,NPTS) 

C 

200  CONTINUE 

210  NT-NT-1 

GOTO  5 
C 

C  WRITE  A  LABEL 
C 

300  CALL  ZERO(ID,100) 

ID(l)-"i77777 

WRITE(7,311) 

311  FORMAT('$LABEL  (LESS  THAN  50  CHARS)  >') 

READ(5,310)  (ID(L),L-2,50) 

310  F0RMAT(49A1) 

CALL  PUT<0,ID,100,MSTAT) 

CALL  WAIT(MSTAT) 

CALL  ZER0(ID,100) 

GOTO  5 
END 
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Appendix  J 


List  of  Subroutine  calls  not  included  In  the  subroutine  section  of  PMSM.FLX. 
A  description  of  the  subroutine  function  and  the  arguments  is  included. 
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Adkill 


-  Abort  an  A/D  Conversion  in  Progress 


Adkill  —  Abort  an  A/D  Conversion  In  Progress 

Format: 

CALL  adkill 

Description: 

The  ADKILL  Subroutine  aborts  any  ADSET  that  might  be  in 
progress  at  the  time.  It  is  useful  only  in  those  cases  where 
the  program  needs  to  abort  any  I/O  activity  that  might  be  in  progress 
before  Exiting.  When  it  is  known  that  no  ADSET  is  in  progress,  it 
is  not  necessary  to  Call  this  Subroutine  before  Exiting.  The  "Stat" 
variable  for  the  ADSET  Call  Aborted  is  not  Set  to  One  by  this  Call. 
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AdSet 


-  Digitize  a  Block  of  Data 


AdSet  —  Digitize  a  Block  of  Data 

Format : 

CALL  adset 

(buf fer,gain, start ,nchan,npts, rate, count ,s tat) 

Where: 

buffer  (Integer  Array,  Nchan  by  Npts) 

The  Buffer  into  which  the  the  Digitized  Data  are 
stored. 

gain  (Int  Val ,  0-3)  Programmable  Gain  Select 

0:  Multiply  Input  by  1  (-10  to  10  Volt  Range) 

1:  Multiply  Input  by  2  (-5  to  5  Volt  Range) 

2:  Multiply  Input  by  4  (-2.5  to  2.5  Volt  Range) 

3:  Multiply  Input  by  8  (-1.25  to  1.25  Volt  Range) 

start  (Int  Val,  1-16)  Channel  to  Start  each  Scan  With 

nchan  (Int  Val,  1-16)  Number  of  Channels  to  Sample 

on  each  Scan 

npts  (Int  Val)  The  Number  of  Scans  to  do,  or,  the  number 

of  Points  in  each  Sweep.  (Npnts  x  Nchan  =  Buffer  Size) 

rate  (Int  Val,  0-7)  Basic  Tick  Rate  for  the  Inter  Scan  Clock. 
Rates  Are: 

0:  1  Microsecond  1  Megahertz 

1:  10  Microseconds  100  Kilohertz 

2:  100  Microseconds  10  Kilohertz 

3:  1  Millisecond  1  Kilohertz 

4:  10  Milliseconds  100  Hertz 

5:  100  Milliseconds  10  Hertz 

6:  1  Second  1  Hertz 

7:  10  Seconds  100  Millihertz 

count  (Int  Val)  Number  of  Units  of  "Rate"  to  count  off 
between  each  scan. 

stat  (Int  Var)  This  Variable  keeps  the  user  informed  of  the 
status  of  the  A/D  Conversion.  When  the  ADSET  is 
executed,  this  variable  is  cleared,  (that  is,  assigned 
the  value  Zero).  When  the  last  Scan  of  the  Sweep  is 
finished  and  the  user's  data  array  is  full,  "Stat"  is 
assigned  the  Value  One... 

0:  A/D  Conversion  in  Progress 

1:  A/D  Converter  Sweep  is  Finished 


Description : 

The  ADSET  Subroutine  permits  its  user  to  Digitize  a  Block  of 
Data  composed  of  "Npts"  sets  of  "Nchan"  Samples.  The  Call  works 
as  Follows.  When  an  ADSET  is  Executed,  "Nchan"  Channels  starting 
at  Channel  "Schan"  are  Digitized  right  away. 

A  Set  of  Samples  across  several  channels  digitized  at  (virtually)  the 
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same  time  such  as  this  will  be  refered  to  herein  as  a  "Scan". 

Then  an  interval 

determined  by  "Rate"  and  "Count"  is  clocked  off  (by  the  A/D  System's 
built  in  clock),  and  another  Scan  is  performed. 

This  process  continues  until  "Npts"  Scans  have  been  performed. 

A  complete  set  of  clock  scheduled  Scans  as  described  above  is 
refered  to  as  a  "Sweep". 

The  Array  into  which  the  Data  will  go  should  be  DIMENSIONed 
to  contain  "Nchan"  by  "Npts"  Elements.  BUFFER(I,J)  will  address 
the  datum  read  from  Channel  "I"  during  the  "J"th  Scan. 

The  Values  returned  for  each  digitized  value  are  12  Bit 
Signed  Integers.  A  value  of  -10  Volts  presented  to  the  A/D  Inputs 
will  result  in  a  value  of  -2048  being  returned  by  Adset  when  that 
channel  is  sampled  and  digitized.  A  value  of  0  Volts  will  cause 
a  Zero  to  be  returned.  A  value  of  +5  Volts  returns  1024  and  a 
value  of  10  Volts  returns  2048.  Other  voltages  vary  the  A/D 
result  proportionally  between  the  two  extreme  value,  of  course. 

The  input  sensitivity  of  the  A/D  System  can  be  altered  with 
the  "Gain"  Parameter.  This  parameter  reduces  the  Voltage  Range 
over  which  the  Digitizer's  Values  are  assigned,  and,  thus,  has 
the  effect  of  increasing  the  resolution  with  which  the  values  in 
the  more  restrictred  range  are  digitized. 

The  "Schan"  Parameter  determines  the  starting  input 
channel  number.  The  "Nchan"  Parameter  indicates  the  number  of 
channels,  starting  with  "Schan"  that  will  be  samples  on  each  Scan. 

"Npts"  determines  the  number  of  Scans  that  will  be  made  by 
the  A/D  System.  Note  that  there  will  be  "Npts"  Samples  for  each 
of  the  "Nchan"  Channels  sampled  in  "Buffer"  when  ADSET  is  finished. 


The  "Rate"  and  "Count"  Parameters,  taken  together,  determine 
the  interval  between  Scans.  The  Rate  Parameter  selects  the  basic 
unit  for  the  A/D's  InterScan  Clock.  For  Example,  a  "Rate"  Parameter 
of  3  and  a  "Count"  Parameter  of  5  will  result  in  sucesslve  scans  being 
separated  by  5  Milliseconds. 

The  "Stat"  Variable  is  Zeroed  by  the  ADSET  Call  Itself. 

As  soon  as  the  digitizer  is  started,  the  ADSET  Call  Returns. 

The  User's  Program  may  go  about  other  business  while  the 
Digitizer  Sweep  is  performed.  When  the  complete  sweep  is  finished, 
the  "Stat"  Variable  is  Set  to  One. 
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Cl init 


-  Enable  the  Hard  Clock  so  that  Software  Timers  can  be  Used 


Clinit  —  Enable  the  Hard  Clock  so  that  Software  Timers  can  be  U^i-d 

Format : 

CALL  clinit 
CALL  clinit(unit) 

CALL  clinit(unit , base, count) 

Where: 

unit  (Int  Val)  Number  of  the  Hard  Clock  Unit  to 
use  with  CLPACK. 

base  (Int  Val)  Basic  Rate  for  the  Programmable 
Clock. 

count  (Int  Val)  Number  of  Units  of  "Base"  to  Count 
Off  per  Each  Interrupt. 

Description: 

The  CLINIT  Subroutine  Starts  Up  the 
Programmable  Clock  that  services  CLPACK's  Software  Timers. 

This  Subroutine  must  be  executed  before  any  CLSET  or  CLSTOP 

Calls  may  be  executed. 

The  Zero  Argument  Form  of  CLINIT  should  be  used 
for  most  applications. 
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Clkill 


-  Turn  Off  CLPACK's  Hardware  Programmable  Clock 


Clkill  —  Turn  Off  CLPACK's  Hardware  Programmable  Clock 

Format: 

CALL  clkill 

Description: 

The  CLKILL  Subroutine  Disables  the  Programmable  Clock 
that  services  CLPACK's  Software  Timers.  This  has  the  effect 
of  turning  off  any  Software  Timers  that  might  have  been  active 
since,  with  no  Hardware  Clock  running,  they  cannot  be  serviced 
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Diinit 


-  Enable  Parallel  Inputs 


Diinit  —  Enable  Parallel  Inputs 

Format : 

CALL  diinit 
CALL  diinit(max) 

Where: 

max  (Int  Val)  Number  of  the  Highest  Bit  for  which  an 
Input  Can  Occur. 

Description: 

The  DIINIT  Subroutine  Enables  the  Parallel  Input 
Interface.  This  Subroutine  must  be  executed  in  order  for  the 
DISET,  DIVAL,  and  DIRESP  Subroutines  to  work. 

The  Optional  "Max"  Parameter  Indicates  the  Highest  Bit  Number 
for  which  inputs  are  expected. 
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Dlkill  -  Turn  off  Digital  Inputs  Entirely 

Dikill  —  Turn  off  Digital  Inputs  Entirely... 

Format : 

CALL  dikill 

Description: 

The  DIKILL  Subroutine  turns  off  Parallel  Inputs, 
thus  effectivly  terminating  all  digital  input  recording 
until  another  DIINIT  is  done. 
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DiResp  -  Report  an  input  for  some  bit  or  other 

DiResp  —  Report  an  input  for  some  bit  or  other 
Format : 

CALL  DiResp(bit ,raode,stat) 

Where: 

bit  (Int  Var)  Place  in  which  the  Bit  Position  of  the 

latest  Digital  Input  "Hit"  is  Stored, 
mode  (Int  or  String  Val)  Single  or  Repeat  Interrupt  Mode: 

0  or  'Single':  Single  Interrupt  Mode. 

1  or  'Repeat':  Repeat  Interrupt  Mode, 
stat  (Int  Var)  Variable  that  is  incremented  whenever 

an  input  occurs.  DiResp  Zeroes  this  Flag  Initially. 

-OR- 

CALL  DiResp(bit ,mode, type,src,dst,stat) 

Where: 

bit  (Int  Var)  Place  in  which  the  bit  number  of  the 

latest  digital  input  "hit"  is  returned. 

mode  (Int  or  String  Value)  Single  or  Repeat  interrupt  mode: 

0  or  'Single':  Single  interrupt  mode. 

1  or  'Repeat':  Repeat  interrupt  mode, 
type  (Int  Value)  Where  to  get  response  time: 

1:  Read  Sofware  Timer  indicated  by  "src". 

2:  Read  Programmable  Clock  indicated  by  "src". 

3:  Copy  variable  given  as  "src"  when  a  response 

occurs. 

src  (Int  Value  or  Variable)  Depends  on  "type"  Above, 

dst  (Int  Variable)  Variable  to  which  the  response  time 

for  the  latest  hit,  as  indicated  by  "mode" 
and  "src"  will  be  copied.  Note  that  this  variable  is 
not  zeroed  or  otherwise  changed  by  the  call  to  DiResp 
unless  an  input  event  occurs.  If  no  input  activity  takes 
place,  "dst"  retains  the  value  it  had  before  the  call  to 
DiResp. 

stat  (Int  Variable)  Variable  that  is  incremented  whenever 
a  digital  input  occurs. 

-OR- 


CALL  direspO 
Description: 

The  DIRESP  subroutine  enables  a  Fortran  program  to 
respond  to  Inputs  for  any  of  the  Parallel  I/O  System's  Bits. 
Whenever  an  Input  Bit  is  Set  (by  the  External  Hardware),  the 
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Bit  Number  on  which  the  Input  occurred  is  copied  into  "Bit"  and 
One  is  added  to  the  value  of  "Flag".  "Flag”  is  Zeroed  by  the 
initial  call  to  DIRESP. 

The  "Mode"  Parameter  spcecifes  whether  DIRESP  will 
stop  reporting  bits  after  the  first  hit,  or  whether  it  will 
continue  to  report  bits  until  it  is  disabled. 

The  "Type",  "Src"  and  "Dst"  Parameters  allow  the 
User  to  record  a  response  time  along  with  the  bit  number. 

The  "Mode"  and  "Src"  Paramters  indicates  where  this  response  time 
should  be  obtained.  When  a  Response  occurs,  this  time 
is  written  to  "Dst". 

If  DIRESP  is  called  with  no  parameters,  the  effect  of 
any  previous  DIRESP  Call  is  canceled. 

DIINIT  must  be  executed  before  DIRESP  is  called  in  order 
for  DIRESP  to  work  effectively.  The  execution  of  a  DIINIT  or 
a  DIKILL,  will  (of  course)  cancel  any  pending  DIRESP, 
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i)  i  s  et 


-  Enable  DigiiaL  InpuL  Reporting  for  a  Given  Bit 


O.i-et  —  Enabit  Digital  Input  R- porting  for  a  Given  Bit. 


'■  ornca  t  : 

CALL  d i  s  e  t  ( b i  t ,  mod  e ,  s  t  a  i_  ) 

UTi  ere: 

bit  (Int  Var,  O-io)  Nuraber  of  the  bit  to  watch. 

mode  (Int  or  Stririg  Val)  Single  or  Repeat  Interrupt  Mode: 

0  or  'Single':  Single  Interrupt  Mode. 

1  or  'Repeat':  Repeat  Interrupt  Mode, 
stat  (Int  Var)  Variable  that  will  be  up'ed  when 

an  input  from  the  Bit  indicated  by  "Bit"  is  detected. 
"Stat"  is  Zero'ed  by  the  Call  to  DISET. 

-OR- 

CALL  diset ( bi t ,mode, type , src ,dst , flag ) 

Wher  e : 

hit  (Int  Val,  0-15)  Bit  Position  to  Watch. 

mode  (Int  or  String  Val)  Single  or  Repeat  Interrupt  Mode: 

0  or  'Single':  Single  Interrupt  Mode. 

1  cr  'Repeat':  Repeat  Interrupt  Mode. 

type  (Int  Val)  Where  to  get  Response  Time: 

I:  Read  Sofware  Timer  indicated  by  "Src". 

2:  Read  Programmable  Clock  indicated  by  "Src". 

3:  Copy  Var  given  as  "Src"  when  a  response 

occurs. 

src  (Int  Var)  Depends  on  "Mode"  Above, 

dst  (Int  Var)  Variable  to  which  the  Response  time 

for  the  latest  hit,  as  indicated  by  "Mode" 
and  "Src"  will  be  copied. 

stat  (Inc  Var)  Variable  that  is  Incremented  whenever 

an  Input  on  the  indicated  bit  comes  in.  This  variable 
is  Cleared  by  the  DISET  Call  Itself. 

Description: 

The  DISET  Subroutine  enables  a  L'ser's  Fortran  Program  to 
respond  to  activity  on  a  given  Parallel  Interface  Input  Bit. 

When  the  DISET  Subroutine  is  executed,  subsequent  inputs  for  the 
"Watched"  Bit  will  cause  the  indicated  "Stat"  Variable  to  be 
incremented.  Since  each  Input  will  cause  "Stat"  to  increase, 

"Stat"  will  be  a  count  of  the  number  of  Inputs  that  have  occurred 
for  Che  Watched  Bit. 

The  "Mode"  Parameter  spcecifes  whether  DIRESP  will 
stop  reporting  bits  after  the  first  hit,  or  whether  it  will 
f'ont-inue  to  report  bits  until  it  is  disabled. 

The  "Type",  "Src"  and  "Dst"  Parameters  allow  the 
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User  to  record  a  response  time  along  with  the  bit  number. 

The  "Mode"  and  "Src"  Paramters  indicates  where  this  response  time 
should  be  obtained.  When  a  Response  occurs,  this  time 
is  written  to  "Dst". 

Inputs  on  the  Watched  Bit  will  be  processed  until  either 
a  DICLR  for  the  Watched  Bit,  or  a  DIINIT  or  DIKILL  is  executed. 

The  DIINIT  Subroutine  must  be  executed  before  DISET  can 
be  used. 
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Doclr 


-  Clear  a  Specified  Output  Bit 


Doclr  —  Clear  a  Specified  Output  Bit 
Format : 

CALL  doclr(bit) 

Where: 

bit  (Int  Val,  0-15)  Bit  to  Turn  Off 
Description: 

The  DOCLR  Subroutine  Turns  Off  the  Digital  Output  Bit  Designated 
by  the  "Bit"  Parameter.  Only  the  designated  bit  is  affected. 

This  Subroutine  clears  only  the  Bit  Indicated. 
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Doinit 


-  Clear  the  Digital  Output  Register 


Dolnlt  —  Clear  the  Digital  Output  Register 

Format: 

CALL  dolnlt 

Description: 

The  DOINIT  Subroutine  Clears  the  Parallel  Interface's 
Digital  Output  Register.  All  16  Digital  Output  Bits  are  Turned  Off. 
It  Is  not  necessary  to  execute  a  DOINIT  to  use  the  other  Digital 
Output  Subroutines,  but  it  is  often  a  good  idea  to  make  sure  that 
all  the  Digital  Output  Bits  are  Off  when  starting  a  program  that  uses 
Digital  Outputs.  A  Call  to  DOINIT  is  equivalent  to  a  Call  to  DOKILL 
or  a  DOVAL(O). 
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Dokill  -  Clear  the  Digital  Output  Register 

Format : 

CALL  dokill 


Description: 

The  DOKILL  Subroutine  Clears  All  the  Digital  Output 
Bits.  It  is  identical  to  the  DOINIT  Subroutine.  The  Two  Names 
are  provided  for  compatibility  with  other  Digital  I/O  Subroutine 
Packages.  It  is  not  necessary  to  perform  a  DOKILL  before  Exiting 
or  Aborting  a  Program  that  has  set  Output  Bits,  but,  depending 
on  the  nature  of  the  equipment  connected  to  the  Digital  Output 
System,  it  is  often  a  good  idea. 


Doset 


-  Turn  on  a  Specified  Digital  Output  Bit 


Doset  —  Turn  on  a  Specified  Digital  Output  Bit 
Format : 

CALL  do8et(bit) 

Where: 

bit  (Int  Val,  0-15)  Output  Bit  to  Set 
Description: 

The  DOSET  Subroutine  turns  on  the  Bit  designated  by  the 
"Bit"  Parameter.  Bits  are  Numbered  from  Right  to  Left  in  the 
Digital  Output  Register  from  0  to  15.  This  is  the  conventional 
way  that  bit  positions  are  numbered  on  the  PDPll.  The  designated 
Output  Bit  remains  set  until  a  DOCLR  is  executed  for  it,  or  until 
a  DOINIT  or  a  DOKILL  is  Called.  The  DOINIT  Subroutine  need  not 
be  executed  to  use  DOSET.  DOINIT  and  DOKILL  are  equivalent  to 
DOCLRs  for  bits  0  through  15. 
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GqDisp 


-  Display  a  waveform 


GqDisp  —  Display  a  waveform 
Format : 

CALL  GqDisp ( loc, chan, data ,nchan,npt s, scale) 

VThere: 

loc  (Int  Value,  1-12)  This  parameter  gives  the  waveform  location 
Waveforms  are  numbered  with  waveform  1  at  the  upper  left, 
waveform  3  at  the  lower  left,  waveform  4  at  the  upper  right 
and  waveform  6  at  the  lower  right.  Waveforms  7  through  12 
overlay  waveforms  1  through  6. 

chan  (Int  Value)  "Channel"  from  the  data  array  which  is  to  be 
displayed.  This  value  should  be  in  the  range  from  1  to 
nchan. 

data  (Int  Array,  nchan  x  npts)  Array  from  which  the  data  we 
will  display  are  to  be  taken. 

nchan  (Int  Value)  Number  of  channels  in  the  data  n-’-ay 

npts  (Int  Value)  Number  of  points  in  the  data  array. 

scale  (Int  Value)  Waveform  scale  factor.  Eight  is  subtracted  from 
this  value,  and  the  result  is  used  as  follows:  If  >  0, 

each  point  in  the  waveform  is  multiplied  by  the  result. 

If  <  0,  each  point  in  the  waveform  is  divided  by  the 
absolute  value  of  the  result.  If  0,  the  data  are  not 

scaled. 

Description: 

The  GqDisp  subroutine  displays  a  waveform  in  an  indicated  display 
frame  on  the  display  created  using  Gqinit.  Previously  written  data 
are  erased  first.  This  subroutine  is  compatible  with  GTDISP,  except 
that  all  six  waveforms  may  be  overlayed.  Hence,  "loc"  values  7  through 
9  will  overlay  the  waveforms  in  the  left  hand  column  of  the  display. 
Instead  of  in  the  right  hand  column  as  with  GTPACK.  To  overlay  the 
right  hand  column,  use  "loc"  values  10  through  12. 
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Gqlrit 


-  Initialize  the  QRGB-Graph  Waveform  Display 


Gqinit  —  Initialize  the  QRGB-Graph  Waveform  Display 
Format: 

CALL  Gqlnlt(flag) 

Where: 

flag  (Int  Value)  Dummy  argument  for  compatibility  with 
GtPack. 


Description: 

The  Gqinit  subroutine  sets  up  the  Matrox  QRGB-Graph  waveform 
display.  The  "flag"  argument  does  nothing,  it  is  provided  for 
compatibility  with  GTPACK.  The  waveform  display  is  always  redrawn 
by  Gqinit. 
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GqLabl  -  Set  a  QRGB-Graph  display  label 

GqLabl  —  Set  a  QRGB-Graph  display  label 
Format : 

CALL  GqLabl(loc, string) 

fcfhere: 

loc  (Int  Value,  1-12)  Label  l.ccation, 
string  (String)  New  label  string. 

Description: 

The  GqLabl  subroutine  resets  a  Matrox  QRGE-Graph  waveform 
display  label.  To  erase  a  label,  an  explicit  null 
string  (either  a  ",  or  a  null  terminated  array)  mutt  be  given. 

No  zero  argument  form  is  allowed.  Note  that  GTPACK  allowed  a 
zero  argument  form  to  clear  a  label.  Labels  are  located  just  below 
the  waveform  frame  corresponding  the  the  "loc"  value  given.  Labels 
for  "loc"  values  1-6  are  written  above  those  for  values  7-12. 
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GqTtl 


-  Set  the  QRGB-Graph  display  title  line 


GqTtl  —  Set  the  QRGB-Graph  display  title  line 
Format: 

CALL  GqTtl(strlng) 

Where: 

string  (String)  New  title  string. 

Description: 

The  GqTtl  subroutine  resets  the  Matrox  QRGB-Graph  waveform 
display  title  line.  To  erase  the  title  line,  an  explicit  null 
string  (either  a  or  a  null  terminated  array)  must  be  given. 
No  zero  argument  form  is  allowed.  Note  that  GTPACK  allowed  a 
zero  argument  form  to  clear  the  title  line. 
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McChar 


-  Write  a  Character  on  the  Matrox 


McChar  —  Write  a  Character  on  the  Matrox 
Format : 

CALL  mcchar(c) 

CALL  mcchar(x,y,c) 

Where: 

c  (Byte  Val)  Character  to  write... 

x,y  (Int  Vais)  New  X  and  Y  Positions... 


Description: 

The  MCCHAR  Subroutine  writes  a  character  on  the  screen  of  the 
current  size  and  color  at  the  given  screen  position.  If  no 
position  is  given,  the  current  cursor  position  is  assumed. 


142 


McDisp 


-  Turn  Matrox  Display  On  or  Off 


McDisp  —  Turn  Matrox  Display  On  or  Off 
Format : 

CALL  mcdlsp(val) 

Where; 

val  (Int  Val)  0;  Turn  Off  the  Matrox,  1:  Turn  On  the  Matrox 
Description: 

The  MCDISP  Subroutine  turns  the  Matrox  On  or  Off. 
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McEras 


-  Make  an  Area  the  Background  Color 


McEras  —  Make  an  Area  the  Background  Color 

Format: 

CALL  me eras 
CALL  mceras(x,y) 

CALL  mceras(a,b,x,y) 

Vfhere: 

a,b  (Int  Vais)  Coordinates  of  Erase  Area  Upper  Left  Hand  Corner. 

x,y  (Int  Vais)  Coordinates  of  Erase  Area  Lower  Right  Hand  Corner. 

Description: 

This  MCERAS  Subroutine  sets  an  area  of  the  Matrox  Screen  to  the 
current  background  color.  If  no  coordinates  are  given,  a  hardware  flood  is 
performed  using  the  current  background  color.  If  two  coordinates  are 
given, 

the  rectangle  determined  by  the  current  position  on  the  upper  left,  and 
the  given  point  on  the  lower  right  will  be  colored  using  the  current 
background  color.  If  four  parameters  are  given,  the  first  two  will 
designate  the  upper  left  hand  corner  of  the  area  to  erase,  and  the 

next  two  will  designate  the  lower  right  hand  corner  of  the  area. 
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Mclnit 


-  Initialize  the  Matrox 


Mclnit  —  Initialize  the  Matrox 

Format : 

CALL  mclnit 
CALL  mcinit(flag) 

Where: 

flag  (Int  Val)  0:  Don't  Erase.  1:  Erase. 

Description: 

The  MCINIT  Subroutine  initializes  the  Matrox.  It  should  be 
called  before  any  other  MCPACK  Subroutines  are  called. 

If  a  "Flag"  value  is  given,  the  screen  is  erased  or  not  erased 
as  indicated.  If  no  value  is  given,  the  screen  is  not  erased. 
Leaving  the  screen  unerased  by  MCINIT  allows  programs  to  generate 
Matrox  images  and  chain  to  other  programs  without  losing  them. 

MCINIT  sets  defaults  as  follows: 


Foreground:  256,  or  White. 

Background:  0,  or  Black. 

Spacing:  Nonproportional. 

Position:  (0,0) 

Origin:  (0,0) 

Scroll:  0 

Video:  On 

Reverse:  Off 

If  an  Initial  erase  is  not  done,  some  of  these  default  values 
will  not  be  established  until  a  subsequent  MCPACK  Subroutine  accesses 
the  Matrox. 
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McSync  -  Wait  for  Blanking  to  Start 

McSync  —  Wait  for  Blanking  to  Start 

Format : 

CALL  mcsync 

Description: 

The  MCSYNC  Subroutine  waits  until  the  Matrox's  Vertical  Blanking 
Bit  goes  from  Off  to  On. 
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Mctext  -  Write  a  Null  Terminated  String  on  the  Matrox 

Mctext  —  Write  a  Null  Terminated  String  on  the  Matrox 
Format ; 

CALL  met ext (string) 

CALL  mctext(x,y, string) 

Where: 

x,y  (Int  Vais)  Place  to  put  the  cursor  at  first 
string  (Byte  Array,  or  Whatnot)  A  null  terminated 
buffer  filled  with  Ascii  Characters. 

Description: 

The  MCTEXT  writes  a  Null  Terminated  String  of  Characters 
out  on  the  Matrox  in  the  current  Color  and  Size.  If  no  position 
is  specified,  the  current  cursor  position  is  used. 
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Pcinit 


-  Initialize  a  Programmable  Clock 


Pcinit  —  Initialize  a  Programmable  Clock 
Format: 

CALL  pcinit(clock) 

Where: 

clock  (Int  Val,  1-6)  Clock  to  Set  Up. 

-  OR  - 

l=lpclnl( ) 

Where: 

i  (Int  Result)  Clock  Allocated  (1-6,  or  0  if  None  are 

available. ) 

Description: 

The  PCINIT  Subroutine  initializes  the  indicated  Programmable 
Clock.  This  Subroutine  must  be  called  before  any  given  clock 
can  be  used.  "Clock"  specifies  which  clock  should  be  initialized. 
The  given  clock  is  marked  as  busy  until  a  PCKILL  is  performed  for 
it.  Attempts  to  Initialize  it  in  the  mean  time  will  fail. 


The  IPCINI  Function  returns  the  number  of  an  available 
Programmable  Clock,  and  initializes  that  clock.  If  no  clocks 
are  available.  Zero  is  Returned. 
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Pckill 


-  Turn  Off  and  Free  One  or  All  Clocks 


Pckill  —  Turn  Off  and  Free  One  or  All  Clocks 

Format: 

CALL  pckill 
CALL  pcklll(clock) 

Where: 

clock  (Int  Val,  1-6)  Clock  to  Kill  and  Free 
Description: 

The  PCKILL  Subroutine  Turns  Off  and  Frees  One  or  all 
of  the  Programmable  Clocks.  If  "Clock"  is  specified,  only 
the  Indicated  Clock  is  Killed.  If  "Clock"  is  not  specified, 
all  the  Clocks  are  Killed. 
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Pcset 


-  Turn  On  a  Programmable  Clock 


Pcset  —  Turn  On  a  Programmable  Clock 
Format: 

CALL  pcset(clock, mode, rate, count) 

CALL  pcset (clock, mode, rate, count ,stat ) 

Where: 

clock  (Int  Val,  1-6)  Clock  to  Use. 
mode  (Int  Val)  Single  or  Repeat  Mode. 

0  or  'Single':  Single  Sweep  Mode 
1  or  'Repeat':  Repeat  Sweep  Mode 
rate  (Int  Val,  0-7)  Basic  Clock  Rate: 

0:  1  Microsecond  1  Megahertz 

1:  10  Microseconds  100  Kilohertz 

2:  100  Microseconds  10  Kilohertz 

3:  1  Millisecond  1  Kilohertz 

4:  10  Milliseconds  100  Hertz 

5:  100  Milliseconds  10  Hertz 

6:  1  Second  1  Hertz 

7:  10  Seconds  100  Mllllhertz 

count  (Int  Val,  0-32767)  Number  of  Units  of  "Rate"  per 
Interrupt.  (0  means  32768  Units.) 
stat  (Int  Var)  Status  Variable.  This  Is  Zeroed  when 
PCSET  Is  Called.  One  Is  Added  to  "Stat"  every 
time  a  Clock  Sweep  Completes. 

Description: 

The  PCSET  Subroutine  starts  one  of  Pearl  II's  programmable 
clocks.  Clocks  may  be  run  either  in  "Single  Sweep  Mode",  in 
which  case  one  interval  of  the  specified  rate  is  timed,  or  in 
"Repeat  Sweep  Mode",  in  which  the  clock  times  intervals  indefinitely 
until  stepped. 

The  "Clock"  Parameter  specifies  which  of  Pearl  II's  programmable 
clocks  should  be  used. 

The  "Mode"  Parameter  indicates  whether  the  Clock  should 
run  in  Single  or  Repeat  Sweep  Mode. 

The  "Rate"  Parameter  specifies  the  Basic  Rate  at  which 
the  Clock  will  run.  "Count"  units  at  the  given  rate  will 
be  timed  for  each  sweep. 

The  "Stat"  Variable  is  zeroed  by  the  call  to  PCSET. 

One  is  added  to  this  Variable  every  time  a  Clock  Sweep  is 
completed. 

If  "Stat"  is  not  given,  the  clock  is  run  without  interrupts. 
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If  a  timer  is  running  when  this  subroutine  is  called,  it 
is  turned  off  first. 


151 


Walt 


-  Wait  for  a  Variable  to  Become  Non-Zero 


Wait  —  Wait  for  a  Variable  to  Become  Non-Zero 

Format: 

CALL  wait(vl) 

CALL  walt(vl ,v2 , . . . ,vn) 

Where: 

vl...vn  (Int  Vara)  The  Variables  to  Test  for  Non-Zero  Values 
Description: 

The  WAIT  Subroutine  Tests  the  Values  of  All  the  Variables  Passed 
to  it.  As  soon  as  ANY  of  these  variables  becomes  something  other 
than  Zero,  the  Subroutine  Returns  to  the  User's  Program.  The 
typical  application  for  this  Subroutine  is  testing  the  "Stat" 

Variables  that  keep  the  user  informed  of  the  progress  of  various 
other  LABPAK  Subroutine  initiated  activity. 

Note: 


The  WAIT  Subroutine  is  Implemented  using  the  PDPll's  WAIT 
Instruction.  This  instruction  causes  the  Central  Processor 
to  relinquish  control  of  the  Bus  to  those  DMA  Devices  that 
might  be  competing  for  it  until  an  Interrupt  Request  has  been 
serviced.  The  Variable  list  is  first  tested  when  this  Subroutine 
is  called,  and  if  the  contents  of  all  the  variables  are  zero, 
a  WAIT  instruction  is  executed.  When  this  instruction 
falls  through,  (meaning  that  an  Interrupt  service  routine 
that  might  have  changed  the  value  of  one  of  the  variables  was 
executed)  the  tests  are  performed  again.  When  One  of  the  Variables 
becomes  Non-Zero,  the  Subroutine  falls  through. 
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Zero 


-  Array  Zeroing  Subroutine 


Zero  —  Array  Zeroing  Subroutine 
Format : 

CALL  zero(buf .size) 

Where: 

Buf  (Array,  "Size"  Bytes  in  Size)  The  Array 
to  Fill  with  Zeroes* 

Size  (Int  Val)  The  Size,  In  Bytes,  of  the  "Buf"  Array. 
Description; 

The  ZERO  Subroutine  Sets  every  element  of  the  indicated 
Array  to  Zero.  It  does  this  in  a  manner  more  efficient  than  the 
corresponding  DO  Loop  Fortran  Code  could  do  it,  and  is  more 
concise.  The  Size  in  Bytes  of  an  Array  Variable  may  be  calculated 
by  multiplying  the  sizes  of  the  Array's  Dimensions  together  and 
multiplying  this  result  by  the  size  of  an  individual  element  of  the 
Array. 

Example: 

The  following  short  Fortran  Program  illustrates  the  use 
of  the  ZERO  Subroutine: 

PROGRAM  zero 
C 

C  Declare  Some  Arrays  Here... 

C 

INTEGER  dog (2, 4, 6)  !2*4*6*2  -  96  Bytes 

REAL  mouse( 100,2)  1100*2*4  »  800  Bytes 

C 

C  Zero  these  and  Exit... 

C 

CALL  zero(dog ,2*4*6*2)  !Dog  is  all  Zeroes  Now... 

CALL  zero(cat,  100*2*4)  .'Cat  is  Zeroed,  too... 

STOP  'Dog  and  Cat  Zeroed...'  !Say  what  this  did... 

END 
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Appendix  K 


F.PEAK  Program  to  perform  single  trial  data  analysis  on  ERPs  collected  by 
PMSM.  Also  generates  are  average  waveforms  for  each  experimental  condition. 
This  program  creates  unadjusted  and  latency  adjusted  waveforms,  and  single 
trial  estimates  of  P300  latency  and  amplitude.  Four  different  algorithms 
are  included;  Vector  cross  correlation,  Pz  Cross  correlation.  Vector  peak 
picking,  and  Pz  peak  picking. 
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c****** ********************************** ********* 


c  * 

C  F.PEAK  * 
C  * 
C  PROGRAMMER;  DAVID  STRAYER  * 
C  UNIVERSITY  OF  ILLINOIS  * 
C  * 


(]************************************************* 

C 

C  IDS  USED  AND  GENERATED  BY  THIS  PROGRAM 
C 

C  ID(1)»SUBJECT 

C  ID(2)=SESSION 

C  ID(3)=BLOCK 

C  ID(4)»EEG  CHANNEL  (EOG,  FZ,  CZ,  PZ,  STM) 

C  ID(5)-DELAY 

C  ID(6)=SOA 

C  ID(7)-MEMSIZ 

C  ID(8)=CMVM 

C  ID(9)=REACTION  TIME 

C  ID(10)=RESPONSE  CLASSIFICATION  (0=N0  RESPONSE,  1=HIT,  2=CR, 
C  3=FA,4=MISS) 

C 

C  ID'S  INCLUDED  IN  THE  OUTPUT  OF  THIS  PROGRAM 
C 

C  SINGLE  TRIAL  ESTIMATES  OF  P300  TO  SI  (300  TO  600  MSEC) 

C 

C  ID(11)»LATENCY  OF  CROSS  CORR  P300  (COMPOSITE) 

C  ID(12)»AMPLITUDE  OF  CROSS  CORR  P300  (COMPOSITE) 

C  ID (13) “LATENCY  OF  CROSS  CORR  P300  (PZ  CHANNEL) 

C  ID(14)=AMPLITUDE  OF  CROSS  CORR  P300  (PZ  CHANNEL) 

C  ID(15)»LATENCY  OF  PEAKPICKED  P300  (COMPOSITE) 

C  ID(16)“AMPLITUDE  OF  PEAKPICKED  P300  (COMPOSITE) 

C  ID(17)“LATENCY  OF  PEAK  PICKED  P300  (PZ  CHANNEL) 

C  IDO 8) “AMPLITUDE  OF  PEAK  PICKED  P300  (PZ  CHANNEL) 

C 

C  SINGLE  TRIAL  ESTIMATES  OF  P300  TO  S2  (800  TO  1700  MSEC) 

C 

C  1D(19)“LATENCY  OF  CROSS  CORR  P300  (COMPOSITE) 

C  ID(20)“AMPLITUDE  OF  CROSS  CORR  P300  (COMPOSITE) 

C  ID(21)“LATENCY  OF  CROSS  CORR  P300  (PZ  CHANNEL) 

C  ID(22)“AMPLITUDE  OF  CROSS  CORR  P300  (PZ  CHANNEL) 

C  ID(23)“LATENCY  OF  PEAKPICKED  P300  (COMPOSITE) 

C  ID( 24) “AMPLITUDE  OF  PEAKPICKED  P300  (COMPOSITE) 

C  ID(25)“LATENCY  OF  PEAK  PICKED  P300  (PZ  CHANNEL) 

C  ID(26)“AMPLITUDE  OF  PEAK  PICKED  P300  (PZ  CHANNEL) 

C 

REAL  ID( 26 ),DATA( 210,5 ),YDATA( 210,5), 
+FZ(210),CZ(210),PZ(210),P3(210),BASE(3),TEMP(50), 
+NTRIAL(24,5),AVEID(26,24,5),AVE(210,3,24,5),BAVE(  ,;:4,5) 

C 
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INTEGER  FNAME( 6 ) , LOW( 2 ) ,HIGH( 2 ) , BIN , IDS( 108 ) 

C 

CALL  0PEN(3,ISTAT) 

DEFINE  WINDOW  ONSET  AND  OFFSET 

DATA  LOW(1),HIGH(1)/40,70/  !  FOR  SI 

DATA  L0W(2) ,HIGH(2)/90,180/  !  FOR  S2 

ASSIGN  FILES 

CALL  ASARGd,  20, "SUB  ".FNAME) 

CALL  ASARG(2,  30, "DATA  ".FNAME) 

CALL  ASARG(3,-40,"PK-0UT  ".FNAME) 

CALL  ASARG(4,-41,"ADJ-OUT  ".FNAME) 

DEFINE  TEMPLET  (POSITIVE  SEGMENT  OF  0.5  HZ  SINE  WAVE) 

FOR  I- 1,50 

TEMP(I)»~C0S(  1*6. 2829/50) 

END  FOR 

SET  UP  FILTER  PARAMETERS  SEE  RUTCHKIN  AND  GLASER  (1979) 
TT«10.0  I  SAMPLING  INTERVAL 

NPTS-210  !  NUMBER  OF  POINTS  IN  WAVEFORM  TO  FILTER 

LGL-6  !  LENGTH  OF  LOW  PASS  FILTER 

NUML=2  !  NUMBER  OF  ITERATIONS  OF  LOW  PASS  FILTER 

LGH=0  !  LENGTH  OF  HIGH  PASS  FILTER 

NUMH-0  !  NUMBER  OF  ITERATIONS  OF  HIGH  PASS  FILTER 

NSTART- 1+ ( LGL*NUML ) + ( LGH*  NUMH ) 

NS  TOP=NPTS- ( LGL*  NUML+LGH*  NUMH ) 

CALL  FILTN( l.NPTS.NSTART, NSTOP.l, LGL, NUML, LGH, NUMH) 

DESCRIBE  FILTER  PARAMETERS 

F0L=1000.0/((2*LGL+1)*TT) 

IF(NUML.EQ.l)  F3DB».44*F0L 
IF(NUML.EQ.2)  F3DB=.31*F0L 
WRITE(3,*)  F3DB,'H2  (-3dB)  low  pass  filter' 

READ  IN  A  TRIAL 

LOOP 

FOR  J-1,5 

RE  AD  (30, 100, END-50)  (ID(K)  ,K-1 ,10) ,(DATA(I,J) ,1-1 ,210) 
END  FOR 

FIND  THIS  TRIAL  IN  THE  SUB  ID  FILE 
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1  READ(20,101)  IDS 

IF(ID(l).NE.IDS(l).OR.ID(2).NE.IDS(2).OR.ID(3).NE.IDS(3))  GOTO  1 
DELAY=IDS(4) 

DETERMINE  NUMBER  OF  CORRECT  RUNNING  MEMORY  TRIALS 
ICORR=0 

IF(DELAY.GT.2)  THEN 
FOR  I»3,ABS(DELAY) 

IF( IDS( 1*5+7) .EQ. IDS( 1*5+8))  IC0RR=IC0RR+1 
END  FOR 
END  IF 

IFdCORR.GE. ((DELAY-2. 0)*0. 50))  THEN  !  ABOVE  CHANCE  ACCURACY 
FOR  EP0CH=1,2 

COMPUTE  BASELINE  (AVERAGE  OF  100  MSEC  EPOCH) 

51  EPOCH  BASELINE:  -100  MSEC  TO  0  MSEC 

52  EPOCH  BASELINE:  500  MSEC  TO  600  MSEC 

FOR  J=l,3 
BASE(J)=0. 

FOR  K» 1+ (( EPOCH-1 )* 59 ),10+(( EPOCH-1) *59) 
BASE(J)»BASE(J)+(DATA(K,(J+1))/10.0) 

END  FOR 
END  FOR 

SUBTRACT  BASELINE 
FOR  IPT-1,210 

FZ(IPT)=DATA( IPT,2)-BASE(1) 

CZ(  IPT)=-DATA(  IPT , 3)-BASE( 2 ) 

PZ(IPT)“DATA( IPT,4)-BASE(3) 

END  FOR 

FILTER  THE  EEG  DATA 

CALL  FILTR(FZ,FZ) 

CALL  FILTR(CZ,CZ) 

CALL  FILTR(PZ,PZ) 

COMPUTE  THE  P300  VECTOR  FILTER 

CALL  VECT(210,300.,15.,FZ,CZ,PZ,P3) 

COMPUTE  CROSS  CORRELATION  ON  P300  VECTOR  (COMPOSITE) 

CALL  CROS S ( 2 1 0 , LOW( EPOCH ) , HIGH( EPOCH ) , TEMP , P3 , CORP 3 , L AT , AMP ) 
ID(3+(8*EPOCH))-LAT 
ID(4+(8*EP0CH))-AMP 
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IF(C0RP3.LE..30)  ID(3+(8*EPOCH))=0.0 
COMPUTE  CORES  CORRELATION  ON  PZ  CHANNEL 

CALL  CROSS(210,LOW(EPOCH),HIGH(EPOCH) , TEMP, PZ,CORPZ,LAT, AMP) 

ID(5+(8*EPOCH))=LAT 

ID(6+(8*EPOCH))=AMP 

IF(CORPZ.LE. .30)  1D(5+(8*EPOCH))=0.0 

PEAK  PICK  ON  P300  VECTOR  FILTER  OUTPUT 

CALL  PKPICK(P3,210,L0W(EP0CH),HIGH(EP0CH),LAT,AMP) 

ID(7+(8*EPOCH))=LAT 

ID(8+(8*EP0CH))=AMP 

IF( ID(7+(8*EPOCH) ) .LE. (LOW(EPOCH)*10-100) .OR.ID(7+(8*EPOCH) ) .GE. 
+( HIGH( EPOCH)* 10-100 ) )  1D( 7+( 8*EPOCH) )=0 .0 

PEAK  PICK  ON  PZ  CHANNEL 

CALL  PKP ICK( PZ , 2 1 0 , L0W( EPOCH ) , HIGH( EPOCH) , LAT , AMP ) 

ID(9+(8*EP0CH))=LAT 

ID( 10+( 8*EP0CH) )=AMP 

IF( ID(9+( 8*EP0CH) ) .LE. (LOW(EPOCH)*10-100) .OR. ID(9+(8*EP0CH) ) .GE. 
+(HIGH( EPOCH)* 10-100))  ID(9+(8*EP0CH) )»0.0 

END  FOR 

WRITE  OUT  ID  FILE  WITH  PEAK  ESTIMATES 

WRITE(40,102)  ID 

NOW  LATENCY  ADJUST  EACH  TRIAL 

FOR  EST=1,5 

CONVERT  P3  LATENCY  TO  NUMBER  OF  DATA  POINTS 

IF(EST.NE.5)  THEN 
IF(ID((EST*2)+17).Eq.O)  GO  TO  40 
NPOINT=(ID((EST*2)+17)+100)/10 
ELSE 

NPOINT=120  !  EST  5  IS  THE  UNADJUSTED  AVERAGE 
END  IF 

WHICH  WAY  SHOULD  WE  SHIFT  THE  WAVEFORM? 

IF( NPOINT.lt.  120)  ISHIFT— 1 
IF(NPOINT.EQ.120)  ISHIFT-0 
IF(NP0INT.GT.120)  ISHIFT-1 

NOW  SHIFT  THE  WAVEFORM 
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C 

IF(ISHIFT.EQ.O)  THEN 

FOR  K-2,4 

FOR  J^1,210 

YDATA( J ,K ) =DATA( J , K) 

END  FOR 
END  FOR 
GO  TO  150 
C 

ELSE  IF(ISHIFT.EQ.-l)  THEN 

SHIFT  THIS  WAVEFORM  TO  THE  RIGHT - > 

NDIFF=120-NP0INT 
JMOVE=210-NDIFF 
FOR  K=2,4 
XBASE*DATA(1,K) 

FOR  J=l,JMOVE 

YDATA(  J+NDIFF  ,K)<=DATA(  J  ,K) 

END  FOR 

NOW  SET  THE  TRUNCATED  REGION  TO  THE  FIRST  DATA  POINT 

FOR  J=1,NDIFF 
YDATA(J,K) -XBASE 
END  FOR 
END  FOR 
GO  TO  150 

ELSE  IF(ISHIFT.EQ.l)  THEN 

SHIFT  THIS  WAVEFORM  TO  THE  LEFT  < - - - 

NDIFF-NP0INT-120 
JMOVE-210-NDIFF 
FOR  K=2,4 
XBASE-DATA(210,K) 

FOR  J»l,JMOVE 

YDATA( J ,K)»DATA( J+NDIFF ,K) 

END  FOR 
C 

C  NOW  SET  THE  TRUNCATED  REGION  TO  THE  LAST  DATA  POINT 
C 

ITEMP-210-NDIFF+1 
FOR  J-ITEMP,210 
YDATA(J,K)-XBASE 
END  FOR 
END  FOR 
GO  TO  150 
END  IF 
C 
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CONTINUE 


150 


IDENTIFY  WHICH  BIN  TRIAL  SHOULD  BE  ASSIGNED  TO 

BIN=0 
C 

IF(ID(10).EQ.l)  THEN  !  TARGETS 
C 

IF(ID(8).EQ.l)  THEN  !  CM  CONDITIONS 
C 

IF(ID(5).EQ.0.AND.ID(7).EQ.2)  BIN=1 
IF(ID(5).EQ.4.AND.ID(7).EQ.2)  BIN=2 
IF(ID(5).EQ.15.AND.ID(7).EQ.2)  BIN=3 
IF(ID(5).EQ.0.AND.ID(7).EQ.4)  BIN=4 
IF(ID(5).EQ.4.AND.ID(7).EQ.4)  BIN=5 
IF(ID(5).EQ.15.AND.ID(7).EQ.4)  BTN=6 
C 

ELSE  IF(ID(8).EQ.2)  THEN  !  VM  CONDITIONS 
C 

IF(ID(5).EQ.0.AND.ID(7).EQ.2)  BIN=7 
IF(ID(5).EQ.4.AND.ID(7).EQ.2)  BIN=8 
IF(ID(5).EQ.15.AND.ID(7).EQ.2)  BIN=9 
IF(ID(5).EQ.0.AND.ID(7).EQ.4)  BIN*10 
IF(ID(5).EQ.4.AND.ID(7).EQ.4)  BIN=11 
IF(ID(5).EQ.15.AND.ID(7).EQ.4)  BIN=12 
END  IF 
C 

ELSE  IF(ID(10).EQ.2)  THEN  !  DISTRACTORS 
C 

IF(ID(8).EQ,1)  THEN  !  CM  CONDITIONS 
C 

IF(ID(5).EQ.O.AND.ID(7).EQ.2)  BIN=13 
IF(ID(5).EQ.4.AND.ID(7).EQ.2)  BIN=14 
IF(ID(5).EQ.15.AND.ID(7).EQ.2)  BIN=15 
IF(ID(5).EQ.0.AND.ID(7).Eq.4)  BIN=16 
IF(ID(5).EQ.4.AND.ID(7).EQ.4)  BIN=17 
IF(ID(5).EQ. 15.AND.ID(7).EQ.4)  BIN=18 
C 

ELSE  IF(ID(8).EQ.2)  THEN  !  VM  CONDITIONS 
C 

IF(ID(5).EQ.0.AND.ID(7).EQ.2)  BIN=19 
IF(ID(5).EQ.4.AND.ID(7).EQ.2)  BIN=20 
IF(ID(5).EQ.15.AND.ID(7).EQ.2)  BIN=21 
IF(ID(5).EQ.0.AND.ID(7).Eq.4)  BIN=22 
IF(ID(5).Eq.4.AND.ID(7).Eq.4)  BIN-23 
IF(ID(5).Eq.l5.AND.ID(7).Eq.4)  BIN=24 
C 

END  IF 
C 

END  IF 

IF(BIN.Eq.O)  GO  TO  40 
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c  SUM  THE  POINTS 

NTRIAL(BIN,EST)=NTRIAL(BIN,EST)+1 
FOR  KID=1,26 

AVEID(K1D,BIN,EST)=ID(KID) 

END  FOR 

FOR  LCHAN=2,4 

FOR  LPT=1,210 

AVE(LPT,LCHAN-1,BIN,EST)= 

+ AVE ( LPT , LCHAN- 1 , BIN , E  ST ) + YD AT A( LPT , LCHAN ) 

END  FOR 
END  FOR 
C 

END  FOR 
END  IF 

40  END  LOOP 

50  CONTINUE 

C 

C  CALCULATE  AVERAGES  FOR  EACH  CONDITION 

C 

FOR  EST=1,5 
DO  300  MBIN=1,24 

IF(NTRIAL(MBIN,EST).EQ.O)  GO  TO  300 
FOR  MCHAN=1,3 
FOR  MPT=1,210 

AVE(MPT,MCHAN,MBIN,EST)=AVE(MPT,MCHAN,MBIN,EST)/NTRIAL(MBIN,EST) 

IF(MPT.LE. 10)  BAVE(MCHAN,MBIN ,EST)=BAVE(MCHAN,MBIN ,EST) 
*+(AVE(MPT,MCHAN,MBIN,EST)/10) 

END  FOR 
END  FOR 
300  CONTINUE 

C  SUBTRACT  BASELINE  FROM  AVERAGE  AND  OUTPUT  AVERAGES 

FOR  JBIN=1,24 
FOR  JCHAN=1,3 
FOR  JPT=1,210 

AVE(JPT,JCHAN,JBIN,EST)=AVE(JPT,JCHAN,JBIN,EST) 

*-BAVE( JCHAN , JBIN ,EST) 

END  FOR  .  „  XX 

WRITE(41 , 100)  EST, REAL (JBIN) ,REAL( JCHAN) , REAL (NTRIAL( JBIN ,EST)) , 

+( AVEID(K, JBIN, EST),K=1 ,26) ,(AVE(N, JCHAN, JBIN.EST) ,N=1 ,210) 

END  FOR 
END  FOR 
END  FOR 
STOP 

100  FORMAT(22F6.0) 

101  FORMAT(16I5) 

102  FORMATC 11F6.0) 

END 
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r* 

SUBROUTINE  VECT(NPT,OR,POL 

,FZ,CZ,PZ,P3) 

L 

c 

THIS  SUBROUTINE  PERFORMS 

A 

VECTOR  FILTER 

c 

FOR  ANY  GIVEN 

POLARITY  AND 

ORIENTATION. 

c 

r* 

THE  DATA  MUST 

COME  FROM 

THREE  SCALP  ELECTRODES. 

Kj 

c 

ARGUMENTS ; 

c 

NPT 

if  OF  POINTS  (MAX=700) 

c 

OR 

ORIENTATION  ANGLE 

c 

FOR  P300  USE  OR=300  OR  330 

c 

POL 

POLARITY  ANGLE 

c 

FOR  P300  USE  P0L=15  OR  30 

c 

FZ 

(210) 

DATA  ARRAYS.  IT  IMPLIES: 

c 

CZ 

(210) 

c 

PZ 

(210) 

c 

EL1=FZ,  EL2=CZ,  EL3=PZ 

c 

n 

P3 

(210) 

VECTOR  FILTERED  DATA 

L 

r 

DIMENSION  FZ(NPT) ,CZ(NPT) ,PZ(NPT) ,P3(NPT) 

VLOAD=COS( POL*6 . 2829/360 

.) 

XLOAD=COS(OR*6 

.2829/360. 

)*VLOAD 

YLOAD=SIN(OR*6 

.2829/360. 

)*VLOAD 

ZL0AD=SIN( P0L*6 .2829/360 

.) 

c 

FOR  K=1,NPT 
C  ROTATE  AXES 


Z=(FZ(K)+CZ(K)+PZ(K))*. 57735 
Y=(2.*CZ(K)-(FZ(K)+PZ(K)))*.40825 
X=(PZ(K)-FZ(K))*. 70711 
C  VECTOR  FILTER 

P3(K)=X*XLOAD+Y*YLOAD+Z*ZLOAD 
END  FOR 
RETURN 
END 
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SUBROUTINE  CROSS 
C 

SUBROUTINE  CROS S( NPTS , ILOW , IHIGH , TEMP , DATA , CORMAX , LAT , BMAX ) 
DIMENSION  DATA( NPTS ) , TEMP( 50 ) ,WIND( 50 ) 

LBEG=ILOW 

LEND=IHIGH 

IF(LBEG.LT.25)LBEG=25 
IF(LEND.GT.(NPTS-25))  LEND=NPTS-25 
CORMAX=0. 

FOR  LAG=LBEG,LEND 
IBEG=LAG-25 
IEND=LAG+24 
FOR  I=IBEG,IEND 

WIND( I-IBEG+1 )=DATA( I) 

END  FOR 

CALL  CORR( WIND, TEMP, 50, R,B) 

I F ( R. GE . CORMAX ) THEN 
CORMAX=R 
BMAX=B 

LAT=LAG* 10-100 
END  IF 
END  FOR 
RETURN 
END 
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SUBROUTINE  CORK 

SUBROUTINE  CORR(X,Y,NPTS,R,B) 
DIMENSION  X(NPTS) ,Y(NPTS) 

R=0. 

IF(NPTS.LT. 2) RETURN 
XSUM=0 . 

YSUM=0. 

X2SUM=0. 

Y2SUM>=0. 

XYSUM=0. 

FOR  I-1,NPTS 
XSUM=X(I)+XSUM 
YSUM=*Y(I)+YSUM 
X2SUM=X(I)*X(I)+X2SUM 
Y2SUM=Y( I)*Y( I)+Y2SUM 
XYSUM=X( I)*Y( I)+XYSUM 
END  FOR 

XMEAN-XSUM/NPTS 

YMEAN»YSUM/NPTS 

XVAR=X2 SUM/NPTS-XMEAN*XMEAN 

YV AR»  Y2  SUM/ NPTS-YME AN* YME AN 

C0V=XYSUM/ NPTS-XMEAN*  YME AN 

IF( XV AR. LE . 0 .. OR. YVAR. LE . 0 .) RETURN 

XSD=SQRT(XVAR) 

YSD=SQRT(YVAR) 

R=COV/(XSD*YSD) 

B=R*XSD/YSD 

RETURN 

END 
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SUBROUTINE  PKPICK  (VECTOR, NPTS,ILOW,IHIGH,LAT, AMP) 
DIMENSION  VECTOR(NPTS) 

PEAK=VECTOR( ILOW) 

LAT=ILOW 
FOR  J=ILOW,IHIGH 
IF(VECTOR(J).LE.PEAK)  GO  TO  10 
PEAK=VECTOR(J) 

LAT=J 
END  FOR 

LAT=LAT* 10-100 

AMP=PEAK 

RETURN 

END 
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C  FILTER  SUBROUTINES  TAKEN  FROM  CPL  LIBRARY  LIB*F.TSLIB 
C 

SUBROUTINE  FILTR(WAV1 ,WAV2) ,FILTN(NEW,LTIME,NSTRT, 

1 NSTOP , NTYPE , LGL , NUML , LGH , NUMH ) 

C  LOW  PASS  AND  HIGH  PASS  FILTER  IN  CASCADE 

C  WAV1=INPUT;  WAV2=OUTPUT 

C  PROGRAMMER  -  D.S.RUCHKIN 

C  VERSION  -  5/25/75 

DIMENSION  WAVl(LTIME) ,WAV2(LTIME) ,TWAV(2001) 

INTEGER  HSTRT.HSTT.HSTOP.HSTP.HRAN 
C 

IF(NTYPE.EQ.O)  RETURN 
IF(LGH.EQ.O.AND.LGL.EQ.O)  RETURN 
IF(LGH.NE.O)  GO  TO  20 
C 

C  ONLY  LOW  PASS 

C 

CALL  BOX(LSTRT,LSTOP, LGL, NUML, DIVL,WAV1 ,WAV2,LTIME) 
CALL  EDGE(NSTT1 ,NSTP1 ,LTIME,WAV1 ,WAV2) 

RETURN 

C 

C  HIGH  OR  BAND  PASS 

C 

20  CALL  BOX(HSTRT,HSTOP, LGH, NUMH, DIVH,WAVI,TWAV,LTIME) 

DO  25  K=HSTRT,HSTOP 
25  WAV2(K)-WAV1(K)-TWAV(K) 

IF(LGL.NE.O)  GO  TO  50 
C 

C  ONLY  HIGH  PASS 

C 

40  CALL  EDGE(NST1,NSTP1,LTIME,WAV1,WAV2) 

RETURN 

C 

C  BAND  PASS 

C 

50  CALL  BOX(LSTRT,LSTOP,LGL,NUML,DIVL,WAV2,WAV2,LTIME) 

CALL  EDGE(NSTT1 ,NSTP1 ,LTIME,WAV1 ,WAV2) 

RETURN 

C 

C  INITIALIZATION  ENTRY 

C 

ENTRY  FILTN 

C  IF  NEW-0;  RETAIN  OLD  PARAMETERS 

C  IF  NEW-1;  READ  NEW  PARAMETERS 

IF(NEW.EQ.O)  GO  TO  280 
NSTTl-NSTRT-1 
NSTPl-NSTOP+1 
IF(LGL.EQ.O)  GO  TO  120 
C 

C  LOW  PASS  PARAMETERS 

C 
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DIVL=2*LGL+1 

DIVL-l.O/DIVL 

LRAN=(NUML-1)*LGL 

LSTRT=NSTRT-LRAN 

LSTT=LSTRT-LGL 

LSTOP=NSTOP+LRAN 

LSTP-LSTOP+LGL 

IF(LSTT.LT.l.OR.LSTP.GT.LTIME)  GO  TO  250 
C 

C  HIGH  PASS  PARAMETERS 

C 

120  IF(LGH.EQ.O)  GO  TO  200 
DIVH=2*LGH+1 
DIVH=1.0/DIVH 
IF(LGL.EQ.O)  GO  TO  140 
C 

C  BAND  PASS-HIGH  PARAMETERS 

C 

HRAN-(NUMH-1)*LGH 
HSTRT=LSTT-HRAN 
HSTOP=LSTP+HRAN 
GO  TO  145 
C 

C  ONLY  HIGH  PASS 

C 

140  HRAN»(NUMH-1)*LGH 
HSTRT=NSTRT-HRAN 
HSTOP=NSTOP+HRAN 
145  HSTT=HSTRT-LGH 
HSTP=HSTOP+LGH 

IF(HSTT.LT. l.OR.HSTP.GT.LTIME)  GO  TO  250 
C 

C  OUTPUT  AND  RETURN 

C 

200  WRITE (6, 501)  NTYPE ,NSTRT,NSTOP,LGL,NUML,LGH,NUMH 

IF(LGL.NE.O)  WRITE(6,502)  LSTRT,LSTO?,DIVL 
IF(LGH.NE.O)  WRITE(6,503)  HSTRT.HSTOP.DIVH 
RETURN 
C 

C  ERROR 

C 

250  NTYPE-0 

WRITE(6,504) 

GO  TO  200 
C 

C  SAME  PARAMETERS 

C 

280  IF( NTYPE. EQ.O)  GO  TO  200 
WRITE(6,505) 

RETURN 


C 
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C  FORMATS 
C 


500 

FORMAT (12 

,614) 

501 

FORMAT (' 

FILTR, 5/25/75 ,NTYPE='I2/ 

,5X,'NSTRT='I4 

1 

'  NSTOP=',I4,'  LGL-'.U/ 

NUML=',I4,' 

2 

NUMH=',I4) 

502 

FORMAT( 5X 

,'LSTRT=',I4,'  LSTOP=',I4,' 

DIVL=',F8.5) 

503 

FORMAT(5X 

,'HSTRT=',I4,'  HSTOP»',I4,' 

DIVH='F8.5) 

504 

FORMAT (' 

FILTER  ERROR') 

505 

FORMAT (' 
END 

NO  FILTER',/) 
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SUBROUTINE  BOX(NSTRT,NSTOP,LG,NUMH,DIV ,WAVE1 ,TWAV1 .NPOINT) 


C 
C 
C 
C 
C 

C 

C  ITERATIVE  BOX-CAR  FILTER 

C  WAVE 1= INPUT;  TWAVE=OUTPUT 

C 

DIMENSION  WAVEl(NPOINT),TWAVl(NPOINT),WAVE3(2O01),TWAVE(2O01) 
C 

C  TRANSMIT  INPUT  TO  WORK  AREA 

C 

NCL=NSTRT-LG 
NCU=NSTOP+LG 
DO  4  J=NCL,NCU 

4  WAVE3(J)=WAVE1(J) 

NUMIT=1 

MSTRT=NSTRT 

MSTOP=NSTOP 

C 

C  BEGIN  ITERATION  LOOP 

C 

10  DO  50  J=MSTRT,MSTOP 
SUM=0.0 
NCL=J-LG 
NCU=J+LG 
DO  40  K=NCL,NCU 
40  SUM=SUM+WAVE3(K) 

50  TWAVE(J)=SUM*DIV 

IF(NUMIT.GE.NUMH)  GO  TO  100 

NUMIT=NUMIT+1 

DO  5  J=MSTRT,MSTOP 

5  WAVE3(J)=TWAVE(J) 

MSTRT=MSTRT+LG 

MSTOP=MSTOP-LG 

GO  TO  10 
C 

C  MOVE  OUTPUT  AND  RETURN 

C 

100  DO  110  K=NSTRT,NSTOP 

110  TWAV1(K)=TWAVE(K) 

C  MODIFICATION  BY  E.  HEFFLEY  9/13/77 

C  MODIFIED  MOVING  AVERAGE  DONE  AT  BEGINNING 

C  AND  END  OF  EPOCH.  THIS  IS  FOR  CASES  WHEN  ITERATIONS-1 . 

C  IT  AVERAGES  ACROSS  LESS  THAN  2*LG+1  POINTS  BY  REDISTRIBUTING 
C  WEIGHTS  ACROSS  POINTS  AT  BEGINNING  (OR  END)  OF  EPOCH 
C  SUCH  THAT  POINTS  CLOSER  TO  THE  BEGINNING  ARE  GIVEN  HIGHER 
C  WEIGHTS  BECAUSE  PRESUMABLY  THEY  ARE  BETTER  PREDICTORS  OF 
C  THE  PRECEDING  POINTS  (WHICH  WERE  NOT  RECORDED  BUT  SHOULD  HAVE 
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C  BEEN  TO  ALLOW  FILTERING  OF  FIRST  LG  POINTS  IN  THE  EPOCH). 
BASE=DIV 
ISTART=NSTOP+l 
DO  200  IPT=I START, NPOINT 
NLEFT=(LG+l)+(NPOINT-IPT) 

NREMN=( 2*LG+1 )-NLEFT 

AREA=NREMN*BASE 

DX=NLEFT 

DY= ( 2.0* AREA) /DX 
SLOPE=DY/(DX-1.0) 

IFIRST=IPT-LG 

SUM=0.0 

WTOT=0.0 

DO  160  IPTR=IFIRST, NPOINT 
WEIGHT=( IPTR-IFIRST)*SLOPE 
WTOT*WTOT+WE IGHT 

SUM=SUM+WAVE1(IPTR)*(WEIGHT+BASE) 

160  CONTINUE 

TWAV1(IPT)=SUM 

C  WRITE(6,1200)  ISTART, NLEFT.NREMN, AREA, DX.DY, SLOPE ,WTOT 
1200  FORMATC'  START*', 14,'  LEFT*', 14,'  NREMN*' , 14 ,/ , 

1  '  AREA=',F8.4,'  DX=',F8.4,'  DY=',F8.4,'  SLOPE*' ,F8.4 , 

2  /,'  WTOT=',F8.4) 

200  CONTINUE 

1END»NSTRT-1 

DO  300  IPT*IEND,1,-1 

NLEFT*(LG+1)+(IPT-1) 

NREMN*(2*LG+1)-NLEFT 

AREA=NREMN*BASE 

DX=NLEFT 

DY=(2.0*AREA)/DX 

SLOPE=DY/(DX-1.0) 

ILAST*IPT+LG 

SUM=0.0 

WTOT*0.0 

DO  260  IPTR=ILAST,1 ,-l 
WEIGHT*( ILAST-IPTR)*SLOPE 
WTOT*WTOT+WEIGHT 

SUM*  SUM+WAVE 1 ( IPTR) * ( WE IGHT+BASE ) 

260  CONTINUE 

TWAV1(IPT)=SUM 

C  WRITE(6,1220)  IEND,NLEFT, NREMN, AREA, DX,DY, SLOPE, WTOT 
1220  FORMATC'  END-', 14,'  LEFT-', 14,'  NREMN*' , 14 ,/ , 

1  '  AREA*',F8.4,'  DX=',F8.4,'  DY*',F8.4,'  SLOPE*' ,F8.4 , 

2  /,'  WTOT*',F8.4) 

300  CONTINUE 

RETURN 

END 
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c 

c 

c 

c 

SUBROUTINE  EDGE(NSTT,NSTP,LTIME,WAV1 ,WAV2) 
C  EDGE  FIXING  ROUTINE  CALLED  BY  FILTR 

DIMENSION  WAVl(LTIME) ,WAV2(LTIME) 

C 

DO  10  K=1,NSTT 
10  WAV2(K)=WAV1(K) 

DO  20  K=NSTP,LTIME 
20  WAV2(K)=WAV1(K) 

RETURN 

END 
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Appendix  L 


TSD.FOR  A  program  to  compute  various  signal  detection  parameters  given  Hit 
and  False  Alarm  Rates. 


172 


c 

C  F.TSD  —  A  PROGRAM  TO  CALCULATE  SEVERAL  MEASURES  OF 
C  SENSITIVITY  AND  BIAS  FROM  SIGNAL  DETECTION 

C  THEORY  GIVEN  THE  FOLLOWING: 

C 

C  PROGRAMMER:  DAVID  STRAYER 

C  UNIVERSITY  OF  ILLINOIS 

C 
C 

C  NUMBER  OF  HITS 

C  NUMB’''?  OF  FALSE  ALARMS 

C 

C  THE  INPUT  CAN  BE  RAW  SCORES,  OR  PROPORTIONS  RELATIVE  TO 
C  CATEGORY  (E.G.,  p(HIT)+p(MISS)=l .0  AND  p(FA)+p(CR)=l .0) 

C 


INTEGER*!  IANS 

REAL  HIT,  FA,  CR,  MISS,  APRIME,  BPMPM,  AG,  E,  PS,  BETA,LNBETA, 
+DPRIME,  ZI,  Z2,  TEST!,  TEST2,  CO,  Cl,  C2,  Dl,  D2,  D3 ,  K, 

+T1,  T2,  T3,  T4,  BI ,  BPRIME 
C 

1  WRITE(3,*)  'Enter  HIT  rate' 

READ(0,*)  HIT 
IF(HIT.GT.I.O)  THEN 
WRITE(3,*)  'Enter  MISS  Rate' 

READ(0,*)  MISS 
HIT=HIT/(HIT+MISS) 

END  IF 

MISS=I. 0-HIT 
C 

WRITE(3,*)  'Enter  False  Alarm  (FA)  rate' 

READ(0,*)  FA 
IF(FA.GT. I.O)  THEN 

WRITE(3,*)  'Enter  Correct  Rejection  (CR)  rate' 

READ(0,*)  CR 
FA=FA/(CR+FA) 

END  IF 
CR=1.0-FA 

WRITE(3,*)  '  Response 


WRITE(3,*)  ' 
WRITE(3,*)  ' 
WRITE(3,98)HIT, 
FORMAT('  Target 
WRITE(3,*)  ' 

URTTKY'?  *■>  ' 

Yes 

No 

98 

MISS 

presented 

:  ',F3.2,' 

:  HIT  : 

:  ',F3 

MISS  : 

99 

WRITE(3,99)  FA,  CR 

FORMAT('  Distractor  presented 
WRITE(3,*)  ' 

:  ',F3.2,' 

:  FA  : 

:  ',F3 

CR  : 

WRITE(3,*)  '  : . . 

WRITE(3,*)  '  ' 

WRITE(3,*)  '  Are  ALL  parameters  correct?  (Y/N)' 
READ(0,2)  IANS 
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non 


2  FORMAT(Al) 

IFdANS.NE.'Y'. AND. IANS. NE.'y')  GOTO  1 
C 

C  FUDGE  FACTORS  TO  PREVENT  DIVISION  BY  ZERO 

C 

IF(HIT.EQ.l.O)  THEN 
HIT=.9999 
MISS=1-HIT 
END  IF 
C 

IF(HIT.EQ.O.O)  THEN 
HIT=.0001 
MISS=1-HIT 
END  IF 
C 

IF(FA.EQ.l.O)  THEN 
FA=.9999 
CR=1-FA 
END  IF 
C 

IF(FA.EQ.O.O)  THEN 
FA=.0001 
CR=1-FA 
END  IF 
C 
C 

C  FORMULA  FOR  CALCULATING  A'  (A  NONPARAMETRIC  MEASURE  OF  EFFICIENCY) 
C 

C  FROM  WICKENS(1984) ,  (P.  502)  SEE  ALSO  POLLACK  AND  NORMAN  (1964) 

C  RECOMMENDED  BY  PARASURAMAN  AND  DAVIES  AS  THE  BEST  MEASURE  OF 

c  SENSITIVITY 

A'= 1-. 25* (( FA/HIT) +( 1-HIT) /(1-FA)) 

APRIME=l-.25*((FA/HIT)+( l-HIT)/( 1-FA)) 

C 

C  FORMULA  FOR  CALCULATING  (ANOTHER  NONPARAMETRIC  AREA  MEASURE) 

C 

C  FROM  CRAIG,  1979  SEE  ALSO  GREEN  AND  SEWTS  (1966) 

C  Ag=  (HlT+(l-FA))/2 

C 

AG=(HIT+( l-FA))/2.0 
C 

C  FORMULA  FOR  CALCULATING  E  (YET  ANOTHER  NONPARAMETRIC  AREA  MEASURE) 
C 

C  FROM  CRAIG  (1979)  SEE  ALSO  MCCORNACK  (1961) 

C 

C  E=l-(  (1-HIT)  /  (  (1-HIT)  *  PS  +  (FA  *  PS)  -  FA)) 

C  WHERE  PS  IS  SIGNAL  PROBABILITY  (BATCH  QUALITY) 

C  (THIS  PROGRAM  ASSUMES  PS-.5) 

C 
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n  o 


PS=.5 

E=  1  -  (  (1-HIT)  /  (  (1-HIT*PS)  +(FA*PS-FA)  )) 

C 

C  FORMULA  FOR  CALCULATING  d' ,  B  AND  In  B 
C 

C  FROM  GARDNER  AND  BOICE  (1986)... 

C 

C 

C0=2. 515517 
C1=0. 802853 
C2=0. 010328 
Dl=l. 432788 
D2=0. 189269 
D3=0.001308 

K=SQRT( 1/(2*3.141592654)) 

Tl=SQRT(LOG( 1/(MISS*MISS) ) ) 

T2=SQRT(LOG( 1/(HIT*HIT) ) ) 

T3=SQRT(LOG(l/(CR*CR))) 

T4=SQRT( LOG( 1 / ( FA*FA) ) ) 

C 

TEST1=0 

IF(HIT.GT. .50)  TEST1=1 

Z1=(T1-(C0+C1*T1+C2*T1*T1)/(1+D1*T1+D2*T1*T1+D3*T1*T1*T1))*TEST1- 

&(T2-(C0+C1*T2+C2*T2*T2)/(1+D1*T2+D2*T2*T2+D3*T2*T2*T2))*(1-TEST1) 

C 

TEST2-0 

IF  (FA.GT..50)  TEST2=1 

Z2=(T3-(C0+C1*T3+C2*T3*T3)/(1+D1*T3+D2*T3*T3+D3*T3*T3*T3))*TEST2- 
&(T4-(C0+C1*T4+C2*T4*T4)/(1+D1*T4+D2*T4*T4+D3*T4*T4*T4))*( 1-TEST2) 
C 

DPRIME=Z1-Z2 

C 

BETA=(K*EXP(-(Zl*Zl)/2))/(K*EXP(-(Z2*Z2)/2)) 

C 

LNBETA=LOG(BETA) 

YET  ANOTHER  MEASURE  OF  RESPONSE  BIAS  B'  FROM  McNICHOL  (1972) 
c  SEE  PARASURAMAN  AND  DAVIES  CHAPTER  3,  P46 
C 

C  B'=BETA/SD(s) 

C  SD(s)  IS  TAKEN  AS  THE  INVERSE  OF  THE  ROC  SLOPE.  THE  SLOPE  OF  THE 

C  ROC  IS  EQUAL  TO  THE  LIKELIHOOD  RATIO  AT  CRITERION  (BETA) 

C  (THIS  SEEMS  TO  BE  AN  UNSTABLE  ESTIMATE) 

C 

BPRIME-BETA/( 1/BETA) 

C 

C 

C  CALCULATE  BETA" 

C 

C  USING  THE  FOLLOWING  FORMULA  FROM  JOHN  POLICH 

C  NOTE:  B"  GETS  WEIRD  IF  P(HIT)  OR  P(FA)  »  1.0 


175 


B"  =  (HIT(  1-HIT) -FA(  1-FA)  )  /(HIT(  1-HIT)+FA(  1-FA)  ) 


BPMPM=(HIT*( 1 H1T)-FA*( 1 .—FA) ) /(HIT*( 1 .-HIT)+FA*( 1 .-FA) ) 

WRITE(3,*)  '  Sensitivity  measures: 

+'Dprime  Aprime  Ag  E 

WRITE(3,100)  DPRIME, APRIME, AG, E 
WRITE(3,*)  '  ' 

WRITE(3,*)  '  Response  Bias  Measures:  , 

+'Beta  LnBeta  Bprime  BPMPM 

WRITE(3,100)  BETA,  LNBETA,  BPRIME,  BPMPM 
WRITE(3,*)  '  ' 

STOP 

F0RMAT(25X,4(F10.6,4X)) 

END 
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Written  instructions  given  to  the  subject  at  the  onset  of  the  experiment, 
which  describe  the  experimental  procedure. 
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The  experiment  which  you  are  about  to  participate  in  has  been  designed  to 
study  aspects  of  the  way  we  use  memory.  It  is  extremely  important  that  you 
follow  the  instructions  that  you  are  about  to  read  very  closely.  Even  if 
it  seems  that  you  could  perform  the  task  better  by  adopting  a  strategy 
different  from  the  one  you  will  be  instructed  to  use,  please  do  only  as  you 
are  instructed.  Please  read  the  instructions  carefully,  and  if  you  have  any 
questions,  please  ask  them  before  beginning  the  experiment. 

During  the  experiment  you  will  work  at  an  IBM  PC.  Stimuli  will  be  displayed 
on  the  computer's  screen,  and  you  will  make  "yes"  and  "no"  responses  by 
pressing  keys  designated  by  the  experimenter.  There  will  be  many  trials, 
each  of  which  will  proceed  as  follows: 

1)  You  will  see  on  the  screen  two  o»-  four  words.  They  will  remain  on  the 
screen  for  three  seconds.  While  the  words  are  on  the  screen,  commit 
them  to  memory  by  repeating  them  over  and  over  to  yourself. 

2)  The  words  will  disappear.  When  they  do,  a  series  of  numbers  may  appear 
on  the  screen.  (On  some  trials  this  part  will  not  happen,  and  we  will 
skip  to  step  3.)  Watch  the  first  two  numbers  carefully,  but  do  nothing. 
When  the  third  number  appears,  compare  it  to  the  first  number  you  saw. 

If  it  is  the  same,  respond  "yes";  if  it  is  not  the  same  as  the  first 
number  press  the  "no"  key.  When  the  fourth  number  appears,  compare  it 
to  the  second  number,  responding  "yes"  or  "no"  as  you  did  for  the  third 
number.  Continue  comparing  every  number  that  appears  on  the  screen  to 
the  number  two  before  it.  I.E.;  The  fifth  number  will  be  compared  to  the 
third  number,  the  sixth  number  will  be  compared  to  the  fourth  number, 
etc..  Respond  "yes"  or  "no"  as  quickly  and  accurately  as  you  can.  This 
part  of  the  experiment  is  difficult.  Do  not  be  discouraged  if  it  takes 
you  a  little  while  to  get  good  at  it.  Do  not  think  about  the  words  you 
saw  at  the  start  of  the  trial  while  you  are  doing  this  task. 

3)  An  asterisk  will  appear  on  the  screen.  This  signals  you  that  it  is  time 

to  remember  the  words  that  you  saw  at  the  start  of  the  trial. 

4)  A  word  will  appear  on  the  screen.  Respond  by  pressing  the  "yes"  key  if 

if  is  one  of  the  words  you  saw  at  the  start  of  the  trial.  Respond  "no" 
if  it  is  not.  Respond  as  quickly  as  you  can  while  maintaining  accuracy. 
If  you  do  not  know  if  the  word  was  presented  at  the  beginning  of  this 
trial,  think  for  a  little  while,  and  then  guess  if  you  still  cannot 
remember.  Try  to  always  make  a  response. 

5)  There  will  be  a  short  pause,  and  then  a  new  trial  will  begin. 

Things  happen  quickly  in  this  experiment  so  please  be  alert.  There  will  be 

a  break  midway  through  each  session.  The  experimenter  will  have  to  use 
your  computer  to  arrange  things  for  the  second  half.  This  is  your 
opportunity  to  rest.  Begin  again  when  you  feel  you  are  ready. 

Before  you  leave,  make  sure  that  you  and  the  experimenter  agree  upon  when 
your  next  session  is  scheduled  to  take  place. 
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Thank  you  for  your  participation! 
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Appendix  N 


Single  subject  scalp  distribution  ERPs  collected  using  program  PMSM.FLX. 

The  solid  line  represents  the  Fz  electrode,  dash  Cz,  and  dotted  Pz.  Ten 
single  subject  plots  are  included.  Note  that  subjects  8  and  9  were  rejected 
due  to  high  error  rates.  Subject  9  also  has  a  large  amount  of  alpha 
activity.  The  baseline  subtracted  from  each  waveform  is  from  -100  msec 
pre-Sl  to  SI.  The  number  of  trials  included  in  each  average  is  displayed  in 
the  upper  left-hand  portion  of  each  plot.  A  description  of  the  ERPS  is 
included  in  the  results  section  of  the  report. 
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